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Abstract 



There has been much recent research into possible methods of polarizing an antipro- 
ton beam, the most promising being spin filtering, the theoretical understanding of 
which is currently incomplete. The method of polarization buildup by spin filtering 
requires many of the beam particles to remain within the beam after repeated in- 
teraction with an internal target in a storage ring. Hence small scattering angles, 
where we show that electromagnetic effects dominate hadronic effects, are impor- 
tant. All spin-averaged and spin-dependent electromagnetic cross-sections and spin 
observables for elastic spin 1/2 - spin 1/2 scattering, for both point-like particles 
and non-point-like particles with internal structure defined by electromagnetic form 
factors, are derived to first order in QED. Particular attention is paid to spin trans- 
fer and depolarization cross-sections in antiproton-proton, antiproton-electron and 
positron-electron scattering, in the low 1 1 | region of momentum transfer. A thor- 
ough mathematical treatment of spin filtering is then presented, identifying the key 
physical processes involved and highlighting the dynamical properties of the phys- 
ical system. We present and solve sets of differential equations which describe the 
buildup of polarization by spin filtering in many different scenarios of interest. The 
advantages of using a lepton target are outlined, and finally a proposal to polarize 
antiprotons by spin filtering off an opposing polarized electron beam is investigated. 



iii 



Summary 



Immense efforts, both theoretical and experimental, have been afforded to gaining a 
better understanding of the spin structure of the nucleon since the startling results 
from the EMC experiment at CERN in 1988 that the intrinsic valence quarks con- 
tribute only a small fraction of the proton's spin. Yet to this day almost nothing is 
known about the transversity distribution of quarks in the nucleon, the last remaining 
leading twist piece of the QCD description of the partonic structure of the nucleon in 
the collinear limit. A high intensity polarized antiproton beam would be required to 
best analyze the transversity distribution function, via Drell-Yan lepton pair produc- 
tion in the scattering of polarized antiprotons off polarized protons. Unfortunately 
no high intensity polarized antiproton beam has been achieved to date. 

Hence there has been much recent research into possible methods of polarizing an 
antiproton beam, instigated by the recent proposal of the VAX (Polarized Antipro- 
ton experiments) Collaboration at GSI, Darmstadt. The most promising method 
under consideration is spin filtering, the theoretical understanding of which is cur- 
rently incomplete. The method of polarization buildup by spin filtering requires 
many of the beam particles to remain within the beam after repeated interaction 
with an internal target in a storage ring. Hence small scattering angles, where we 
show that electromagnetic effects dominate hadronic effects, are important. The 
theoretical background to this effort is investigated in this thesis. 

We derive fully relativistic expressions for all spin-averaged and spin-dependent 
electromagnetic cross-sections and spin observables for elastic spin 1/2 - spin 1/2 
scattering, for both point-like particles and non-point-like particles with internal 
structure defined by electromagnetic form factors, to first order in QED. Particular 
attention is paid to spin transfer and depolarization cross-sections in antiproton- 
proton, antiproton-electron and positron-electron scattering, in the low 1 1 | region 
of momentum transfer. Of the spin-averaged formula derived we highlight that a 
generalization of the Rosenbluth formula is presented in a new compact Lorentz 
invariant form. It is a two-fold generalization in that the masses of both particles 
are included and both particles are taken to have internal structure determined by 
electromagnetic form factors. While these results are eventually applied to spin 
filtering later in the thesis they are not limited to this application. The complete set 
of spin 1/2 - spin 1/2 helicity amplitudes and spin observables should prove useful 
to many other areas in particle physics. 
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The complete set of spin - spin 1 electromagnetic helicity amplitudes are also 
presented to first order in QED. These are useful in describing the spin-dependent 
scattering of deuterons off carbon nuclei for example. 

A thorough mathematical treatment of spin filtering is then presented, identifying 
the two key physical processes involved: (a) selective scattering out of the ring and 
(b) selective spin flip while remaining in the ring. The dynamical properties of the 
physical system under investigation are highlighted. Sets of differential equations 
are presented and solved which describe the buildup of polarization by spin filtering 
in many different scenarios of interest. These scenarios are: 1) spin filtering of a 
stored beam, 2) spin filtering while the beam is being accumulated, i.e. unpolarized 
particles are continuously being fed into the beam at a constant rate, 3) unpolarized 
particles are continuously being fed into the beam at a linearly increasing rate, i.e. 
the particle input rate is ramped up, 4) the input rate is equal to the rate at which 
particles are being lost due to scattering beyond the ring acceptance angle, the beam 
intensity remaining constant, 5) increasing the initial polarization of a stored beam 
by spin filtering, 6) the input of particles into the beam is stopped after a certain 
amount of time, but spin filtering continues. 

The depolarization of a polarized beam on interaction with an unpolarized target 
or beam, as in the important case of electron cooling, is also investigated and shown 
to be negligible. 

We show that there are advantages of using a lepton target instead of an atomic 
gas target for spin filtering, principal amongst them that antiprotons will not anni- 
hilate with the target as they do with the protons in the atomic targets, leading to a 
loss of beam intensity. After showing that electrons in an atomic target are not mas- 
sive enough to scatter antiprotons beyond the acceptance angle of any storage ring 
we propose using an opposing polarized electron beam, of momentum large enough 
to provide scattering of antiprotons beyond ring acceptance, as a possible method to 
polarize antiprotons by spin filtering. This is presented as a practical application of 
the theoretical work presented throughout the thesis. The areal density of the polar- 
ized electron beam is identified as the key parameter limiting the rate of antiproton 
polarization buildup in this proposal. 
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Chapter 1 
Introduction 



The great Danish physicist Niels Bohr, it is 
said, had a good-luck horseshoe hanging in 
his office. "You don't believe in that non- 
sense, do you?" a visitor once asked, to 
which Bohr replied, "Of course not, but they 
say it brings you good luck whether you be- 
lieve in it or not. " 




Figure 1-1: The great physicists Wolfgang Pauli (left) and Niels Bohr musing over the spin 
of a spinning top toy, trying to gain insight into the nature of spin in particle and nuclear 
physics. 
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1.1 Spin and polarization 



Spin is a fundamental property of elementary particles. It was introduced theoreti- 
cally by Wolfgang Pauli in the 1920's to explain how two electrons can exist in the 
ground state of an atom, while not violating his famous exclusion principle, for which 
he was awarded the 1945 Nobel Prize. This principle states that two particles sat- 
isfying Fermi-Dirac statistics (later called fermions and defined by their half-integer 
spin) cannot exist in the same state at the same time. Thus the two electrons in 
the inner shell of an atom must somehow be different, Pauli hypothesized that they 
have some differentiating characteristic called spin. He theorized one electron to be 
in a 'spin up' state and the other electron to be in 'spin down' state, thus they do 
not violate the exclusion principle. Uhlenbeck and Goudschmidt [1] also introduced 
the concept of spin around the same time as Pauli's work. Pauli's theory of spin was 
non-relativistic, Paul Dirac developed the relativistic theory of spin in 1928 with the 
famous Dirac equation of a relativistic electron [2J. 

Spin was discovered experimentally by the famous Stern-Gerlach experiment in 
1922 [3]. Since then it has been an integral part of the Quantum Field Theories that 
describe particle interactions. 

A beam of spin 1/2 particles will have the spins of each of the particles in either 
the 'spin up' or 'spin down' state. For a beam of spin 1/2 particles the polarization 
is defined as 

N+ _ N _ 

iV + + JV_ ' 1 1 

where N + and iV_ are the number of particles in the 'spin up' and 'spin down' states 
respectively. 

An unpolarized beam (V = 0) is one where the spin states are randomly dis- 
tributed, thus for a beam with a large number of spin 1/2 particles, half of the 
particles will be in the 'spin up' state and half in the 'spin down' state, as seen in 
Figure 11-21 A polarized beam is one where more of the particles are in one spin 
state than the otheiQ. For example, a 100 % polarized beam (V — 1) has all of the 
particles in one of the spin states. 

Originally particle physics experiments used unpolarized beams and targets, thus 
completely overlooking the spins of the particles. This way only total cross-sections 
instead of spin- dependent cross-sections can be measured. Only a small portion of the 
reaction can thus be investigated. In the words of the originator of the spin filtering 
method of polarization buildup, which much of the investigations in this thesis are 



1 For further discussion on this see Figure 
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Figure 1-2: The upper beam is unpolarized, with equal number of particles in the 'spin 
up' and 'spin down' states, which in reality are randomly distributed. The lower beam is 
100% polarized, with all particles in the 'spin up' state. In practice the maximum beam 
polarizations achievable are about 90%. 

based on, P. L. Csonka [4] ; "One could, perhaps, say that the physicist who is able to 
measure only total cross-sections, is like the man in an art gallery who is only told the 
total weight of each statue, but is kept in ignorance of all other parameters specifying 
their shapes. Most of us would agree that he is missing something". Nowadays highly 
polarized beams of certain particles are possible. Polarized electrons and positrons 
have been used for many decades. Baryons have proved more difficult to polarize. 
The Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory, 
New York is the first accelerator to use a high energy polarized proton beam. 

1.2 A "spin crisis" in the parton model 

Prior to the European Muon Collaboration (EMC) [6] and the Spin Muon Col- 
laboration (SMC) [7] experiments at CERN it was assumed that all the spin of the 
nucleon was carried by its three constituent valence quarks. The startling results 
of EMC in 1988 [5] and 1989 \6\ showed that the three constituent valence quarks 
contribute very little to the spin of the nucleon. This caused, what was dubbed "The 
spin crisis in the parton model" [8j |9], prompting a new theoretical investigation into 
the spin structure of the nucleon, which continues to this day. The phrase "spin 
crisis" which endures to this day was coined in a beautifully titled paper "A crisis 
in the parton model: where, oh where is the proton's spin?" by Mauro Anselmino 
and Elliot Leader [8], and presented at the SPIN 1988 Symposium in Minneapolis, 
USA. The fact that the two original EMC papers were the most cited experimental 
papers in the field for three years and have a combined total of over 2500 citations 
shows the immense effort that has been afforded to solving the "spin crisis". It is 
now proposed that the spin of the nucleon is made up of the helicity of the con- 
stituent quarks Aq, the helicity of the gluons AG, the orbital angular momentum of 
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Figure 1-3: These diagrams describe the internal structure of protons and neutrons, accord- 
ing to the theory of Quantum Chromo- Dynamics (QCD). The diagram on the left shows 
the naive expectation that the spin of the nucleon is entirely constituted by the three valence 
quarks. The EMC and SMC results proved that this was not correct. The diagram on the 
right shows the current more complex view of nucleon spin structure, with contributions to 
the nucleon spin coming from the valence quarks, sea quarks, gluons and orbital angular 
momentum. 

the quarks L q , the orbital angular momentum of the gluons L g and the transversity 
of the quarks [referred to in the literature as either Sq, At q or hi q , we shall use the 
latter notation], as seen in Figure [U3l Gluons, being massless spin-1 bosons, cannot 
be transversely polarized, hence there is no gluon transversity. This gives us the 
longitudinal spin sum rule as follows: 

^Nucleon = \ = \ Aq + AG + + (1.2.1) 

where the superscript L refers to Longitudinal, and the transverse spin sum rule [TO] : 

^Nucleon = g = 2 ^ 9 + ^ + ^ ' (1-2.2) 

where the superscript T refers to Transverse. 

The current knowledge of these constituents is summarized in Figure 11-41 The 
contributions L q , L g and Aq are known from experiment [TTJ, [J_21 [J_3l E] and from 
Lattice QCD studies [15]. There are currently many theoretical models [16j and ex- 
perimental programs obtaining information on AG, these include HERMES, COM- 
PASS, JLAB and RHIC. But the last piece of the puzzle, the transversity distribution 
function hi q is to date almost completely unknown. In order to best measure hi q , a 
beam of polarized antiprotons would be required as we explain in the next section. 
Much of this thesis is devoted to a theoretical investigation of possible methods to 



Figure 1-4: The contributions to the spin of the nucleon: the helicity of the constituent 
quarks Aq, the orbital angular momentum of the quarks L q , the orbital angular momentum 
of the gluons L g all of which are known; the partly known helicity of the gluons AG and 
the unknown transversity of the quarks h\ q . 

polarize an antiproton beam in a storage ring. 

1.3 Antiprotons 

Antiprotons are the anti-particles of protons, which in turn are the core of the hy- 
drogen atom, the most abundant element in the Universe. The proton was shown 
to have an internal structure, i.e. not be a point particle, during seminal elastic 
electron-proton experiments in the Stanford Mark 3 accelerator, from 1954 to 1957 
[17] . Robert Hofstadter was awarded the 1961 Nobel prize for this ground-breaking 
discovery which ushered in a new era of investigation into the structure of the nu- 
cleon. A decade later, the much higher energy SLC accelerator was built at Stanford 
to investigate Deep Inelastic Scattering (DIS) experiments, showing the proton to be 
made up of point like quarks [T5] . Again this achievement warranted the Nobel prize, 
in 1990 to Jerome Friedman, Henry Kendall and Richard Taylor. The proton consists 
of uud valence quarks, hence the antiproton consists of uud valence quarks. Above, 
as in the rest of the thesis, we denote antiparticles by an over-bar. In shorthand 
notation protons are denoted p and hence antiprotons are denoted p. Antiparticles 
have the opposite electromagnetic charge of their corresponding particle, thus for 
chargeless particles (e.g. the photon) the antiparticle is the same as the particle. 
A proton has electromagnetic charge +1 in units where the electron charge is — 1, 
thus an antiproton has electromagnetic charge —1. While the concept of anti- matter 
often seems mysterious at first glance, it should be remembered that the positron 
(the anti-particle of the electron) was the third elementary particle discovered [T9"] . 
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after the electron and the photon. Hence anti-matter has been an integral part of 
physical theories since 1932, the same year the neutron was discovered. 

Antiprotons were discovered in 1955 by Owen Chamberlain and Emelio Segre [20J, 
who were awarded the 1959 Nobel prize for this ground-breaking discovery. Professor 
Chamberlain, who recently passed away, spent the rest of his life contributing great 
efforts to the investigation of polarization phenomena and spin physics in general 
[2T] . He was the first to investigate a possibility of polarizing antiprotons, and he 
co-organized the first workshop on polarizing antiprotons at Bodega Bay, California 
in 1985. 

The conclusions of this workshop |22j, were that a high intensity polarized an- 
tiproton beam was not achievable at that time. Another 22 years passed before the 
International community felt a sequel to this workshop was necessary, during which 
time interest in polarizing antiprotons grew steadily. There has been much recent 
interest in producing a high intensity beam of polarized antiprotons, starting in 2004 
with a proposal by the Polarized Antiproton experiments {VAX) Collaboration at 
GSI Darmstadt [23J. Since then many theories have been put forward on how to pro- 
duce such a beam. So in 2007 a sequel to the Bodega Bay workshop was organized in 
the newly founded Cockcroft Institute for accelerator research at the Daresbury Lab- 
oratory, UK [23] . A thorough investigation of the theoretical aspects of producing a 
polarized antiproton beam is presented in this thesis. 

Antiproton-proton colliders have played an important role in the advancement 
of High Energy Physics. In particular they led to the discovery of the W and Z 
bosons, and thus to the verification of the Weinberg, Glashow, Salam (1979 Nobel 
Prize) unified theory of electroweak interactions. This was done by the UAl and 
UA2 experiments at the SPS (Super Proton Synchrotron) collider in CERN in 1982 
and led to the 1984 Nobel Prize to be awarded to Carlo Rubbia and Simon van 
der Meer. One hopes that future polarized antiproton-proton colliders will lead to 
further epoch-making discoveries. 

1.4 Outline of the thesis 

The major theme of the thesis is a theoretical investigation of the spin filtering 
method of polarization buildup, and an application of this to producing a high in- 
tensity polarized antiproton beam. There is much debate in the International com- 
munity as to the correct theoretical description of spin filtering. We hope that the 
thorough analysis of spin filtering presented here will clarify some of this confusion. 
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No high intensity polarized antiproton beam has ever been achieved, and since a 
high intensity polarized antiproton beam could be used to measure many important 
quantities in particle physics, it is a main goal of the International community. We 
first calculate all polarization dependent cross-sections in QED for the processes of 
interest, then we develop a set of differential equations using these polarization de- 
pendent cross-section to describe spin filtering; finally numerical results are obtained 
from this formalism. 

In Chapter [2] the motivation for the thesis is outlined. The benefits of a polarized 
antiproton beam are described, as are all possible methods to produce such a beam. 
The methods to polarize bunches of other particles and atoms are also presented, 
such as electrons, positrons, protons, hydrogen and deuterium; and it is explained 
that none of these can be applied to the elusive case of polarizing antiprotons. It 
is concluded that spin filtering is the most promising method to produce a high 
intensity polarized antiproton beam and the chapter concludes with an overview of 
spin filtering and a description of how it has been verified experimentally. 

In Chapter [3] all electromagnetic helicity amplitudes and spin observables, ac- 
counting for polarization effects in spin 1/2 - spin 1/2 elastic scattering are cal- 
culated. Many of these results will be utilized in later chapters when providing a 
mathematical description of spin filtering, although their use is certainly not lim- 
ited to this. The spin 1/2 electromagnetic currents are introduced, both for point 
particles and particles with internal structure determined by electromagnetic form 
factors. A generic equation is derived that can be used to calculate all polarization 
phenomena in elastic spin 1/2 - spin 1/2 electromagnetic scattering to first order in 
QED. We then present results for all electromagnetic helicity amplitudes and spin 
observables for elastic spin 1/2 - spin 1/2 scattering. The spin- averaged differential 
cross-section for spin 1/2 - spin 1/2 scattering is also presented in a new compact in- 
variant form. These results are then presented in Chapter H] for the specific cases of: 
antiproton-proton, antiproton-electron and positron-electron scattering. Then the 
cross-sections and spin observables needed for spin filtering are explicitly presented, 
which will be utilized in the polarization evolution equations developed in Chapter 
The chapter concludes with a calculation of all spin - spin 1 helicity amplitudes, 
which describe the scattering of deuterons off a carbon nucleus for example. 

The theory of spin filtering is developed in Chapter A mathematical description 
of the related but simpler process of polarization buildup by the Sokolov-Ternov effect 
is first presented. The ideas presented are utilized in the mathematical descriptions 
of spin filtering which follow. The rates of change of the number of particles in each 
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spin state are combined into a set of polarization evolution equations which describe 
the process of polarization buildup by spin filtering. This set of polarization evolution 
equations is then analyzed and solved, emphasizing the physical implications of the 
dynamics. The chapter concludes with an investigation of spin filtering of a stored 
beam. 

Chapter O presents a thorough investigation of spin filtering under various alter- 
nate scenarios, which would be of interest to any practical project to produce a high 
intensity polarized antiproton beam. These scenarios are: 1) spin filtering while the 
beam is being accumulated, i.e. unpolarized particles are continuously being fed into 
the beam at a constant rate, 2) unpolarized particles are continuously being fed into 
the beam at a linearly increasing rate, i.e. the particle input rate is ramped up, 3) 
the particle input rate is equal to the rate at which particles are being lost due to 
scattering beyond the ring acceptance angle, the beam intensity remaining constant, 
4) increasing the initial polarization of a stored beam by spin filtering, and 5) the 
input of particles into the beam is stopped after a certain amount of time, but spin 
filtering continues. 

As an application of the theoretical work presented throughout the thesis a pos- 
sible method to produce a high intensity polarized antiproton beam by spin filtering 
off an opposing polarized electron beam is presented in Chapter [JJ It is also outlined 
how this work can be applied to polarizing antiprotons by spin filtering off a polarized 
hydrogen target. Firstly a description of the electron cooling technique to refocus 
the beam after scattering off the target each revolution in order to maintain high 
beam density is presented. Then the various experimental input parameters, such 
as revolution frequency, target areal density, target polarization and the effective 
acceptance angle; needed to obtain realistic numerical estimates from our mathe- 
matical formalism are each described. The benefits of using a lepton target are then 
described, before analyzing the case of spin filtering off an opposing polarized elec- 
tron beam. Finally spin filtering off a polarized hydrogen target is discussed, in the 
three cases of hydrogen with only electrons polarized, hydrogen with only protons 
polarized and finally hydrogen with both electrons and protons polarized. It is shown 
that electromagnetic effects dominate hadronic effects in pp scattering in the region 
of low momentum transfer of interest in spin filtering. 

In Chapter [S] some concluding remarks are presented. 
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1.5 Notation and conventions 



The conventions will mainly follow the book of Peskin and Schroeder [25]. Ratio- 
nalized units, where H — c — 1, will be used throughout the thesis unless otherwise 
stated. Units in this system are as follows: [length] = [time] = [energy] _ 1 = 
[mass] -1 . The usual convention of Greek characters representing four dimensional 
space-time indices {0,1,2,3} and Latin characters representing three dimensional 
space indices {1,2,3} is used throughout the thesis. The end of a proof is repre- 
sented by the symbol □, and = means that new objects are being defined. 

Lorentz 4-vectors are written 3 ), where x° is the time 

component and x 1 , x 2 , and x 3 are the x, y, and z space components respectively. 
Three dimensional vectors in Euclidean space are displayed in boldface, so that we 
can also write x = (x°, x). The Einstein summation convention, where repeated 
indices are summed over, is used throughout the thesis, unless otherwise specified. 
Our conventions for Dirac spinors are presented in Appendix [A] 

The Feynman slash notation = j^p^, and the Minkowski metric tensor 77 ^ = 
diag(+l, — 1, — 1, — 1) are used. The spin four vectors are normalized such that 
S^S^ = — 1. We use the shorthand notation A ■ B for the scalar product in 4- 
dimensional Minkowski space, where A-B = A^B^ = rj Ufl A u B fl = A B° — A-B. 
The totally antisymmetric tensor ^ v p° ' ( also known as the Levi-Civita symbol, is de- 
fined such that e 0123 = + 1 and e i23 = — 1, as seen in Appendix [Al 

Antiparticles are denoted by an over-bar. In shorthand notation protons are 
denoted p and hence antiprotons are denoted p. Electrons and positrons are denoted 
by e~ and e + respectively. Time increases from left to right in all Feynman diagrams 
throughout the thesis. Arrows on particle lines in Feynman diagrams denote the flow 
of particle number, which is forwards for particles and backwards for antiparticles. 

We denote the time variable in each of the dynamical systems by r to avoid confu- 
sion with the squared momentum transfer (Mandelstam t variable) used throughout 
the thesis. 

The scattering processes investigated in the thesis are always 2 — » 2 elastic pro- 
cesses, with the momentum and spin 4-vectors of each particle labeled as: 

A ( Pl , S,) + B (p 2 , S 2 ) — ► A (p 3 , S 3 ) + B (p 4 , S 4 ) , 

with the particles above being the beam (1), target (2), scattered (3) and recoil (4) 
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particles respectively. The four momentum transfer is defined as 

q = P3 ~ Pi = P2 - P4 ■ 

The helicity amplitudes are represented by 

M, (scattered, recoil ; beam, target ) = M. ( A 3 , A4 ; Ai, A2 ) . 

The arguments are to be read from right to left, as Ai and A 2 correspond to the 
incoming particles in the reaction and A3 and A4 correspond to the outgoing par- 
ticles in the reaction. For spin 1/2 - spin 1/2 scattering the helicities \ = ± for 
i E {1,2,3,4} are + if the particles spin vector points in the direction of its mo- 
mentum vector and — if the particles spin vector points in the opposite direction to 
its momentum vector. The ± in the helicity amplitudes are shorthand for ±1/2 , the 
helicity of a spin 1/2 particle. 

The spin observables are represented by K ab for the polarization transfer observ- 
able and ( 1 — D ab ) for the depolarization observable, where a, b G { X, Y, Z } for the 
direction of the particles spin vector where its momentum is along the Z direction. 
The subscripts are read from right to left, in e.g. K ab where b is the direction of the 
spin vector of the incoming particle and a is the direction of the spin vector of the 
outgoing particle. 
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Chapter 2 
Motivation 



"Polarization data has often been the grave- 
yard of fashionable theories. If theorists had 
their way, they might just ban such measure- 
ments altogether out of self-protection. " 

J. D. Bjorken 

In this chapter the motivation for the present work is discussed. The benefits to 
the high energy physics community of a high intensity polarized antiproton beam 
are first presented in section 12.11 by describing the important parameters in parti- 
cle physics that could be measured and investigated with such a beam. Section 12.21 
describes how bunches of other particles and atoms are polarized, such as electrons, 
positrons, protons, hydrogen and deuterium. Unfortunately none of these tried and 
tested techniques can be applied to the elusive case of polarizing antiprotons. Sec- 
tion 12.31 describes and compares some possible methods to polarize antiprotons. It is 
concluded that spin filtering is the most promising technique to produce a polarized 
antiproton beam as it is the only technique that has been experimentally verified. 
The chapter concludes with an overview of the theory of spin filtering, and a section 
showing how spin filtering was verified for polarizing a proton beam by repeated scat- 
tering off a polarized hydrogen target in a storage ring by the FILTEX experiment 
in 1993. 
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2.1 Motivation for a polarized antiproton beam 



A high intensity polarized antiproton beam would provide the unique possibility to 
measure many very important quantities in particle physics. The most important 
quantity that could be measured is the transversity distribution of quarks inside 
protons, which has eluded direct measurement thus far. Two other important in- 
vestigations, into single spin asymmetries and nucleon electromagnetic form factors, 
can also be greatly advanced if a high intensity polarized antiproton beam was avail- 
able. These motivations for producing a high intensity polarized antiproton beam 
are described in detail below. 



2.1.1 The transversity distribution function 

The transversity distribution function is the last leading twistQ piece of the QCD 
description of the partonic structure of the nucleon, in the collinear limit f|, that has 
not been directly measured. It describes the quark transverse polarization inside a 
transversely polarized nucleon. In fact, to date, almost nothing is known about the 
transversity distribution, except for the recent work of Anselmino et al. [261 EZ] • Un- 
like the other leading twist distributions [the unpolarized quark distribution q (x, Q 2 ) 
and the helicity distribution Aq (x, Q 2 )} which have been measured, the transversity 
hi q (x, Q 2 ) [sometimes referred to in the literature as A T q (x, Q 2 ) or 5q (x, Q 2 )] can 
neither be accessed in deep inelastic scattering of leptons off nucleons, nor can it 
be reconstructed from the knowledge of q(x,Q 2 ) and Aq (x, Q 2 ) [23]. In a trans- 
versely polarized hadron, h\ q (x,Q 2 ) is the number density of quarks with momen- 
tum fraction x and polarization parallel to that of the hadron, minus the number 
density of quarks with the same momentum fraction and antiparallel polarization, 
i.e. h lg (x,Q 2 ) = q^ (x,Q 2 ) — q^(x,Q 2 ) [25J. One cannot claim to understand the 
spin structure of the nucleon until all three leading twist structure functions have 
been measured. 

In order to best access the transversity distribution function, the double spin 
asymmetry Att in the Drell-Yan production of lepton pairs must be measured; thus 
both initial particles in a reaction must be transversely polarized. It could in future 
be done for p^ p^ scattering at RHIC, but this asymmetry is expected to be small 
from theory [29J , as explained below. Also the cross-section for Drell -Yan lepton pair 



1 Leading twist means that in the factorization of a physical process the parton distribution 
function appears in the leading order of 1/Q 2 . 

2 The collinear limit is where the intrinsic transverse motion of the quarks is averaged over. 
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Figure 2-1: The Drell-Yan lepton pair production process pp —tllX, via a virtual photon 
7*- 

production is much higher for pp scattering than for pp scattering (a D ^f 3> c^y); 
because in the former case valance quarks in the proton annihilate with valance 
antiquarks in the antiproton, as opposed to with sea antiquarks in the second proton 
in the latter case. The Drell-Yan lepton pair production process [31] is shown in 



Figure 12-11 The double spin asymmetry, an experimentally measurable quantity, is 
defined as 



Aj = j_-m- , T~+\ ' (2.1.1) 



da n - d<7^ 
dcxtt + dcxti 

where % and j can be either L for longitudinal, or T for transverse. For p^ p^ Drell-Yan 
processes the transverse double spin asymmetry is [29] 

pp = dAa E q e 2 q [hf q (x u M 2 ) hf q (x 2 ,M 2 ) + hf q (x u M 2 ) K q {x 2 ,M 2 )] 
TT da J] ? e 2[gP(xi,M 2 ) gf(x 2 ,M 2 ) + gf(xi,M 2 ) gP(x 2 ,M 2 )] ' 

w dAa hf u ( Xl ,M 2 ) K u (x 2 ,M 2 ) + hf a ( Xl ,M 2 ) hf u (x 2 ,M 2 ) 

~ da u*(x u M 2 ) u p {x 2 ,M 2 ) + u*( Xl ,M 2 ) up{x 2 ,M 2 ) ' 1 ' ' ' 

where dAa and da are the polarized and unpolarized cross- sections of the elementary 
QED process qq — > l~ l + respectively, M is the invariant mass of the lepton pair, 
e q is the electromagnetic charge of the quarks and x\ and x 2 are the fraction of 
their respective nucleon momentum carried by each of the interacting partons. The 
leading term in the approximation comes from the fact that the u quark dominates 
at large x [22] • Whereas for p* p^ Drell-Yan processes 
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dAa T. q e 2 q [K q {x u M 2 ) hf q (x 2 ,M 2 ) + hf q (x u M 2 ) hf q (x 2 ,M 2 )} 
d ^ E g e 2 [qP(x u M 2 ) qP(x 2 ,M 2 ) + qP(x u M 2 ) qP(x 2 ,M 2 )] 

dAa h p lu (x u M 2 ) h? u (x2,M 2 ) 



da up(x 1 ,M 2 ) up(x 2 ,M 2 ) 



(2.1.3) 
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the latter of which is much larger since there are more antiquarks in antiprotons mak- 
ing hf u (x,Q 2 ) = hf u (x,Q 2 ) » h? u (x,Q 2 ) = h( u {x,Q 2 ) [29]. Thus A!&, which 
can only be measured using a polarized antiproton beam, is expected to be much 
bigger than A^. Note in the literature all quantities in these equations are often 
written with respect to the proton using the fact that the distribution of antiquarks 
in a proton is equal to the distribution of quarks in an antiproton etc. but here we 
want to keep the antiproton distribution functions explicit. Using the above fact, 
and at x\ = X2, eq. (12.1. 3p reduces to 



APP 

yirpr-j-i 



dAa 



da u(x u M 2 ) 



(2.1.4) 



providing a unique direct way to measure a single transversity distribution func- 
tion [22]. 

Also q (x, Q 2 ) and q (x, Q 2 ) decrease with increasing x, so to measure A TT large 
x\ and X2 is favoured [291 ED] • Interestingly this happens for lower energy scattering 
again making a low/medium energy facility, such as that proposed by the VAX 
Collaboration, more suited than RHIC^l At RHIC energies, even though A TT could 
be detected it only measures the transversity of the sea quarks, at the lower VAX 
energies we could investigate the transversity of the valence quarks 



2.1.2 Single spin asymmetries 

Single Spin Asymmetries (SSA), where one of the initial particles in the reaction is 
polarized in the direction of the arrows below, are defined as 

dtr 1 " - da± . 

A = — — , 2.1.5 

J do-t + da^ ' y J 

where j can be either L for longitudinal, or T for transverse. 

Some data on SSA in p Drell-Yan lepton pair production was obtained by the 
E704 experiment at Fermilab [321 [37], but because the collisions were in the energy 
region of J/ip production^ it was difficult to distinguish the Drell-Yan signal from 



3 Note there is a proposal to run RHIC at y/s — 50 GeV instead of their usual y/s = 200 GeV 
which would make it suitable in this regard, but the problem of no antiproton beam would still 
remain. In Run 6 (2006) RHIC used </s = 62.4 GeV. 

4 The J/tp particle is the first excited state of charmonium, a meson consisting of one charm 
quark and one charm antiquark. Two papers by separate experiments announcing its discovery 
were published on the same day, one group naming it the J particle and the other group naming it 
the ip particle. It has come to be known as the J/tp particle. 
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the large J /if) production background. The low intensity polarized antiproton beam 
used in E704 is described in section 12.3.11 

Importantly analyzing charm production in p^ p scattering will make it possible 
to disentangle the Sivers [331 EI] and the Collins mechanisms [35] , of which there is 
great theoretical interest. In general, both effects contribute to the measured SSA, 
but in the case of charm production the Collins mechanism drops out. 

A polarized antiproton beam would allow further analysis of single spin asymme- 
tries in p^ p scattering, augmenting the brief Fermilab data on this [32], and adding 
to the current data on single spin asymmetries which have been observed in p p^ 
and p^ p reactions [36, E7J EFJ and the double spin asymmetries observed in p^ p^ 
reactions at RHIC [39|. These observed asymmetries are very large, up to 40% [36J, 
prompting Stan Brodsky to call them "the greatest asymmetries ever seen by a hu- 
man being" constituting "one of the unsolved mysteries of hadron physics". There is 
much current interest in the theoretical community to try to achieve a satisfactory 
understanding of these large single spin asymmetries [381 STil El H21 H31 144"] . 

2.1.3 Electromagnetic form factors of the proton 

The fact that nucleons (protons and neutrons) are not point particles and have an 
internal structure, is parameterized into electromagnetic form factors, as treated 
later in section I3TT1 The Sachs electric and magnetic form factors Ge and Gm con- 
tain information on the finite charge radius of the proton, thus are very important 
components of a complete understanding of particle physics. They can be measured 
experimentally but there is not, to date, complete agreement between the experi- 
mental results and theoretical models of the form factors [HI HBJ . 

There is much current theoretical interest in nucleon time-like form factors [47] . A 
polarized antiproton beam would enable the first measurement of the moduli and the 
relative phase of the time-like electric and magnetic form factors Ge and Gm of the 
proton. An unexpected Q 2 = — q 2 dependence of the Ge(o 2 ) / Gm(q 2 ) ratio of the 
electric and magnetic form factors of the proton, has been observed at the Jefferson 
laboratory (JLAB), the ratio decreasing monotonically with increasing Q 2 [451145] . It 
would be possible to clarify this unexpected Q 2 dependence by a measurement of the 
relative phases of Ge{q 2 ) and Gm{q 2 ) in the time-like region, which would constrain 
and discriminate strongly between the models for the form factors. This phase can 
be measured for the first time in the reactions p^ p — > e + e~ and pp^ — )■ e + e~ [47], the 
former of which is uniquely possible with a polarized antiproton beam. The JLAB 
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data was obtained by analyzing the polarization transfer reaction pe~' — > p'e~ } this 
data could be augmented and checked by analyzing the polarization transfer reaction 
pe~^ — > e~ , which is at the heart of the spin filtering technique discussed in detail 
throughout this thesis. 

The relative phase ambiguity can also be addressed by measuring the double spin 
asymmetry in the reaction p^ p^ — > where I is any lepton. This reaction can 

also be used to analyze the Ge — Gm separation, thus serving as a check of the 
Rosenbluth separation in the time-like region [50J. 

2.2 Polarizing bunches of particles or atoms 

High intensity beams of polarized electrons, positrons and protons, as well as polar- 
ized atomic gas targets have been used in high energy physics laboratories throughout 
the world. We now briefly describe how they are polarized. 

• Sokolov-Ternov effect ('radiative' or 'self '-polarization): A beam of charged 
particles circulating in a storage ring is automatically polarized because of 
a difference in the spin-flip transition rates due to emission of photons by 
synchrotron radiation induced by bending in the magnetic field of the ring. 
This method works well for polarizing electrons and positrons, but not for 
heavier particles such as protons or antiprotons, as explained in section f)2.3.3p . 

• Atomic hydrogen and deuterium are polarized by removing atoms in certain 
hyperfine states, and inducing angular momentum conserving transitions be- 
tween hyperfine states. 

• Once hydrogen is polarized the electrons can be stripped off in a magnetic field 
leaving polarized protons. 

Unfortunately it is not possible to produce a high intensity polarized antiproton 
beam using any of these tried and tested methods, as to do so one would need a 
large supply of antihydrogen atoms. 

The Sokolov-Ternov effect will be described in section (I2.3.3|) and we describe 
how polarized hydrogen is obtained in section (12.2.ip . 
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2.2.1 Polarizing hydrogen gas 

Unpolarized hydrogen atoms in a strong magnetic field equally populate each of four 
hyperfine states: 

I tp te ) | tp le ) | Ip te ) | tp te ) 

An inhomogeneous magnetic field acts as a Stern-Gerlach apparatus separating the 
atoms in the states | tp te ) an d | | p te ) from those in the states | | p | e ) and | tp le ) • 
A sextupole magnet focuses the atoms in one pair of states while defocusing the 
others. Thus one can extract atoms in the states | tpte) an d | Ipte), ^.e. hydrogen 
atoms in which the electrons are totally polarized, but protons unpolarized. If one 
then requires hydrogen atoms in which the protons are totally polarized but the 
electrons unpolarized, angular momentum conserving transitions from | tp te ) to 
| tp le ) can be induced by a radio frequency field. 

Hydrogen with both the electrons and protons polarized can be obtained by 
isolating the | tpte) state, but with only half the intensity of hydrogen with either 
electrons or protons polarized. 

In summary there are three types of polarized hydrogen, with all atoms in the 
hyperfine states as follows 

I tpte) + Itple) V p = 1 and V e = (2.2.1) 

I tpte) + I tpte) V p = and V e = 1 (2.2.2) 

I tpte) V p = 1 and V e = 1 (2.2.3) 

where we denote the polarization of the protons in the hydrogen by V p and the polar- 
ization of the electrons in the hydrogen by V e ■ In practice the atoms are not perfectly 
isolated in certain hyperfine states, thus the electron and proton polarizations in po- 
larized hydrogen are less than one. The HERMES Collaboration at DESY have 
utilized polarized hydrogen and polarized deuterium targets with V v = 0.9 and/or 
V e = 0.9 [51], and these targets are now being used by the VAX Collaboration in 
COSY Jiilich for preliminary tests on spin filtering [521 03] • 

Other atoms, such as deuterium, can be polarized analogously. Stripping these 
atoms of electrons in a magnetic field leaves a polarized ion beam. 

It is not possible to generate a beam of polarized antiprotons by this means as it 
is, thus far, not possible to accumulate large numbers of antihydrogen atoms. 
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The RHIC collider in Brookhaven National Laboratory, New York is the world's 
first high intensity polarized proton accelerator. A high intensity polarized antiproton 
collider would greatly supplement and add to results obtained at RHIC. 

2.3 Methods to polarize an antiproton beam 

Now that we have demonstrated the incredible potential of a high intensity polarized 
antiproton beam, let us investigate the various methods of generating such a beam. 
Physicists have been trying to produce beams of polarized antiprotons for over 25 
years, a great summary of proposed methods is given in Ref. [22J. Atomic beam 
sources, used in the production of polarized protons and heavy ions will not work 
because of the annihilation of antiprotons with matter. The E704 experiment at 
Fermilab has produced polarized antiprotons from the decay of polarized A hyperons, 
but the intensities achieved were too low for current needs. Storing antiprotons in a 
storage ring would help build up to a high luminosity beam. 

Spin filtering has been proven to work for protons scattering off a polarized in- 
ternal hydrogen target in the FILTEX experiment at the TSR ring in Heidelberg 
in 1992-1993 [53]. Thus spin filtering is the only plausible experimentally tested 
technique for generating a high intensity polarized antiproton beam. In light of this 
we devote much of this thesis to the theoretical understanding of spin filtering in 
general. As an application of our theoretical work we propose a method to polarize 
antiprotons by spin filtering off an opposing polarized electron beam, and calculate 
the polarization buildup time and maximum polarization possible in this case. 

2.3.1 Antihyperon decay 

Antihyperons are produced when a Multihundred-GeV proton beam strikes a target. 
The antihyperons decay into antiprotons, which should have the same polarization 
as the protons from hyperon decay. A polarized antiproton beam of this type was 
produced at Fermilab's E704 experiment [55]. The low intensity (because it is a 
tertiary beam) and large phase space made it difficult to store and accelerate these 
polarized antiprotons; however it was possible to scatter them off a polarized or 
unpolarized proton target. The polarization of the antiprotons comes from parity 
violating decays of antilambdas, and the measured polarization was as high as 64%. 
The target they used to produce the antilambdas was Beryllium, and their polarized 
antiproton beam intensities were up to 1.5 x 10 5 s _1 [55]. This method of producing 
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polarized antiprotons is not suitable for high luminosity experiments, such as the high 
intensity beam a storage ring could provide, which are needed to access transversity 
and other measurements as outlined in section 12.11 

2.3.2 Stern-Gerlach separation 

A possible method to produce a polarized antiproton beam from an unpolarized 
antiproton beam is based on the Stern-Gerlach effect. In an inhomogeneous magnetic 
field the spins of particles, aligned parallel or antiparallel to the field, are deflected 
in opposite directions and become spatially separated. For this reason this method, 
proposed by Niinikoski and Rossmanith in 1985 [56J, is also called the spin-splitter 
technique. A major advantage of this method is that the beam can first be accelerated 
to any desired energy and then polarized, thus avoiding the loss of polarization 
associated with accelerating polarized beams f|. 

In a typical storage ring inhomogeneous magnetic fields are provided by the 
quadrupoles. It was hoped that the spatial separation of the particles in the two 
spin states would add up on passing through many quadrupoles, and further over 
many revolutions in the storage ring; eventually leading to a macroscopic separation 
of the particles in opposite spin states [56]. One spin state can then be dumped, or 
flipped, and one is left with a polarized antiproton beam. 

Unfortunately, after much interest in this technique [5TI 158] , the International 
Community has doubts as to whether effects in successive quadrupoles will add up 
coherently [59]. The effects may continuously cancel each other out and one will be 
left with no net separation of particles in the two spin states. At the very least this 
method would have to be experimentally verified before being considered a practical 
method of producing a polarized antiproton beam. 

2.3.3 Spontaneous synchrotron radiation emission 

Charged particles emit synchrotron radiation, in the form of photons, when bent in a 
magnetic field. There is a slight difference in the spin-flip transition cross-sections due 
to this photon emission: a ( — ?• 7 ) ^ a ( — > 7 ) , thus over time the beam 
of charged particles acquires some polarization. The cross-section for a particle in the 
'spin up' state to flip to the 'spin down' state on emitting a photon by synchrotron 

5 Beams tend to lose some of their polarization at certain depolarization resonance energies during 
acceleration. This problem can be circumvented by utilizing Siberian Snakes [BD], devices which flip 
the polarization vector of each beam particle by 180 degrees each revolution. Thus any deflections 
from the polarization axis are canceled out every two revolutions. 



19 



+ 




7 




7 
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Figure 2-2: The two spin-flip Feynman diagrams that contribute to the Sokolov-Ternov 
effect. In (a) an electron in the 'spin up' state gets flipped to the 'spin down' state, while 
in (b) an electron in the 'spin down' state gets flipped to the 'spin up' state; after emitting 
a photon due to synchrotron radiation induced by bending in a magnetic field. The cross- 
sections for these two processes are not equal and as such there will be a gradual buildup of 
polarization in the beam, known as the Sokolov-Ternov effect. 

radiation is different to the cross-section for a particle in the 'spin down' state to 
flip to the 'spin up' state. The 'spin up' and 'spin down' states are defined as the 
particle's spin being aligned parallel and antiparallel to the magnetic field of the 
storage ring respectively. This 'self polarization' is called the Sokolov-Ternov effect 
after the Russian theorists who discovered it around 1963 |61j . 

In a perfect storage ring an equilibrium polarization of Vst = 8/(5 \/S ) ~ 0.924 
is reached [H2]- In practice the maximum polarization achievable is slightly less 
than this ideal value due to imperfections of the magnetic fields in the storage ring. 
However in less than one hour electron beams at TRISTAN in Japan and HERA 
in Germany acquired polarizations of about 80% or more [62]. The Sokolov-Ternov 
radiative polarization is along the vertical direction perpendicular to the storage ring 
plane. 

The effect is much stronger for electrons than for protons as the rate of syn- 
chrotron radiation (number of photons emitted per second) is related to the velocity 
of the particle not its energy. Because (anti)protons are approximately 1800 times 
more massive than electrons, at a given energy electrons are traveling at a much 
higher velocity, i. e. much closer to the speed of light (j e = ( m p / m e ) 7 p ^ 1800 7 P ). 
If (anti)protons were moving this close to the speed of light they too would become 
self polarized by the Sokolov-Ternov effect. Thus at a given energy the time taken to 
polarize electrons by the Sokolov-Ternov effects is much shorter than the time taken 
to polarize (anti)protons. Even at the 20 TeV of the proposed Superconducting Su- 
per Collider (SSC) it would take antiprotons or protons about 10 7 years to acquire a 
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Electron storage 


rings 


Proton storage rings 




LEP 


TRISTAN 


HERA 


LHC 


SSC 


E [GeV] 


46.5 


30 


27.521 


7000 


20000 


7 = V l ~P 2 


90999 


58710 


53858 


7462 


21317 


R [m] 


4243 


480 


1008 


4243 


12096 


p H 


3096.2 


246.5 


575 


3096.2 


10108 




6039 


3896 


3574 


567 


1429 


r ST [min] 


308 


2 


35 


2.8 x 10 14 


4.5 x 10 13 



Table 2.1: Properties of some high energy electron storage rings, and proposed proton 
storage rings. E is the kinetic energy of the beam, 7 = yT — /3 2 is the relativistic Lorentz 
factor, R the mean radius of the storage ring consisting of identical bending magnets of 
bending radius p separated by straight sections combining to give a total circumference 
2 ttR, iV 7 is the average number of photons emitted per particle per revolution and tst is 
the time taken to reach the equilibrium Sokolov-Ternov polarization. Parts of this table are 
reproduced from Refs. and fU3jj . As one can see this method of polarization buildup 
takes too long for (anti)proton rings. 



useful polarizatioro, and much higher kinetic energies would be required to provide 
a practical method of polarizing antiprotons by the Sokolov-Ternov effect. 

The time taken to reach the equilibrium polarization is given by [6T| l62| 164"] : 

8 mp 2 R 

^st = — -7= z 5 , (2.3.1) 

5V3 rh>y 5 

where m is the particle's mass, r the classical radius of the particle (electron or 
proton), h = h / 2 7r the reduced Planck's constant, R the mean radius of the storage 
ring consisting of identical bending magnets of bending radius p separated by straight 
sections combining to give a total circumference 2 ttR and 7 = 1/ yl — (3 2 is the 
relativistic Lorentz factor, where /3 = v / c is the ratio of the particles velocity to that 
of light. Some properties of current and proposed future synchrotrons are presented 
in Table EU 

The Sokolov-Ternov effect is similar to systems investigated later in this thesis, 
and it can be described by systems of differential equations similar to ones we develop 
to describe spin filtering. Hence to provide a comparison we present and solve a set of 
differential equations describing the Sokolov-Ternov effect in section I5TT1 Synchrotron 
radiation is the physical principle behind the antenna, emitting photons in the form 



6 See Table O 
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of radio waves, and some lasers generated by wiggler magnets. The light produced by 
synchrotron radiation is used by biologists and many storage rings have a second life 
as intense light sources after high energy physics experiments have ceased. Electrons 
lose much of their energy in a storage ring due to synchrotron radiation, typically 
emitting hundreds to thousands of photons per revolution 5 , and as a consequence 
very high energy electron/positron accelerators, such as the International Linear 
Collider (ILC), must be linear to avoid this problem. 

We conclude that this method of polarization buildup would take too long for an 
antiproton beam to be considered practical at present energies. 

2.3.4 Polarization of directly produced antiprotons 

It is well known that the particles produced when a high energy proton beam strikes 
a target have some polarization at certain production angles. Some of the parti- 
cles produced will be antiprotons; in fact this is how antiproton beams are obtained 
[65] . Unfortunately the polarization generally seems to be larger at larger production 
angles where the cross-sections are smaller [22]. Thus it appears difficult to simulta- 
neously obtain antiprotons with a high polarization and a high beam intensity using 
this method. 

2.3.5 The theory of spin filtering 

The spin filtering method of polarization buildup [UEUEE], described schematically 
in Figure !2-3[ consists of a circulating beam repeatedly interacting with a polarized 
internal target in a storage ring. Originally proposed by P. L. Csonka in 1968 [I], it 
is based on the selective removal of particles from the beam, and selective spin-flip 
while remaining in the beam, due to spin-dependent scattering off a polarized target. 
Many particles are scattered at small angles but remain in the beam after refocusing 
each revolution. This introduces a characteristic laboratory frame acceptance angle 
# acc , scattering above which causes particles to be lost from the beam. There is also 
a minimum laboratory frame scattering angle 9 m \ ni corresponding to the Bohr radius 
of the atoms in the target, below which scattering is prevented by Coulomb screen- 
ing |J7] , as described in section 17.2.51 The two physical processes that contribute to 
polarization buildup by spin filtering, as described in Figure !2-4[ are: 

(a) spin selective scattering out of the beam, and 

(b) selective spin-flip while remaining in the beam 
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Figure 2-3: This diagram describes the spin filtering technique. Beam particles travel 
along the beam axis and scatter off the target. Particles scattered at angles greater than the 
acceptance angle 9 acc are lost from the beam, while particles scattered at angles less than 
$acc P a ss through a beam focuser and remain in the beam. In this simplistic diagram the 
beam focuser is represented by a lens, but in reality the beam is focused by electron cooling 
as explained in section \n\ 

Thus particles in one spin state may be scattered out of the beam, or have their 
spin-flipped while remaining in the beam, at a higher rate than particles in the other 
spin state. Hence over time one spin state is depleted more than the other leading 
to a beam polarization. The beam will diverge slightly after many interactions with 
the target, but can be refocused by beam cooling, as explained in section 17.11 We 
prove later in the thesis that beam cooling does not depolarize a stored antiproton 
beam. 

As the beam polarization increases the beam intensity decreases, when there is 
scattering out of the beam. So one can obtain beam polarization at the expense of 
losing beam intensity. Low beam intensity means low event rate, hence low statistics 
in a measurement, which is never desired. This trade-off between beam polarization 
and beam intensity is characteristic of spin filtering and must be optimized to produce 
a sufficient beam polarization while maintaining reasonable beam intensity. 

An advantage for the spin filtering method is that polarized hydrogen and deu- 
terium jet targets have already been developed for other projects. Highly polarized 
high density gas jet targets have been used in the HERMES and COMPASS ex- 
periments. The HERMES experiment has been decommissioned since the shutdown 
of the HERA accelerator complex in DESY, and the polarized gas target has been 
transfered to COSY in Julich, Germany to be used in spin filtering studies. It is likely 
that the HERMES polarized gas target will be used by the VAX Collaboration in a 
future spin filtering Antiproton Polarizer Ring at FAIR, GSI Darmstadt. 
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Figure 2-4: The following two diagrams provide a schematic representation of the 
two physical processes, selective scattering out of the ring (left) and selective spin- 
flip (right), that contribute to polarization buildup by spin filtering in a storage ring. 
Particles in the 'spin up' state are represented by blue squares and particles in the 
'spin down' state are represented by yellow squares, while the grey box represents a 
polarized target. In both cases the beam is initially unpolarized with equal numbers 
of particles in the 'spin up ' and 'spin down' states. 
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Selective scattering out of the ring (left): When interacting with the polarized 
target at certain energies particles in the 'spin up ' state are scattered out of the beam 
at a higher rate than particles in the 'spin down' state, hence the larger blue arrow 
than yellow arrow. Thus one is left with a beam that has more particles in the 'spin 
down' state, i.e. the beam is now polarized, represented by the excess of yellow squares 
in the final beam. Note that since particles have been scattered out of the ring there 
are less particles in the beam after interaction than were in the beam initially, this is 
represented by the smaller final beam. If the target was unpolarized particles in both 
spin states would be scattered out of the beam at equal rates, thus no polarization 
buildup would occur via this process. 

Selective spin-flip (right): On interaction with the polarized target at certain 
energies, the 'spin up' to 'spin down' spin-flip cross-section is larger than the 'spin 
down' to 'spin up' spin-flip cross-section. We represent this by different size arrows 
with colours fading from blue to yellow and from yellow to blue respectively. Thus 
after interaction with the target the beam will have more particles in the 'spin down' 
state than in the 'spin up' state, i.e. the beam is now polarized, represented by the 
excess of yellow squares in the final beam. Note that the beam intensity is the same 
after interaction with the polarized target in this process since particles are not lost 
from the beam, they are just flipped from one spin state to the other. If the target 
was unpolarized particles in both spin states would have their spins flipped at equal 
rates, thus no polarization buildup would occur via this process. 
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There has been much debate amongst theorists as to what mechanisms are re- 
sponsible for the polarization buildup in spin filtering. Contributions come from 
the electromagnetic scattering of beam antiprotons off the electrons in the hydrogen 
target and from the electromagnetic and hadronic scattering of beam antiprotons 
off the protons in the hydrogen target. Horowitz and Meyer, in 1994, were the first 
to highlight the importance to spin filtering of the electrons in the hydrogen target 
[BSl EH] • They claim that electrons in a hydrogen target are not massive enough to 
deflect antiprotons beyond the acceptance angle of any storage ring, a fact which we 
demonstrate later in the thesis. Thus scattering of the antiprotons off the electrons 
in a hydrogen target causes no beam losses and any polarization buildup must be due 
to spin-flip transitions [70]. In 2005 two groups from the Budker Institute for Nuclear 
Physics, Russia and the Institute for Nuclear Physics, Julich, Germany claimed that 
such spin-flip effects are small thus spin filtering off polarized electrons in a hydrogen 
target will lead to a negligible rate of polarization buildup [TTJ (72J E3]- An exper- 
iment has been proposed to test this claim [71], by investigating the converse case 
of whether unpolarized electrons in a helium-4 target depolarize a stored polarized 
proton beam. The helium-4 target is chosen because its nuclei is spin - 0, hence any 
polarization transfer must come from scattering off its electrons. There are currently 
two schools of thought regarding spin filtering of antiprotons off a polarized hydro- 
gen target, (1) proposal [SS] building on the work of Horowitz and Meyer, which 
advocates using a hydrogen target with high electron polarization and low proton 
polarization; and (2) the Budker/ Julich proposal [TTJ [72] to use a hydrogen target 
with low electron polarization and high proton polarization. As is often the case the 
matter must be resolved by an experiment [74] to see which method is preferable. 

There are many advantages of using a lepton target instead of an atomic target, 
the foremost of which is that antiprotons cannot be absorbed by a lepton target 
as they are in a atomic target due to annihilation with the protons in the atomic 
target. This fact has led to two proposals for spin filtering off polarized lepton 
beams: one off a co-moving polarized positron beam [75] by a group in Mainz, 
Germany and the other, presented in this thesis, off an opposing polarized electron 
beam [76]. The momentum of an opposing electron beam causes antiprotons to be 
scattered beyond acceptance, hence allowing contributions from both of the physical 
process, selective scattering out of the beam and selective spin-flip, of spin filtering. 
A thorough treatment of the dynamics of spin filtering has been presented recently 
by the present author [TTJ [78] and forms much of the later chapters of this thesis. 



25 



Spin filtering is the only method to produce a polarized antiproton beam in a 
storage ring that has been successfully tested, by the FILTEX experiment in 1993 
[53] , as described in section 12.41 As a result, much of this thesis is devoted to a 
theoretical understanding of the spin filtering process, under various scenarios. 



2.4 Verification of spin filtering 

Polarization buildup by spin filtering has been proven to work in the FILTEX ex- 
periment at the Test Storage Ring (TSR) at the Max Planck Institute for Nuclear 
Physics in Heidelberg, Germany [53] . We summarize their results below. 
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Figure 2-5: The results of the FILTEX experiment, showing polarization buildup over 
time. The solid lines show the best fit to the data with a rate of polarization buildup of 
1.24 x 10 ~ 2 h' 1 . The dashed lines are based on the expected buildup rate from the model 
presented in Ref. \5$ , from where this plot has been reproduced with permission from the 
authors. 

In the TSR a 23 MeV proton beam was stored and repeatedly made to interact 
with a polarized internal hydrogen gas target. A polarized hydrogen target with 
atoms in the hyperfine state | tpte) wa s used, i.e. where both the protons and 
electrons are polarized. The target density was 6 x 10 13 polarized hydrogen atoms 
per cm 2 , and the frequency of revolution was 1.177 MHz [M]. The beam was left 
to orbit in the ring passing through the target each revolution for times between 30 
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and 90 minutes; then the polarization was measured. The proton beam was initially 
unpolarized and over time it gained a small amount of polarization as shown in 
Figure [223 

The polarization buildup rate of the proton beam at FILTEX was [51] : 

— = 0.0124 ± 0.0006 per hour. (2.4.1) 
d r 

After 90 minutes the polarization had increased to 1.86% and the beam intensity 
had decreased to 5 % of its original value [51] . But with a better configuration of 
the experiment, and a dedicated spin filtering polarizer ring, the rate of polarization 
buildup could be greatly increased. The TSR ring had an acceptance angle measured 
to be # acc = 4.4 ± 0.5 mrad, which we could optimize for our needs. 

The beam lifetime r*, which we discuss later in the thesis, is the time taken 
for the number of particles in the beam to decrease by a factor of e = 2.78. The 
beam lifetime in the TSR during the FILTEX experiment, with the polarized internal 
target in the ring, was 30 minutes. We show in section 15.3.11 that the polarization 
achieved after two beam lifetimes is an important measure of a spin filtering scheme. 
At FILTEX this value was measured to be Pbeam ( 2 ) = 0.0124. 

This was just a feasibility test for spin filtering, and while it verifies that the 
method works, the polarization buildup rate was small. In order to maximize the 
effect of spin filtering a dedicated spin filtering polarizing ring would need to be built. 
The VAX Collaboration has recently proposed the construction of such a ring called 
the Antiproton Polarizer Ring (APR) inside the HESR at FAIR in GSI Darmstadt, 
Germany. 
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Chapter 3 

Generic helicity amplitudes and 
spin observables 

"The most incomprehensible thing about the 
world is that it is comprehensible. " 

Albert Einstein 

In this chapter all electromagnetic helicity amplitudes and spin observables, ac- 
counting for polarization effects in spin 1/2 - spin 1/2 elastic scattering are calculated. 
Many of these results will be utilized in later chapters when providing a mathemat- 
ical description of spin filtering, although their use is certainly not limited to this. 
We begin in section [37X1 by introducing the spin 1/2 electromagnetic currents, both 
for point particles and particles with internal structure determined by electromag- 
netic form factors. A generic equation is derived in section 13.21 that can be used to 
calculate all polarization phenomena in elastic spin 1/2 - spin 1/2 electromagnetic 
scattering to first order in QED. The spin- averaged differential cross-section for spin 
1/2 - spin 1/2 scattering is presented in a new compact invariant form in section 1331 
We then present results for all electromagnetic helicity amplitudes for elastic spin 
1/2 - spin 1/2 scattering in section [374]. and for electromagnetic spin observables for 
elastic spin 1/2 - spin 1/2 scattering in section l3~5l 
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3.1 Spin 1/2 electromagnetic currents 



In this section the spin 1/2 electromagnetic currents of both point-like particles 
and non-point- like particles, with internal structure defined by electromagnetic form 
factors, are introduced. We investigate 2 particle — > 2 particle elastic scattering 
processes in the space-like region, with the mass and the momentum and spin 4- 
vectors of each particle labeled as: 

A ( M, pi, ft) + B (m,p 2 , S 2 ) — ► A (M,p 3 , S 3 ) + B (m, p 4 , S 4 ) , 

with the particles above being the beam (1), target (2), scattered (3) and recoil (4) 
particles respectively. The four momentum transfer is defined as 

Q = P3 ~ Pi = P2 - Pi- (3.1.1) 

For 2 particle — > 2 particle elastic scattering the spin-averaged differential cross 
section is related to the helicity amplitudes -M(A 3 A 4 ; Ai A 2 ) by 

v ' Ai A2 A3 A4 

where Ai, A2 and A3, A4 are the helicities of the initial and final particles respectively, 
sa and sb are the spins of the two particles in the elastic process, and the s and t are 
Mandelstam variables [75] defined in Appendix [B] We label the mass of particle A, 
taken to be an antiproton, as M and the mass of particle B, taken to be an electron 
or a proton, as m. Define electromagnetic form factors Fi(q 2 ) and ^(g 2 ), with 
normalization -F\(0) = 1 and 1*2(0) = K p — im v — 1, the anomalous magnetic moment 
of the proton, where q 2 = t in the t-channel case that we are solely interested in. 
Form factors are empirical quantities, obtained from experiment, which describe the 
fact that protons are not point-like particles and have an internal structure. They 
include all effects of the strong nuclear interaction inside the proton, hence are very 
difficult to calculate theoretically. The Sachs electric G E (t) = F 1 (t) + F 2 (t) tj (AM 2 ) 
and magnetic Gm(*) = Fi(i) + F 2 (t) form factors are used. In the t-channel, also 
known as the space-like region, the form factors are real functions of t. Although not 
treated in this thesis it is worth mentioning that this is not true in the s-channel, 
also known as the time-like region, where the form factors are complex functions of 
s. For a treatment of polarization observables in the time-like region see Ref. [80J. 
From Figure I3-1[ and using the Feynman rules for QED presented in Appendix |Dl 
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Figure 3-1: The electron, proton and antiproton electromagnetic currents; , Jp and Jp 
respectively. Time increases from left to right. These Feynman diagrams are converted into 
the mathematical expressions for the currents presented in this section using the Feynman 
rules of Appendix\I^ The shaded circle in the proton and antiproton currents describe that 
these are not point particles, and have an internal structure described by form factors. 

the electron current for particle B is 



where e is the electron charge and u(p2, X2) an d u (p±, A4) are the spinors of the 
incoming and outgoing electron respectively. 

Generalizing this to a non-point-like particle, such as the proton, with a finite 
extent defined by the Pauli and Dirac electromagnetic form factors F\(q 2 ) and i^g 2 ), 
one obtains the most general (anti)protorj*] electromagnetic current for particle A: 



where e p = — e p = — e is the charge on the proton and consequently the upper sign 

1 We are only interested in t-channel elastic scattering, i.e. not s-channel antiproton-proton 
annihilation which is suppressed in comparison to the i-channel amplitude in the low 1 1 | region of 
interest in a storage ring, as explained in Figure FT1 In this case an antiproton can be treated as 
a negatively charged proton, with internal electromagnetic structure described by the same form 
factors as the proton. Hence, as is customary in the phenomenology literature [551 EI], u and u 
spinors can be used for the antiproton current in the t-channel instead of the v and v anti-spinors 
required in the treatment of annihilation. The treatment presented here is in exact agreement 
with a treatment where the proton current involves spinors and the antiproton current involves 
anti-spinors. The minus signs introduced into eq. p. 2. lip because of the anti-spinor completeness 
relations of eqs. (|A.8l and lA.10j) only contribute to terms that vanish when the traces are evaluated 
using the trace theorems of eq. (|A.14|) . Hence all cross-sections and spin observables for i-channel 
elastic antiproton-antiproton, antiproton-proton and proton-proton scattering are equal to first 
order in QED. The anti-spinor formalism is less general then the one presented here as one is 
forced, due to the Gordon decomposition identity for anti-spinors derived in Appendix [Cl to specify 
that the particle and antiparticle have opposite anomalous magnetic moments. Hence the results 
could not be applied to antiproton-neutron scattering for example, as they can be in the present 
formalism. 




(3.1.3) 



Jjy = ±ie p u(p 3 ,\ 3 ) F^q 2 ) 7^ + 



%^ V { V3 -vx) v u(p 1 ,\ 1 ), (3.1.4) 



2M 
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is for the antiproton current and the lower sign is for the proton current, and where 



= ±[>f,>f] = 1(7*7" -7Y). (3.1-5) 



Using the Gordon decomposition identity [82, 83J : 



u (p') 7 M u (p) — u (p') 



.P + v' ) %o iXV (p' — p) v 



2M 2M 



u{p) , (3.1.6) 



which is derived in Appendix [C] the (anti)proton electromagnetic current can be 
written as 

JX = ±ie p u(p 3 ,Xs) (G M r ~ F 2 P \ + M P " ) u( Pl , AO , (3.1.7) 

which has a simpler gamma matrix structure than eq. (13.1.41) . hence will allow for 
easier computations. 

It is important to note that the structureless limit (point-like particles) is ob- 
tained when Fi(q 2 ) = 1 ancjfl F 2 (g 2 ) = 0, hence G7 A /(g 2 ) = ^(g 2 ) = 1- Applying 
this condition to the proton current in eq. ( I3.1.7P one obtains the electron current 
presented in eq. f)3.1.3p . Because we are interested in both antiproton-proton and 
antiproton-electron scattering, we shall calculate the generic case of the elastic scat- 
tering of two structured spin 1/2 particles. One can then take the one particle point- 
like limit to account for antiproton-electron scattering, or the two particle point-like 
limit to account for positron-electron scattering. Hence we generalize the electron 
current to account for non-point-like particles by using 

Jg = -ieu(p A ,\ A )(g Ml » - f2 P \ +P * )u(p 2 ,\ 2 ), (3.1.8) 

where we label the Dirac, Pauli and Sachs electromagnetic form factors of particle 
B by lowercase hit), f 2 (t), g M (t) = f\{t) + f 2 (t) and g E (t) = f x {t) + f 2 (t)t / (Am 2 ) 
respectively, with the usual normalizations /i(0) = 1 and f 2 {0) = k = ^ — 1 therefore 
9m{0) = At and ^(0) = 1. 

The above electromagnetic currents, and hence all results derived using them, 



2 More correctly for a point-like electron ^(0) = a/(2w) + 0(a 2 ) [84]. This is a result of the 
anomalous magnetic moment of the electron and was first calculated by Julian Schwinger shortly 
after the famous Shelter Island Conference in 1947 [5S]. The anomalous magnetic moment of the 
electron has now been calculated to order a 4 in QED and agrees with the experimentally measured 
value to 10 decimal places [55], making QED one of the most accurately verified theories in the 
history of physics. 
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can be applied to electrically neutral particles which nevertheless have an anoma- 
lous magnetic moment, such as the neutron, by encompassing the charge e into the 
definitions of the form factors and using the new normalizations -^"(O) = and 
F 2 n (0) = K n = —1.913 in units of nuclear magnetons. 



3.2 Generic elastic spin 1/2 - spin 1/2 calculation 

The helicity amplitudes can be related to the currents, to first order in QED, by the 
relation 

z^(A 3 A 4 ;A 1 A 2 ) = J/(M, A 3 , Ai) ( Z ^ SL J Jb{™, A 4 , A 2 ) . (3.2.1) 

Therefore the generic electromagnetic amplitude for a structured spin 1/2 particle of 
mass M scattering elastically off a structured spin 1/2 particle of mass m, via single 
t-channel photon exchange, is 



iM(X 3 X 4 ;X 1 X 2 ) = ±(-tYe 2 u(p 3 ,X 3 )[G Ml » - ^ R» )u(p 1 ,X 1 ) (3.2.2) 



T]fll> \ -I \ \ { v /' 



x( ^ ) u(p 4 , A 4 ) ( gul v ~ TT~r v )u(p 2 ,X 



q J \ 2m 



u(P3, A 3 )(g'm7 m - -^T7 RI *) u (Pu Ai) u(p 4 , X^)(gMlit ~ -fi-rJufa, A 2 ) 



g ^ V 2 M J \ 2m 



' v ' ' v ' 

where the upper sign is for p p oi pp scattering and the lower sign is for pp scattering, 
and where we have defined 

R" = pf + ^ and r v ee p£ + p v A , (3.2.3) 

and used the fact that = — e p = e. Note that in the t-channel the hermiticity of 
the electromagnetic currents implies that the form factors are real functions of q 2 = t 
{i.e. F x = F*, h = /*, F 2 = F 2 \ f 2 = f* and hence G M = G* M and g M = g* M ) 0. 
Now one must obtain 

4 4 

\M\ 2 = MM* = ^—r (AB) (AB)* = — A ABB*A\ (3.2.4) 

g 4 g 4 
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Figure 3-2: The Feynman diagrams for single photon exchange in antiproton-proton, 
antiproton- electron and positron- electron scattering respectively. As elsewhere in the thesis 
the time axis increases from left to right. The antiproton-proton case is generic and encom- 
passes the other two cases. In this section we calculate all spin- dependent phenomena for 
electromagnetic spin 1/2 - spin 1/2 elastic scattering, of particles with internal structure 
described by form factors, to first order in QED. 



where, using relations in Appendix [AJ the complex conjugates of A and B, defined 
in eq. f !3.2.2j> . are 



A* = u{ Vi Ai)(g m1 » ~ ^R^u(p 3 ,X 3 ), (3.2.5) 

B* = u{p 2 A2)\gMlv - -p- r v ) w(p4, A4) • (3.2.6) 
\ 2m J 

Substituting the above into eq. (13.2. 4p . requires using the completeness relations: 

u(pi,\i)u(pi,\i) = -(^p. + mi^(l + j 5 $i), (3.2.7) 

where A, and Si are the helicity and spin four vector of the particle with momentum 
Pi where i G {1, 2, 3, 4}, normalized such that Si^ = — 1 and constrained by the 
orthogonality condition Si^ = 0. The mass of the particle with momentum pi is 
denoted by m ; where we have that mi = m 3 = M and m 2 = m 4 = m. 

The result is a generic equation for all polarization phenomena in elastic spin 1/2 
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spin 1/2 electromagnetic scattering to first order in QED: 

16 ( f ) 4 1 1 2 = (3 ' 2,8) 



Tr 



Tr 



(f A + m) (1 + 75 &) [9m1 v - ^ J [f 2 + m) (l + 75 &) (<7m7" 



2m 

F<) R,, \ ( , A a. \ I „ F 2 R u 

2M 



i+ M ) (1 + 75 &) (Gm7m - ) (^ 3 + M) (1 + 75 &) (Gm7, " 



This generic equation can thus be used to calculate all helicity amplitudes and 
spin observables by substituting specific values for the spin (Si) and momenta (pj) 
four vectors, and can describe equal particle scattering in the case /1 — > Fi, f 2 — > F 2 
and m — > M. This also applies to antiproton-electron scattering by setting one 
particle to be point-like using fi — > 1 and f 2 — > and hence gu — > 1, in which case 
the first trace simplifies to the familiar electron trace: 



Tr 



+ m) (1 + 75 7" (V 2 + mj (1 + 75 $2) 7 M ] , (3.2.9) 



which when polarization effects are averaged over gives the familiar spin-averaged 
electron trace 



Tr 



[f 2 + m) 7m (f> 4 + m) 7,] =4 U 2/ ,p 4 , + p 2 ,P4^ + ^""J , (3.2.10) 



which has been evaluated using the trace theorems presented in Appendix |Aj 

Equation (I3.2.8P can be generalized to directly evaluate the squares of all elec- 
tromagnetic helicity amplitudes. We introduce the constants e, multiplying each Si, 
where % G {1,2,3,4}. 

16 ( - Y \M | 2 = (3.2.11) 



Tr 



Tr 



f,+ m){i + e 4 75 $,){g M i v - f j£)(t 2 + + $ 2 )(g M r - f -£ 

(^+M)(1 + e l75 $i)(g m ^ - f5 i )(^3+ M )( 1 + £ 3 75^)^7, " \ 



x 



Later we will set 6j = ±1 to account for different helicity states. Equation (13.2. lip is 
used to derive the helicity amplitudes and spin observables throughout the remainder 
of this chapter. 

The spin four vectors are now normalized so that all 6j = +1 corresponds to the 
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helicity amplitude 4>\ = + ; +, + ), and the ±1 in the helicity amplitudes now 

relate to the signs of the 

The momenta and longitudinal, transversal and normal spin four vectors in the 
Centre-of-Mass frame are presented in Table 13.11 Substituting combinations of these 
four vectors into eq. (13.2. lip , and computing the traces, will provide expressions for 
all spin-dependent cross-sections for elastic spin 1/2 - spin 1/2 scattering to first 
order in QED. The traces were computed using the computer algebraic program 
Mathematica and its add on package Tracer (87]. The Mathematica code for a 
typical calculation is presented in Appendix [Ej 



Centre-of-Mass Momenta vectors 


pi = (E A , 0, 0, k) 
p 2 = (E B ,0,0,-k) 


P3 = (Ea, k sin 9, 0, k cos 9) 
p 4 = (Eb, —k sin 9, 0, — k cos9) 


Centre-of-Mass Normal spin vectors 


S? = (0,0,1,0) 
S 2 N = (0,0,1,0) 


s£ = (0,0,1,0) 
S? = (0,0,1,0) 


Centre-of-Mass Transverse spin vectors 


Sf = (0,1,0,0) 
S 2 T = (0,1,0,0) 


^ 3 T = (0, cos^, 0, -sin0) 
S 4 T = (0, -cosfl, 0, sin^) 


Centre-of-Mass Longitudinal spin vectors 


S 2 L = — (— k, 0, 0, Eg ) 
m 


S 3 L = ( k, Ea sin 9, 0, Ea cos 9 ) 

SP = — (- k, E B sin 9, 0, E B cos 9 ) 
m 



Table 3.1: Momenta and spin ^-vectors in the Centre-of-Mass frame. The Centre-of-Mass 
energies of particles A and B are Ea = y/k 2 + M 2 and E B = y/k 2 + m 2 respectively, 
where k is the modulus of the Centre-of-Mass 3-momentum. The Centre-of-Mass scattering 
angle is denoted by 0. 



The spin 4- vectors of Table 13.11 satisfy the general expression: 

Ml E + M 



(3.2.12) 



3 Transverse to the direction of motion but still in the scattering plane, sometimes called Sideways 
and denoted S. 
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where p, E and M are the momentum 3-vector, the energy and the mass of the 
particle in question and s is a unit 3-vector identifying a generic spatial direction 
[SSI EH] ■ When s is parallel (or antiparallel) to p then the particle is longitudinally 
polarized while if § is perpendicular to p then the particle is transversely polarized. 



From this equation one can easily verify that S^S, 



•1 and p^S n = 0, where 



= (E , p), as is seen to be satisfied by all of the vectors in Table I3~T1 



3.3 Spin-averaged cross-section 

The spin- averaged differential cross-section for t-channel elastic spin 1/2 - spin 1/2 
scattering can be obtained by setting each Si = in eq. (13.2. lip , and multiplying 
by 2 4 = 16 to counter the four factors of 1/2 from the spin-dependent completeness 
relations that are absent in the spin-averaged completeness relations. One obtains 



M 



q 



Tr 



2 m 



ti + m ) ( 9 Ml 



f 2 r» 



2 m 



Tr 



{A + M ){ G Ml,-^){h + M ){ G Ml 



F2 Rv 



(3.3.1) 



These traces can be evaluated using the trace theorems of Appendix [A] to obtain 

F2 R p \ ( , , „ r \ I ^ E 2 R U 



Tr 



^ + M)(«^-ff)(^ + M)(« M> -f^) 



(3.3.2) 



4G^[p lfl p 3v +p lu p 3tl + -rf v ) -AG M F 2 R,R U + 2F 2 2 R,R U [ 1- 



t 



AM- 



The result follows from eq. (I3.3.ip . after including a factor of 1/4 from averaging 
over initial spin states and summing over final spin states, as presented in our recent 
paper [90] : 



s da 



Am 2 g 2 - tg 2 M \ ( AM 2 G\ - tG 2 M 



Am 2 — t 



AM 2 - t 



M'< 



m 



2\2 



s u 



2mMg E G E \ 1 2 2 
^ / 2 ^ M ' 



(3.3.3) 
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where we have used the electromagnetic coupling constant (fine structure constant) 
a = e 2 /4 7r and the Sachs electric form factor of each particle: 

G E (q 2 ) = F^ + ^Ftf) and g E (q 2 ) = f x (q 2 ) + -^—^ f 2 (q 2 ) . (3.3.4) 

Equation (I3.3.3p . describing the spin-averaged elastic scattering of any two spin 1/2 
particles or antiparticles, is a generalization of the famous Rosenbluth formula for 
elastic electron-proton scattering [91] . It corresponds to a twofold generalization of 
the Rosenbluth formula in that the mass of neither particle has been neglected and 
the internal structure of both particles is included, and is expressed here in a new 
invariant form. 



3.4 Helicity amplitudes 

Helicity is defined as the projection of the particles spin 3-vector in the direction 
of its momentum 3-vector. Helicity is a discretized quantity, having values of either 
±h j 2 for a spin 1/2 particle, because the spin of a particle with respect to an axis 
is quantized. Helicity states are always longitudinally polarized, i.e. either along 
direction of motion, which we denote by + for + 1/2, or opposite direction of motion, 
which we denote by — for —1/2 in spin 1/2 - spin 1/2 scattering^. 

The helicity amplitudes for 2 particle — > 2 particle scattering processes were 
introduced by Jacob and Wick [92J, and are represented by 

M. ( scattered, recoil; beam, target) = M. (A 3 , A 4 ; Ai, A 2 ) . (3.4.1) 

The arguments are to be read from right to left, as Ai and \ 2 correspond to the 
incoming particles in the reaction and A3 and A4 correspond to the outgoing par- 
ticles in the reaction. For spin 1/2 - spin 1/2 scattering the helicities \ — ± for 
i G {1,2,3,4} are + if the particles spin vector points in the direction of its mo- 
mentum vector and — if the particles spin vector points in the opposite direction 
to its momentum vector. The ± in the helicity amplitudes are shorthand for ±1/2 , 
the helicity of a spin 1/2 particle. For 2 particle — > 2 particle spin 1/2 - spin 1/2 

4 Remembering we set h = 1 throughout the thesis. 
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scattering there are 16 helicity amplitudes: 
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(3.4.2) 



Parity Invariance and Time-reversal Invariance, which are strictly satisfied in both 
electromagnetic and hadronic reactions, provide relations between these 16 helicity 
amplitudes. Parity changes the direction of motion but does not change the spin, 
hence it flips the helicity of a particle. Therefore Parity Invariance acts on the helicity 
amplitudes as followf] [921 [93]: 

M ( -A 3 , -A 4 ; -A 1; -A 2 ) = (-1) M ( A 3 , A 4 ; X 1 , A 2 ) , (3.4.3) 

where A = Ai — A 2 and fi = A3 — A4 , keeping in mind that each Aj is ±1/2 . Thus 
this reduces the 16 helicity amplitudes to eight independent ones: 
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(3.4.4) 



Time-reversal Invariance means that the amplitude for a reaction is unchanged if 

5 In the general case there would be a quantity 77 = ^ ^° (— jJ^+sb-sc-sd^ where s and 77 
correspond to the spin and intrinsic parity of each particle in the reaction, as a factor on the right 
hand side of eq. (|3.4.3[) . But for elastic scattering, which we are solely interested in, particles A 
and C are the same and particles B and D are the same; therefore 77 = (—1)° = 1. 
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the direction of time is reversed. It acts on the helicity amplitudes by interchanging 
incoming and outgoing particles as follows [92j [93] : 



M(\i, A 2 ; A3, A 4 



-1 
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(3.4.5) 



This reduces the eight remaining helicity amplitudes to six independent ones, by the 
relations, 
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(3.4.6) 



Hence there are six independent helicity amplitudes, 0i, 2 , 3 , 4 , 5 and 6 , for 
2 particle — >■ 2 particle spin 1/2 - spin 1/2 scattering, as follows: 
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(3.4.7) 



As can be seen from eq. (13.4. 7p . <pi and (f) 3 are non-spin-flip amplitudes, 2 and </>4 
are double-spin-flip amplitudes and 5 and 6 are single-spin-flip amplitudes. By 
double-spin-flip we mean that both particles in the reaction have their spins flipped, 
i.e. their helicities reversed, and by single-spin-flip only one particle in the reaction 
has its spin flipped, while non-spin-flip means none of the particles in the reaction 
have their spins flipped. The non-spin-flip amplitudes, 0i and 03, are also known as 
spin- elastic amplitudes. 

Note that for identical particle scattering one also has that 06 = — 05 an d hence 
there are only five independent helicity amplitudes in this will be shown. 

The multiplicity of an independent helicity amplitude is defined as the number of 
times it appears in the set of total helicity amplitudes, as seen in eq. (13.4.71) . Hence 
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0i, 02, 03 an d 04 each have multiplicity 2 while 05 and 06 have multiplicity 4. Thus 
one can write the spin-averaged differential cross-section in terms of the independent 
helicity amplitudes as 



da 



~T2 Yl 7 l^( A 3- A 4; Ai, A 2 



^ Ai A2 A 3 A4 



1 



2(8t0 



(|0i| 2 + |0 2 | 2 + |0s| 2 + |0 4 | 2 + 2|0 5 | 2 + 2|0 6 | 2 ). (3.4.8) 



We now calculate the modulus of all six independent helicity amplitudes for elastic 
scattering of spin 1/2 - spin 1/2 particles or antiparticles to first order in QED, 
by substituting the four longitudinal spin vectors (four pure helicity states) from 
Table I3TTI into eq. (13.2. lip and setting specific values (±1) for the e's to obtain 
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(3.4.9) 
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in agreement with those found by other methods [941 1951 196] . The linear combinations 
+ = 0! + 03 and 0_ = 0i — 03 appear often in the observables 



>1 + 03 

±a 



9m G m + 1 + 



t 



4/c 1 



5 _ m 2 _ j\//" 2 

2/iFi + f 2 F 2 ) , (3.4.10) 



— 9 m G 



m 



(3.4.11) 



along with the simpler combinations 2 + 04 = and 2 — 04 = 2 2 
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On comparison to previously published proton-proton helicity amplitudes [9U [951 
196] . and using the fact that ep = —e p , one finds that the ±a factors above are plus 
for like-charged particles (e.g. proton-proton or antiproton-antiproton scattering) and 
minus for unlike-charged particles (e.g. antiproton-proton scattering). 

The above t-channel helicity amplitudes can be transformed into the correspond- 
ing s- and w-channel helicity amplitudes using crossing symmetry [3U], I9"7]. 

3.5 Spin observables 

The spin observables for a 2 particle — > 2 particle elastic reaction, as introduced in 
[94], 196] , are described in Table 13.21 In this section we present expressions for them 
to first order in QED for t-channel electromagnetic scattering. In particular the spin 
transfer and depolarization observables will play an important role in the remainder 
of the thesis. 

3.5.1 Polarization transfer observables 

Setting Si = S4 = in eq. (13.2. lip and subtracting the spin-averaged differential 
cross-section gives a generic equation for spin transfer from initial particle B to final 
antiparticle A : 

= P%Tls { F > t «' S * [^fi P3 .S 2 + f 2 (2(mi + M*-s) q .S 2 
ail 8m Ms t l 1 L v v ' 

+ p 3 ■ S 2 t - q ■ S 2 t)] - lQm 2 M 2 Fi g E Pi ■ S 2 q ■ S 3 
+ AM 2 Ge [4m 2 fi (p 3 -S 2 q-S 3 - S 2 -S 3 t) 

+ t f 2 (p 3 ■ S 2 q ■ S 3 + q- S 2 (p 2 + Pi) ■ S 3 - S 2 - S 3 t)}} , (3.5.1) 

into which specific vectors S2 and S3 will be inserted to give the various polariza- 
tion transfer observables. Scattering is in the XZ plane, so the coordinates are 
X (Transverse), Y (Normal) and Z (Longitudinal), where in the above equation 

i,je{X,Y,z}. 
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Unpolarized: 




1 particle polarized: 



Single spin asymmetries 




2 particles polarized: 



Polarization transfer: 



Depolarization: 



Double spin asymmetries: 



3 particles polarized: 



4 particles polarized: 




Table 3.2: The 16 spin observables of a 2 particle —> 2 particle elastic reaction. An 
upward pointing triangle represents that the beam is polarized, while the absence of a triangle 
represents an unpolarized beam. The right hand column shows the symbols for the spin 
observables that will be used throughout the thesis, as defined in Ref. \9$ . The subscripts 
i,j,k,l G {X,Y,Z} correspond to the direction of the polarization of each particle. Time 
increases from left to right as elsewhere in the thesis. 
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The spin observables Kxx, Kyy an d K Z z 

The spin observable Kxx is obtained by inserting the transverse polarized spin vectors 
£2 = S 2 and S3 = S 3 T from Table 1331 into the generic spin transfer equation (13. 5. lft : 



. a 
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xx 
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{4mVi (M 2 F X - k 2 F 2 ) + f 2 [-U?M*Fr 



+ (4A; 4 + {Ah 2 + t) Vk 2 + mWk 2 + M 2 ) F 2 } . (3.5.2) 

Inserting the normal polarized spin vectors S 2 = S 2 and S3 = S 3 N from Table 13.11 
into the generic spin transfer equation (13.5.11) gives 
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(3.5.3) 



Inserting the longitudinally polarized spin vectors S 2 = S 2 and S3 = S 3 L from 
Table 13.11 into the generic spin transfer equation (13.5. ip gives 
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The spin observables Kxz and Kzx 

When the spin four vectors S 2 = S 2 and S3 = S 3 from Table 13.11 are substituted 
into the generic spin transfer equation (13.5.11) we obtain 
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Inserting the spin four vectors S2 = 5* 2 and 5*3 = S 3 from Table 13.11 we obtain 
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As expected by parity and time-reversal invariance we confirm that 

K XY = K YX = K YZ = Kzy — . (3.5.7) 

3.5.2 Depolarization spin observables 

The observable Dij is the polarization remaining after interaction with the target, 
while of interest here is the loss of polarization after interaction with the target, 
i.e. (1 — D^). We present results to leading order in small 1 1 |. Here setting 5 2 = 
S4 = in the generic eq. (13.2.1 ip and subtracting the spin-averaged differential cross- 
section, then subtracting the result from the spin-averaged equation gives a generic 
depolarization equation for antiparticle A, into which the various vectors Si and 
S3 can be substituted. The Mathematica code to derive this generic depolarization 
equation is included in Appendix [Ej but the equation itself is too long to be practical 
to include here. The ~ sign means to first order in small \t\. 

The spin observables ( 1 — Dxx ) ? ( 1 — -Dyy ) and ( 1 — Dzz ) 

Substituting the transverse spin vectors Si = S± and 5*3 = S 3 T from Table 13.11 into 
the generic depolarization equation gives 

- k 2 m' (s + M 2 - m 2 ) /, / 2 ]} . 

Inserting the normal polarized spin vectors S\ = Si and S3 = S^ from Table 13.11 
into the generic depolarization equation gives 

d <j o? 

—— ( 1 — D YY ) = — gl 4 Glj , complete to all orders in t. (3.5.9) 

CI u ' L. "S 

Inserting the longitudinal spin vectors S\ = S± and S3 = S 3 L from Table I3TT1 into the 
generic depolarization equation gives 

^d-^zz) - Z*^ {[M *p +m > )jf + ak * fi (3.5.10) 
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The spin observables ( 1 — D X z ) an d ( 1 — D Z x ) 

Substituting the spin vectors Si = S± and 5*3 = S 3 T from Table 13.11 into the generic 
depolarization equation gives 

« ^ j JJ" [2k* + m 2 M*+ e{s-2e)], (3.5.11) 

d (X 

and -j-^ ( 1 — Dzx ), found by inserting into the generic depolarization equation the 
spin vectors Si = S^ and S3 = S 3 L from Table I3.1[ is the same as above to leading 
order in t. Note this is just the leading t part of the spin-averaged case. 
As expected by parity and time-reversal invariance we confirm that 

Dxy = D YX = D YZ = Dzy = . (3.5.12) 



3.5.3 Spin asymmetries 

For electromagnetic scattering, to first order in QED, all single and triple spin asym- 
metries are zero, 

Ai = C ijk = where i, j, k G{X,Y,Z}, (3.5.13) 
and the double spin asymmetries equal the polarization transfer spin observables: 

Aij = Kij where i, j G { X, Y, Z } . (3.5.14) 



relicity amplitudes all being real, 
There are also 'four spin mea- 



as a consequence of the tree-level electromagnetic 
and the fact that 02 = — $4 to first order in QED@. 
surement' spin observables C^i as described in Table 13.21 where the polarization of 
the beam, target, scattered and recoil particles are measured. These 'four spin mea- 
surement' spin observables are not needed for spin filtering as the polarization of the 
recoil particle will not be measured, hence they will not be treated in this thesis. For 
a treatment of this see Refs. [96j |98j [99] . 

Radiative corrections are not treated in this thesis. In elastic e~ p scattering 
radiative corrections are estimated to give corrections of 1% - 3%, and are only 
considered in very high precision experiments. This is treated in Refs. |1UU[ I101[ 

[1021 rrna] . 



6 See expressions for the spin observables in terms of helicity amplitudes presented Ref. [104] and 
in Table A. 10.5 of Ref. EHJ. 
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Chapter 4 



Specific helicity amplitudes and 
spin observables 

"Mathematics, rightly viewed, posses not 
only truth, but a supreme beauty - a beauty 
cold and austere, like that of a sculpture. " 

Bertrand Russell 

The expressions for the generic spin 1/2 - spin 1/2 electromagnetic helicity am- 
plitudes and spin observables derived in Chapter [3] are presented in this chapter for 
the specific cases of: antiproton-proton scattering in section 14.2} antiproton-electron 
scattering in section 14.31 and positron-electron scattering in section 14.41 Then the 
specific cross-sections and spin observables needed for spin filtering are explicitly 
presented in section 14.51 These results are of importance to many areas of particle 
physics, and will be utilized throughout the remainder of this thesis. The chapter 
concludes with a calculation of all spin - spin 1 helicity amplitudes in section I4.6[ 
which describe the scattering of deuterons off a carbon nucleus for example. 
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4.1 Introduction 



In this chapter we calculate all spin-averaged and spin-dependent electromagnetic 
cross-sections for elastic antiproton-proton, antiproton-electron and positron-electron 
scattering, to first order in QED. These are required in the descriptions of spin 
filtering that follow later in the thesis. In the region of low momentum transfer of 
importance in storage rings electromagnetic effects dominate over hadronic effects. 
So we calculate the electromagnetic contribution to these cross-sections, and focus on 
the low 1 1 1 behaviour. As explained in section 17.2.51 there is a minimum momentum 
transfer q m m (and correspondingly a minimum scattering angle m i n ), below which 
scattering is prevented by Coulomb screening, corresponding to an impact parameter 
of the Bohr radius as = 52900 fm of the atom [57]. Antiprotons flying past an atom 
at impact parameters greater than clb see the atom as neutral and do not interact with 
the atom [72]. For momentum transfers of q > q m \ n = 1/ob the energy transfered is 
greater than the binding energy of the hydrogen atom, hence to a good approximation 
the antiproton scatters from free protons and electrons in the hydrogen atom [68, 69J, 




(4.1.1) 



where H represents a hydrogen atom. 

Since the objective is to polarize the antiproton beam the transfer of polarization 
from initial target particles to the final antiproton beam is of principal importance. 
Hence a thorough investigation of the elastic polarization transfer reactions 

p + p r — > p t + p (4.1.2) 
p + e _t — > p T + e" (4.1.3) 

is presented. The effects of depolarization on the antiproton beam are also calculated. 

We are also interested in the interactions of an antiproton beam with an opposing 
polarized electron beam, which eq. ( I4.1.3P equally describes. Coulomb screening 
also introduces a minimum scattering angle, 9 m i n , in this case; due to the impact 
parameter being limited by half of the average separation of electrons in the beam, 
as explained in section 17.2.51 
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4.2 pp helicity amplitudes and spin observables 



In this section we present the expressions for the electromagnetic helicity amplitudes 
and spin observables in the specific case of antiproton-proton scattering. We have 
only derived expressions for t-channel scattering, while for antiproton-proton scat- 
tering there is also a contribution from s-channel (annihilation) scattering, as shown 



in Figure 14-11 However the t-channel results are dominant in the low 1 1 1 region of 
interest in a storage ring, as explained in Figure 14- 1[ and electromagnetic effects also 
dominate over the hadronic effects in this low 1 1 | region. The results in this section 
are important in the region where \ t\ < \t c \ for antiproton-proton collisions with 
total cross section of ^, defined by [95, 105J 



1 7r a 



Aab C 



0.001 (GeV/c) 



(4.2.1) 



where the electromagnetic interaction dominates the hadronic interaction, as derived 
in section [7.4.11 Here the laboratory velocity is /3j a b = ^ s(s — AM 2 ) / (s — 2 M 2 ) 
and p = TZe{<j)+}/Tm{(f)+} the ratio of real to imaginary parts of the hadronic 
non-flip amplitudqj. 

The electromagnetic helicity amplitudes and spin observables for pp scattering 
can be obtained by setting equal masses and form factors (/i = F%, f 2 = F 2 , Qe = Ge, 
qm = Gm and m = M) in the expressions provided in sections (13.41 and l3.5p . These 
are required by the VAX Collaboration to analyze the buildup of polarization of an 
antiproton beam by interactions with the protons in a hydrogen target. The helicity 



amplitudes in section 13.41 now become 
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1 The total pp cross-section behaves like of t ~ 75.5//?i a b mb for energies up to 1 GeV [5511134) . 
and p 2 ~ at all energies measured thus far (see Figure 4-12 of Ref. |140j ). Therefore t c is energy 
independent, at least up to energies of 1 GeV. 
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where 6 = — 5 in this expected. The generic spin transfer equation for 

this case is 



da 



a 2 G 2 M 



dfl"* " 8M , a1 ,U-S, q -S,{4M>-2s-t)tF> (4.2.3) 

+ A[2p 3 -S 2 q-S 3 + q-S 2 (p 2 + Pi ) ■ S 3 ] M 2 t F 2 G E 

+ 16M 4 G E (q-S 2 q-S 3 - S 2 -S 3 tG E )} , 



and the spin transfer observables now become 
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For pp elastic scattering Kxz = Kzx, and we obtain 
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The depolarization spin observables to leading order in small 1 1 | become 
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da. ^ 



da. ^ , da 4a 2 F 1 4 , , 9 „,.2 



Note that to first order in small 1 1 | the antiproton and proton form factors can be 
approximated as F 1 m 1, F 2 ~ /i p — 1 (i.e. Gm ~ A*p and G E ~ 1), where // p is the 
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magnetic moment of the proton. 

The t-channel spin-averaged differential cross-section for this case simplifies to 



.s da (AM 2 G 2 -tG 2 M \ 2 (su\ (2M 2 G 2 \ 2 1 , 



The above helicity amplitudes and spin observables are seen to be correct by a 
comparison to the relations between the helicity amplitudes and spin observables 
presented in Ref. |104] and in Table A. 10. 5 of Ref. [98] . 

In addition to antiproton-proton scattering the results of this section equally 
apply to any spin 1/2 (anti)baryon-(anti)baryon elastic electromagnetic scattering 
in the t-channel. Also, given that quarks are spin 1/2 particles, the above expressions 
could be applied to t-channel (anti) quark- (anti) quark electromagnetic scattering if 
in future quarks are found to have an internal structure. 



4.3 pe helicity amplitudes and spin observables 

The spin observables for p e scattering can be obtained by setting f\ = 1 and /2 = 
(i.e. gE = Qm = 1) in the expressions from sections (13.41 and l3.5[) . to account for the 
elastic scattering of a structured particle off a point-like particle. These are required 
by the VAX Collaboration to analyze the buildup of polarization of an antipro- 
ton beam by interactions with the electrons in a hydrogen target, or interactions 
with a polarized electron beam as treated in Chapter [7J The helicity amplitudes in 
section I3T41 for antiproton-electron elastic collisions become 

^ Ak 2 ) t 1 2 ' 

m F 2 
~ 2 M ' 

i (s -m 2 + M 2 )Fi 

4 k 2 s J ' 

M (s + m 2 -M 2 )Fi F 2 ' 
4 k 2 s 2M ' 
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The generic spin transfer equation for this case is 



dff -2a 2 mM G M 

dQ Kij ~ 7t 



F 2 F ]_ j 
G E S 2 ■ S 3 + ( YJ^p Pa- S 2~ —q - S 2 )q- S 3 



, (4.3.2) 



which is a generalization of equation (3) of Ref. [68J, which was at the root of the 
initial interpretation of the FILTEX results on spin filtering in 1994. We do not 
neglect the 1/M 2 term, or make the non-relativistic approximations of eqs. (I4.3.6p . 
as done in Ref. [68J. 

The spin transfer observables for antiproton-electron elastic scattering are 
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The depolarization spin observables to leading order in small 1 1 | for this case are 



da 

cm 

dff 
dm 

dff 
cm 

dff 

dm 



1 — D 



xx 



-m 2 a 2 F/ , 2 „^ 2 \2 

t_ ( e i hjt * i 

2 k 2 s 2 t 



1-Dyy) 



a 



r 2 



1-D 



zz 



(1-D 



2s 

—M 2 a 2 F 2 
2 k 2 s 2 t 

dff 



(s-m 2 + M 2 )' 
complete to all orders in t . 
( s + m 2 -M 2 ) 2 , 



(4.3.4) 



xz 



dtt 



1-D 



zx 



dff 4 a 2 F? , , 9 9 , , 9 \ 
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The spin-averaged differential cross-section from section 13.31 for this case simplifies 
to 



s da / 4 M 2 G\ — t G\ { \ (M 2 -m 2 ) 2 - su (2mMG E V 1 2 

^dti ~ V 4MM J ^ + V t ) + 2 ' 

(4.3.5) 

the familiar Rosenbluth formula [91] , where we have not neglected the mass of the 
electron. 

The low energy (non-relativistic) limit of the above equations can be obtained 
using 



■ E cm + £cm) 2 ~ ( m + M) 2 « M ! 



Ge ~ 1 , (4.3.6) 
Rj // p = 1 + k p . 

Of particular importance is the non-relativistic limit of (da/dfl) Kyy in eq. (I4.3.3[) 
for polarization transfer in pe~^ — > e~ scattering, which using t = —Ak 2 sin 2 | 
where 9 is the Centre-of-Mass scattering angle gives 

da a 2 m (1 + k p ) 

dfi ^ YY ~ " 2PM sin 2 (f) ' (43 ' 7) 

This is precisely equation (4) of Ref . [68] , used extensively throughout their paper and 
early VAX calculations. The equations presented in this thesis generalize this work 
to relativistic energies. Non-relativistic expressions for the spin transfer observables 
have also been presented recently in Ref. [SSI EH], in the context of a proposal to 
polarize antiprotons by repeated interaction with a co-moving polarized positron 
beam at very low relative velocities [75] . 

In addition to antiproton-electron elastic scattering the results of this section 
apply also to any (anti)baryon-(anti)lepton elastic scattering. In particular they 
apply to antiproton-muon elastic scattering which may be of use in spin filtering, 
where the fact that muons have approximately 200 times the mass of electrons will 
greatly enhance the polarization transfer cross-sections Kxx and Kyy in eq. (I4.3.3p . 
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4.4 e + e helicity amplitudes and spin observables 



The electromagnetic helicity amplitudes and spin observables for t-channel elastic 
positron-electron scattering can be obtained by making the transformations f% = 
Fx — >• 1, / 2 = F% — > 0, qe = Ge —> 1, 9m = Gm — > 1 an d M — > m in the expressions 
provided in sections (13.41 and 13. 5p . The helicity amplitudes from section I3~41 now 
become 



s + 4 k 2 m 
+ 



— a 



2t 



2 k 2 



_2_ 

■ a 

- a 



— a 



m 
2k 2 

2 m 2 



t 



1 + 



Ak 2 



(4.4.1) 



05 

— a 



m 
Ik 2 



s (4 k 2 + t) 



where again 6 = — 5 in this case, as was found for pp scattering. 
The generic spin transfer equation for this case is 



— K 



2a 2 m 2 
s t 



5*2 ' 5*3 ~ 



q- S 2 g- S 3 
t 



(4.4.2) 



and the spin transfer observables for e + e elastic scattering are 
da 



da 

dn 

da 



K 



xx 



A" 



YY 
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a m 
2sA; 2 ' 

2a 2 m 2 
st 

a 2 



K 77 = 

dtt ZZ 8 k 2 s t 



[s(4k 2 + t) + 16 fc 4 ] 



For e + e elastic scattering Kxz = Kzx hence one obtains 



(4.4.3) 
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The depolarization spin observables to leading order in small 1 t 1 for this case are 



da 



2a 2 m 2 
k 2 st 



1 — D 



xx 



(k 2 + m 2 ) 



a 



— complete to all orders in t . 



da 
dU 

da 
dti 



— 2a 2 m 2 
k 2 st 

da 



1 — D 



zz 



1 — D 



(k 2 + m 2 ) 



(4.4.4) 



xz 



1 — D 



zx 



The t-channel spin-averaged differential cross-section for this case simplifies to 



s da 

a* cm 



-s u 



2 m 



2 \ 2 



+ r 



(4.4.5) 



which can be written in the more familiar form 11081 



da 

cm 



a 

7i? 



(2 



m 



,2 t Z 
s) + St + - 



(4.4.6) 



In elastic electron-positron scattering, also known as Bhabha scattering [109J, for 
very low momentum transfer the t-channel part of the spin-averaged differential 
cross-section dominates over the s-channel part. But the complete spin- averaged 
Bhabha differential cross-section must include both t-channel and s-channel Feynman 
diagrams in the amplitude, as shown in Figure 14-11 to give 



t-channel part 



da 

cm 



a 

+ 



(2 



m 



a 



(2 



m 



s) + st + - 



t) 2 + St 



(4.4.7) 



a 

7t 



[(s + t) 2 - 4m 4 ] . 



s-channel part 



cross-term part 



where one sees that the t-channel part transforms into the s-channel part, and vice- 
versa, by the well known crossing symmetry obtained by interchanging s and t. 

The results of this section also apply to any combination of t-channel elastic 
(anti)lepton-(anti)lepton, (anti)lepton-(anti)quark and (anti) quark- (anti) quark elec- 
tromagnetic scattering. This is particularly useful since the masses of each particle 
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M t oc - M s oc - 

t s 

Figure 4-1: The two tree-level Feynman diagrams for elastic positron-electron scattering, 
also known as Bhabha scattering. The antiproton-proton case is analogous. Contribu- 
tions to the full amplitude M = Mt + M s come from the t-channel (left diagram), being 
proportional to 1/ 1, and s-channel (right diagram), being proportional to 1/s, as shown. 
For low momentum transfer (small \ t \), and also at high energies (large s), the t-channel 
contributions dominate. As always time increases from left to right. 

have been retained in the equations. 



4.5 Observables needed for spin filtering 

In this section we present the leading t approximation of all spin observables needed 
for spin filtering. These expressions will be integrated over a range of t later in the 
thesis. All expressions are written in terms of invariants using eqs. (14.5.21 and 14.5.41) 
above. The low 1 1 | approximation of the form factors^, ~ 1 and F^it) ~ fj, p — 1 
therefore ~ 1 and Gm(^ ~ are used, as is seen to be valid in the dipole 

model for the Sachs form factors 

G E (t) = ^ = Sl (4.5.1) 

W \h (1-t/A 2 ) 2 1 ^ 

with A 2 = 0.71 (GeV/c) 2 obtained from a best fit to experimental data |110[ 1111] . 

Spin filtering is an azimuthally symmetrical process, as will be explained in Chap- 
ters and El hence the transverse contributions will be averaged in what follows. 



2 Writing F\(€) and F^ii) in terms of GE(t) and Gw(t) and then using a Taylor expansion of 
eq. (|2XH one obtains the Taylor expansions F^t) = 1 + 2.30777 t + 4.37246 t 2 + O (t 3 ) and 
F 2 (t) = 1.79285 + 5.5594 t + 12.2483 t 2 + O (t 3 ). 
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The Centre-of-Mass expressions presented throughout this chapter and the previous 
chapter can be transformed into invariants, or into the Laboratory reference frame 
using the following relations |112j : 

A = Ak 2 s = [s — (m + M) 2 ] [s — (m — M) 2 ] , (4.5.2) 
= s (s - AM 2 ) whenm = M, (4.5.3) 



da 7i da Atx s da 



dt k 2 dtt X dfl 



(4.5.4) 



M 

k = —= pi a b , (4.5.5) 



Q 

t = -2k 2 (1 -cosfl) = -4£; 2 sin 2 -, (4.5.6) 

where 9 is the Centre-of-Mass scattering angle and A is a Lorentz invariant. Of 
interest in eq. (I4.5.6P are the particular cases that t = when 6 = and t = —4 k 2 
when 8 = it, corresponding to total backward scattering. One sees that for elastic 
scattering t < for all Centre-of-Mass scattering angles 9 ^ 0. 



4.5.1 Antiproton - proton scattering 

In this section, where the masses of the two particles are equal, A is the invariant 
defined in eq. f)4.5.3p . The ~ sign refers to the first term in the expansion in t. To 
leading order in small 1 t 1 the relevant spin observables for single photon exchange 
antiproton-proton scattering are: 

^xx + ^yy da Ana 2 M 2 fi 2 
2 clt ~ AT ' 

d-^ + d-^ d, K zi^lMM [2 „., + M1 -„)]', 

/, n da -2na 2 (X + AM 2 s) r ~ W1 ,, 2 

The leading t approximation of the spin-averaged cross-section for this case is: 

da ^ W (s-2M 2 ) 2 

dt ~ A t 2 y } 
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4.5.2 Antiproton - electron scattering 

In this section, where the masses of the two particles are not equal, A is the invariant 
defined in eq. (I4.5.2p . The leading t terms of the relevant observables for antiproton- 
electron scattering are: 

K X x + K YY da Ana 2 m M fj, p 
2 cU ~ At ' 

(1- D xx ) + (1- D YY ) da -Airm 2 a 2 (s - m 2 + M 2 ) 2 



2 dt A 2 t 

zz dt ~ At ( s ~ m M ) ' ( 45 - 9 ) 
. d a -8 7T M 2 a 2 (s + m 2 - M 2 ) 2 

(1 " Dzz) dt * >Tt • 

The leading t approximation of the spin-averaged cross-section for this case is: 

do- Ana 2 ( s - m 2 - M 2 ) 2 

-t X —~ 13 (4 ' 5 ' 10) 

4.6 Spin - spin 1 helicity amplitudes 

In this section we calculate the electromagnetic helicity amplitudes for spin - spin 
1 scattering f|, generalizing the spin - spin 1 helicity amplitudes presented to leading 
order in the low 1 1 1 approximation in Ref. |113j . While not being directly used later 
in the thesis this calculation uses the formalism developed earlier, and would be 
applicable to the scattering of deuteronsj (which are spin 1) off a carbon nucleus 
(spin 0). 

A spin 1 particle has three possible spin states, represented by —1, and +1, 
while a spin particle has only one spin state. We represent the helicity amplitudes 
for spin - spin 1 scattering as M. ( A;, ; A a ) where A a and \ are the helicities of the 
incoming and outgoing spin 1 particle respectively, hence \ a , \ b e {—, 0, +}. Hence 



3 This is the only section of the thesis where the particles involved are not spin 1/2. 

4 The nucleus of a deuterium atom, a bound state consisting of a proton and a neutron, is called 



a deuteron. 
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the complete set of spin - spin 1 helicity amplitudes are: 



1. M(+;+) 
2- M(+;-) 
3. M(+;0) 



4. M(—,+) 

5. M(-;-) 

6. M(-;0) 



7. M(0;+) 

8. M(0;-) 

9. M(0;0) 



(4.6.1) 



All definitions from section 13.41 are equally valid here and one obtains A — /i = 
( Ao — A a ) — ( Ao — Xt ) = At — A a . Hence Parity Invariance defined in eq. (13.4.31) 
gives the following relations between the helicity amplitudes 



M(-;-) 
M(-;0) 
.M(0;-) 
M(-;+) 



(-1) 1 - 1 M(+;+) 
(-1) 1 - M(+;0) 
(_l)0-i M (0;+) 

(-l)l-(-l) A^(+; 



M(+;+) 
-M(+;0) 
-M(0;+) 
M(+;-) 



(4.6.2) 



Thus reducing the nine helicity amplitudes to five independent ones. Applying Time- 
reversal Invariance, defined in eq. (I3.4.5[) . further reduces to four independent helicity 
amplitudes by the following relations: 



M(+',0) 
M(-;0) 



-l) 1 - M(0;+) 
-l)- 1 - M(0;- 



-M(0;+) 
-M(0;-) 



(4.6.3) 



Now we can present the four independent helicity amplitudes for spin - spin 1 
scattering: 



H, = M(+;+) = M{-;-) 

H 2 ee M(+;0) = M{0;-) 

H 3 ee M(+;-) = M{-;+) 

H 4 ee M(0;0), 



M(-;0) 



(4.6.4) 



where it is seen that Hi and H 4 are the non-spin-flip amplitudes, while H2 and 
i?3 are the spin-flip amplitudes. The non-spin-flip amplitudes, Hi and if 4, are also 
known as spin- elastic amplitudes. It can be seen from eq. (14.6.41) that Hi and if 3 
have multiplicity 2, H 2 has multiplicity 4 and if 4 has multiplicity 1. Thus one can 
write the spin-averaged differential cross-section in terms of the independent helicity 
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amplitudes as 
da 



S dQ ~ {87r f {2sA + l)(2s B + l) lM{Xb;Xa 



' (2 1^1 2 + 4|# 2 | 2 + 2|# 3 | 2 + l^ 4 | 2 ) , (4.6.5) 



3 f8 7r) 2 



where sa = and Sb = 1 are the spins of the two particles in the elastic scattering 
process. The helicity amplitudes for elastic spin - spin 1 electromagnetic scattering 
can be found by calculating 

M(Xb ; A ) = J " JM(Ab ; Aa) , (4.6.6) 

where and J ^ are the spin and spin 1 electromagnetic currents respectively, as 
defined below. 

The spin current is very simple as it has no helicity structure. The electro- 
magnetic current for a spin particle of charge Z e and form factor Fo(q 2 ) is simply 

isaini] 

j, = ZeF (q 2 ) (p + pV (4.6.7) 

Since the deuteron is a spin 1 object its electromagnetic structure is described by 
three form factors, charge monopole Gc, charge quadrupole Gq and magnetic dipole 
Gm, assuming parity and time- reversal invariance. 

The most general form of the deuteron electromagnetic current, assuming Parity 
and Time-reversal invariance is |115L 11161 1117j : 



J"(X b ; K) = ee*(p b , X b 



e a ( Pa , X a ), (4.6.8) 



where M d is the mass of the deuteron, R^ 1 = Pa + p b and e CT (p a , A a ) is a po- 
larization 4-vector restricted to three independent components by the condition 
Pa e ti{Pa, X a ) = 0. The quantities F*, F 2 d and Gf are the electromagnetic form 
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factors of the deuteron with normalizations 



F x d (0) =1, 



F 2 d (0) = Q + m - 1, 



Qi(0) = fi d . 



(4.6.9) 



where Q is the quadrupole moment of the deuteron in units of e/Mj and /i^ is the 
magnetic dipole moment of the deuteron in units of e/2 Md . 

The complete set of initial and final deuteron polarization vectors, in the Centre- 
of-Mass (CM) frame, where the initial momentum of the deuteron isp a = ( E d , 0, 0, k 
where Ed = a/A; 2 + M d 2 is the energy of the deuteron, follows [TPS] : 

1 



Initial < 



e"(p.,0) 



0, -1, -i, 



— ( fe, o, o, £ d 



V2 

1 



7! (0, 



(4.6.10) 



Final < 



—j= (0, — cos 9, i, sin 9) 

— ( k, E d sin 0, 0, E d cos ' 

—j= (0, cos^, i, — sin 9) 
v2 



(4.6.11) 



where 6 1 is the CM scattering angle and i 



•1. The momentum transfer is 



0, — ksin9, 0, k — kcos9) . 



(4.6.12) 



Combining all of the above into eq. (I4.6.6P gives the helicity amplitudes for a deuteron 
of mass Md colliding with a spin zero nucleus of charge Z e and electromagnetic form 
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factor F (t), as follows: 



Ze 2 F 



(— ) I + is? ) + G ? 



1 + 



4k 2 J ' 



(4.6.13) 
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Mj k 2 
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2 1 



-f 2 /c 2 



, (4.6.14) 



Ze 2 F 



= f; 



4 k 2 



2 AM 2 



+ Gf 



1 + 



4P 



(4.6.15) 



H A 



Ze 2 F 



d (s-u) 



2 1 1 

t + mF + k 2 



+ F 9 d 



t(s-u) ( 1 



8M 2 



Mi 



1 

fc2 



+ 2G? 



s — u) t 

+ TTTT + 1 



AM 2 ' 4P 



(4.6.16) 



In the above expressions the dependence on the mass of the spin-0 nucleus is in the 
( s — u) terms. 
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Chapter 5 

Polarization buildup by spin 
filtering 

"If, as I have reason to believe, I have dis- 
integrated the nucleus of the atom, this is of 
greater significance than the war. " 

Ernest Rutherford, apologizing for absence 
from a meeting of the International Anti- 
submarine Warfare Committee. 

The theory of spin filtering is developed in this chapter. A mathematical descrip- 
tion of the related but simpler process of polarization buildup by the Sokolov-Ternov 
effect is first presented in section 15.11 The ideas presented are utilized in the mathe- 
matical descriptions of spin filtering which follow. In section 15.21 the rates of change 
of the number of particles in each spin state are combined into a set of polarization 
evolution equations which describe the process of polarization buildup by spin fil- 
tering. This set of polarization evolution equations is then analyzed and solved in 
section I5.3[ emphasizing the physical implications of the dynamics. 
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5.1 The Sokolov-Ternov effect 



The fact that an electron beam acquires a 'self-polarization' due to the emission 
of synchrotron radiation in a storage ring is called the Sokolov-Ternov effect [61], 
and has been described in section 12.3.31 of this thesis. It turns out that one can 
describe the Sokolov-Ternov effect by a system of polarization evolutions equations 
very similar to a scenario of spin filtering when there is no scattering out of the ring. 
The physical principles behind both systems are identical, the polarization buildup 
in both being induced by a discrepancy in the spin-flip transition rates. The spin-flip 
in the Sokolov-Ternov effect is induced by synchrotron radiation as a result of the 
charged particle being bent in a magnetic field, whereas in spin filtering the spin-flip 
is induced by interactions with a polarized internal target. 

In order to introduce the mathematics of systems of polarization evolution equa- 
tions we now present and solve a set of polarization evolution equations that describe 
the Sokolov-Ternov effect. This will provide a stepping-stone to the mathematical 
description of spin filtering which follows later in the chapter. The intermediate de- 
tails of the calculations will be presented here, but omitted in later sections. Please 
note as described in section 12.3.31 that this effect is much stronger for electrons than 
(anti)protons and the Sokolov-Ternov 'self-polarization' is not a practical method to 
produce a polarized antiproton beam at present energies. 

We denote the transition rates in this section by W a b instead of cr a &, where a, b G 
{+, — }, to distinguish from the spin filtering spin-flip case. It is important to note 
that there is no target in the Sokolov-Ternov effect, which makes it physically different 
to spin filtering, although it can be described similarly. The number of particles in 
the 'spin up' and 'spin down' states can change by two mechanisms while emitting 
synchrotron radiation: (1) 'spin up' particles can be flipped to 'spin down' particles, 
the cross-section for which we label as W+- ; and (2) 'spin down' particles can be 
flipped to 'spin up' particles, the cross-section for which we label as W_ + . Mechanism 
(1) constitutes a decrease in the number of 'spin up' particles (N + ) and an increase in 
the number of 'spin down' particles (iV_), while mechanism (2) constitutes a decrease 
in the number of 'spin down' particles and an increase in the number of 'spin up' 
particles. This explains the signs of the coefficients in eq. (j5.1.ip which follows. 
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Therefore the Sokolov-Ternov effect can be described by the polarization evolution 
equations 



d 


" N + ~ 




' -W+- 






' N + ~ 


d7 


. N - . 






-w_ + _ 




. N - . 



(5.1.1) 



where r is the time variable 3. 

As long as the two spin-flip transition rates W+- and W-+ are not equal, i.e. 
W + _ 7^ W_ + , there will be a buildup of beam polarization over time. It has been 
found that there is a slight difference between these rates, and this is the basis of the 
original Sokolov-Ternov idea [6T] . 

This system is identical to the system presented later in eq. (I5.2.8p . which de- 
scribes spin filtering when there is no scattering out of the ring, except for the 
values of the matrix entries. Solving this system gives eigenvalues Ai = and 
A2 = — ( W- + + W + - ), leading to eigenvectors 



Vi 



and 



v 2 



f 



And the solution to the system is 

N+(t) 
N-(r) 



c lVl e AlT + c 2 v 2 e A2T 





" W-+ ' 




f 


Cl 




e° + c 2 










- I 



-(W-+ + W+-)t 

1 

where C\ and c 2 are constants to be determined from the initial conditions, hence 

N + (t) = Cl W^ + + c 2 e- {w - + + w+ - )T , 

iV_( r ) = Cl W + - - C2e -(w- + + w + .)r ? (5>L2) 

where one sees that N + (t) + AT_(r) = c\ ( + W+- ) = constant = N , and 



1 We denote the time variable in each of the dynamical systems by r to avoid confusion with the 
squared momentum transfer (Mandelstam t variable) used throughout the thesis. 
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2 4 6 8 10 12 

T 

Figure 5-1: A schematic graph of the Sokolov-Ternov effect, showing that the number of 
particles in the 'spin up' state (N + ) decreases with time while the number of particles in 
the 'spin down' state (N-) increases with time. Hence \P\ = \{N + — N_)/{N + +iV_)|, the 
absolute value of the beam polarization, increases with time. One sees that the relations 
&N+/&T = — d-/V_/dr and N + + iV_ = Nq are satisfied throughout. The graph just shows 
general trends, therefore we do not specify the units of the time axis. 



also that d N + / d r = — d N_ / d r as it must be if there is no scattering out of the 
ring. Imposing the initial conditions N + (0) = AL(0) = Nq / 2, corresponding to a 
beam that is initially unpolarized, gives the constants 

ci = — — — and c 2 = — 



w_ + + iy+_ 2 iy_ + + iy+_ ' 

and thus the complete solution 

v(t) - N +~ N - - w -+ ~ w +- n _ e -(w. ++ w + . ) r ] f513 x 

The spin-flip transition rates and W-+ are defined from the theory of syn- 
chrotron radiation as 



where EI 



5y/3mp 2 R 

W ° = —Z FT' 5 - L5 

8 r 
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all quantities denned as in eq. (12.3. ip which in fact, as we shall show, is the reciprocal 
of the above equation. Adding and subtracting the spin-flip transition rates gives 



W_ + -W + _ = f^=W . (5.1.6) 
5 v<J 

One can now present the complete solutions for the number of particles in both the 
'spin up' and 'spin down' states as a function of time r : 

N + (r) = f + ^ (.->*' - 1) , (5.1.7) 

JV-(r) = f - ^T 3 (e- W ° T " 1). (5.L8) 

which can trivially be seen to satisfy N(r) = N+(r) + N-(r) = Nq = constant. 

The steady state polarization (Sokolov-Ternov polarization) Vst is reached when 
r > 1 / ( W_ + + W+- ) = W^ 1 = r ST where r ST is presented in eq. (12.3. ljl R 

W , - W+ -8 

^ = wttwt- - ^ - ^ 

and one has the complete solution 

V(r) = V ST [1 - e- WoT ] = [l - e - T / rsT ] . (5.1.10) 

5 v3 

In a perfect ring one obtains an equilibrium polarization of 92.4% after time tst- 
In practice the maximum polarization achieved is slightly less than this due to im- 
perfections in the magnetic fields of real synchrotrons. Some parameters, including 
Tst, of current and proposed future synchrotrons are presented in Table 12.11 The 
general behaviour of the number of particles in the 'spin up' and 'spin down' states, 
along with the polarization buildup, due to the Sokolov-Ternov effect are plotted in 
Figure 15-11 



2 For times well above r = tst, (l — e~ r / TST ) w 1. The time tst depends strongly on the 
Lorentz 7 factor and the mean radius R of the storage ring but is typically of the order of minutes 
or hours for electron storage rings, see Table ETT1 
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5.2 Polarization evolution equations 



a out 




a out 

Figure 5-2: These diagrams describing the two physical processes, selective scattering out 
of the beam (left) and selective spin-flip while remaining in the beam (right), that contribute 
to polarization buildup by spin filtering in a storage ring have been explained in Figure \2-4\ 
Here we label the cross-section for a particle in the 'spin up ' state to be scattered out of the 
beam as <r^ ut , the cross-section for a particle in the 'spin down' state to be scattered out of 
the beam as <7° ut ; the cross-section for a particle in the 'spin up' state to be flipped to the 

'spin down' state while remaining in the beam as cr_| and the cross-section for a particle 

in the 'spin down' state to be flipped to the 'spin up ' state while remaining in the beam as 

a In order for each of these processes to contribute to polarization buildup of the beam 

we must have o~° ut ^ o~° ut and ^ a y respectively. These cross-sections will be used 

below in the mathematical evolution equations to describe the rate of buildup of polarization 
by spin filtering. 

In this section we develop sets of differential equations that describe the buildup 
of polarization of an antiproton beam by spin filtering. Consistency checks are then 
performed on the systems of equations, which provides a chance to highlight the 
underlying physical phenomena under investigation. The method of polarization 
buildup by spin filtering has been outlined in section I2.3.5[ which the reader may 
wish to read again before continuing here. Since the cross-sections for an interaction 
between a beam particle and a particle in the target are low we neglect the effects of 
multiple scattering, which has a cross-section orders of magnitude smaller than single 
scattering. Hence we consider two possibilities each time the beam passes through 
the target: (1) a beam particle can pass through the target without interaction, or 
(2) a beam particle can scatter off at most one of the target particles. 

Interactions with residual gas due to a non-perfect vacuum in the storage ring can 
also be neglected as the density of the internal target is much higher than the density 
of the residual gas. The high density of the target, and the fact that it is a constantly 
replenished gas jet, ensures that there is no significant target depolarization. 
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Recall the two physical processes that contribute to spin filtering shown in Fig- 



ure 15-21 The number of particles in the 'spin up' state can change by three means: 
(1) 'spin up' particles being scattered out of the beam, the cross-section for which 
we label as cx° ut , (2) 'spin up' particles being flipped to 'spin down' particles while 
remaining in the beam, the cross-section for which we label as and (3) 'spin 
down' particles being flipped to 'spin up' particles while remaining in the beam, the 
cross-section for which we label as cr_ + . Mechanisms (1) and (2) constitute a decrease 
in the number of 'spin up' particles and (3) constitutes an increase in the number 
of 'spin up' particles. This explains the signs of the coefficients in eq. (15.2. ip which 
follows. Correspondingly the number of particles in the 'spin down' state can also 
change by three means: (1) 'spin down' particles being scattered out of the beam, 
the cross-section for which we label as cr° ut , (2) 'spin down' particles being flipped 
to 'spin up' particles while remaining in the beam (er_ + ) and (3) 'spin up' particles 
being flipped to 'spin down' particles while remaining in the beam (cr + _). All of this 
can be expressed in the following set of polarization evolution equations: 
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' (7 + out + a+_ - a_+ 
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d7 




= — n v 








AL 




- a + _ a_ out + _ 




AL 



(5.2.1) 



where r is the time variable, n is the areal density of the target, v is the revolution 
frequency of the beam and N + (t) and AL(r) are the number of beam particles in 
the 'spin up' and 'spin down' states at time r respectively. 

For a beam that is initially unpolarized one imposes the following initial condi- 
tions 

N + (0) = AL(0) = ^, (5.2.2) 



where Nq is the total number of particles in the beam initially. We define the beam 
intensity N(t) = N + (r) + AL(r) as the total number of particles in the beam at 
time r, and the beam total spin at time r as J(r) = N + (r) — AL(r) so that the 
polarization of the beam at time r is simply given by 



P(r) 



N+(t) - AL(r) 



m_ 

N+{t) + N_(r) N(t) 



(5.2.3) 



The change in beam polarization as the number of particles in each of the spin states 
changes is very important throughout the thesis, so we plot this dependence and 
highlight a few points in Figure 15-31 

A problem with spin filtering where particles are scattered out of the beam is that 
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Polarization versus N+l AL ratio 
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Figure 5-3: This graph, from equation \5.2.3l shows how the polarization (P) changes as 
the ratio (N+ / N-) of the number of particles in the 'spin up' state to the number of 
particles in the 'spin down' state changes. The horizontal axis is a log scale. The change in 
polarization as the number of particles in each of the spin states changes is very important 
throughout the thesis, so we highlight a few points in the table. The second last entry in 
the table is particularly relevant as the target used in our numerical calculations is 90% 
polarized. Note that < \V\ < 1 and when iV_ > N + the polarization is defined to be 
negative, i.e. —1 < V < 0. 

while the beam polarization increases the beam intensity decreases. We propose 
possible solutions to this problem in Chapter O The behaviour of the number of 
particles in the 'spin up' and 'spin down' states, along with the polarization, as time 
increases is shown in Figure 15-41 

Note some treatments of spin filtering investigate a scenario where no particles 
are scattered out of the beam, i.e. the maximum scattering angle for the process is 
less than the ring acceptance angle, which is the case for antiprotons scattering off 
electrons in an atomic target [681 ESI EH E21 EE] and for antiprotons scattering off a 
co-moving beam of electrons or positrons |75j. In these scenarios only selective spin- 
flip can contribute to polarization buildup, and one avoids the problem of decreasing 
beam intensity. The low density of the targets currently available causes the rate of 
polarization buildup using these methods to be slow, but the enhanced cross-sections 
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Figure 5-4: A schematic graph showing that the number of particles in the 'spin up' (N + ) 
and 'spin down' (N-) states each decrease with time, but at different rates. Hence the beam 
polarization V = (N + — N-)/(N + + N-) increases with time. The time axis is scaled by 
the beam lifetime r*, as described in section \5.3.1\ 

at low energies suggested in Refs. [73 ESI 1107] may counteract this difficulty. We 
analyze such systems later in the thesis. 

Before solving this system of polarization evolution equations we shall prove a 
number of short lemma's providing a consistency check that the equations accurately 
describe the physical phenomena we wish to model. This also provides a chance to 
highlight the dynamical properties of the physical system, as this plays a major role 
in the rest of the thesis. 

Lemma 1 

If o+ ut = o"° ut and = there will be no buildup of beam polarization, but 

there will still be loss of beam intensity N(t). 

Proof: 

When <7° ut = <r° ut and o~ + _ = a h the polarization evolution equations reduce 

to 
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i.e. 



nv [(a + out +a + -) N + - <t+_ AL ] , 

nv [-a + ^N + + (a + out + <r + _ ) N_ ] , (5.2.5) 
which can be added and subtracted to give 

nua° ut N{T), 

n v ( (j + out + 2 a + _ ) J(r) , (5.2.6) 

two uncoupled first order separable ODE's which can be integrated to give the solu- 
tions 

N(t) = N(0) e- nua + utr = N e- nua + utr , 
J( T ) = J(0) e -^( CT + out + 2 -+-h, (5.2.7) 

and we see that N(t) will decrease exponentially, provided that cr + out ^ 0, and J(r) 
which is zero initially will always remain zero (i.e. if J(0) = then J(r) = for 
all r). So there will be no polarization buildup. Also in the case when the beam 
is initially polarized ( J(0) ^ 0) its polarization will decrease exponentially to zero, 
remembering the cross- sect ions are positive quantities. □ 

We later show that when the internal target is not polarized <r° ut = a° ut and 
c + _ = (J_ + . Thus there will be no polarization buildup by spin filtering if the in- 
ternal target is unpolarized. 

Lemma 2 

If c+ ut = °"° ut = there will be no loss of beam intensity N(t) = Constant = N , 
but there may still be polarization buildup. 

Proof: 

When a° ut = <x° ut = (this happens when the maximum scattering angle for the 
process is less than the ring acceptance angle) the polarization evolution equations 



dN + 
dAL 



djV(r) 
dr 

dJ(r) 
dr 
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Figure 5-5: A schematic graph of the system treated in Lemma 2 and Corollary 1, showing 
that the number of particles in the 'spin up' state (N + ) increases with time while the 
number of particles in the 'spin down' state (N-) decreases with time. Hence the beam 
polarization V = (N+ — N-)/(N+ + iV_) increases with time. One sees that the relations 
diV_|_ /dr = — diV_ / dr and N = N + + iV_ = Nq are satisfied throughout. The graph just 
shows general trends, therefore we do not define the units of the time axis. 

reduce to 
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(5.2.8) 



i.e. 

dN_ 

and adding these gives 
dN d 



-nu ( cr + _ N + - o-_ + N_ ) , 
-nu (-a + _N + + o-_ + AL) 



dr 



dr dr 



+ 



= —nu{ (T + _ Nj, 



dAL 
"dr"' 

o+-N + + cr_ + iV_ 



(5.2.9) 



0. 



So we have d/dr [iV + (r) + iV_(r)] = which implies N + (r) + iV_(r) = 
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constant = N(0). Thus there will be no loss of particles, as expected. Subtract- 
ing eqs. (I5.2.9P from each other gives 



dJ 

d7 



— ( N + - N. 
dr 

- 2 n v ( a. N, 



— a. 



dN + 
~dr~ 



diV_ 

~dr~ 



which leads to a non-zero J{r) {i.e. non-zero polarization) provided that er + _ 7^ 0"-+- 
□ 

The system without scattering out of the ring described in eq. (I5.2.8P is very 
similar to the system which describes the Sokolov-Ternov effect presented in sec- 
tion 15. 1[ The systems describing these two physical processes should be similar as 
they are both governed solely by spin-flip transitions. In spin filtering the spin- 
flip transitions are induced by scattering off the polarized internal target while in 
the Sokolov-Ternov effect the spin-flip transitions are induced by spontaneous syn- 
chrotron radiation emission of photons while the charged particles of the beam are 
being bent in the magnetic field of the ring. In fact these systems are identical except 
for the interpretations of the matrix entries, and that because there is no target in 
the Sokolov-Ternov effect, the system of equations describing it does not depend on a 
target areal density n. The solution of the system presented in eq. (I5.2.8P is identical 
to the solution of the Sokolov-Ternov system presented in detail in section 15 . 11 



Corollary 1 



When there is no scattering out of the beam, i.e. a° ut = <x° ut = ; the condition 

diV, dAL 



dr dr 

must be satisfied. 



(5.2.10) 



Proof: 



We have shown that when there is no scattering out of the ring the equations reduce 
to eqs. (I5.2.9P which can immediately be seen to satisfy eq. (15.2. 10p . □ 
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Lemma 3 



When there is no spin-flip, i.e. cr + _ = = 0, the change in one spin state 

should not depend on the number of particles in the other spin state, i.e. the equa- 
tions should decouple. 



Proof: 



When there is no spin-flip o~ + _ = o~_ + = and thus the polarization evolution 
equations reduce to 



diV + 

~d7 



nva° ut K 



nva° ut N_, (5.2.11) 



which is an uncoupled system of equations as required. □ 

The solutions of the above equations, for an initially unpolarized beam, are easily 
found to be 

N + (r) = ^ e -n- + ° ut - and N _ {t) = 1^ e -nua J^r _ (5^3) 

One sees that if in addition er° ut = a° nt in eqs. (15.2.121) then no polarization buildup 
occurs, as shown in Lemma 1. It is claimed by the Budker-Julich groups that the spin- 
flip transition rates are negligible for antiprotons scattering off polarized electrons in 
a hydrogen target [TT], [T2]. As we show in Chapter the maximum scattering angle 
of antiprotons scattering off atomic electrons is 0.54 mrad, below the acceptance 
angle of a typical storage ring, thus there is no scattering out of the beam. Since 
there is no scattering out of the beam, and spin-flip transitions are negligible the 
Budker-Julich groups conclude that polarized electrons in an atomic target are not 
effective in transferring polarization to an antiproton beam by spin filtering [71, 72J. 
To force some antiprotons to be scattered out of the beam, and to avoid the problem 
of loss of beam intensity due to antiprotons annihilating with the protons in an 
atomic target, we suggest to use an opposing polarized electron beam of sufficient 
energy to scatter some antiprotons beyond the ring acceptance angle [75]. Such a 
system is treated in Chapter [71 The rate of polarization buildup using this method 
is slow due to the low densities of polarized electron beams currently available [76J, 
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but the enhanced cross-sections at low energies suggested in Refs. [7J3 
counteract this difficulty. 

Differentiating eq. (I5.2.3P and rearranging gives 
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(5.2.13) 



and on substituting in eqs. (15.2.111) one obtains 
dV 



— nv 



dr 



out 
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out 
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(5.2.14) 



which is exactly the important equation (3) from Ref. [69], the first theoretical de- 
scription of the FILTEX results. Integrating the above equation, and using a result 
from the next section eq. (I5.2.2ip . leads to 



V(t) = tanh 



— nv 



a 



out 



a _out j 



tanh ( — nv Vt A out r 



(5.2.15) 



Which was proposed initially as a model of the rate of polarization buildup in spin 
filtering [54, 69J. Since then the importance of scattering within the beam has been 
highlighted and a more complex treatment of spin filtering is required, involving the 
polarization transfer to and depolarization of particles scattering within ring accep- 
tance. Therefore it is seen that when spin-flip effects are neglected the theoretical 
treatment of spin filtering presented in this thesis reduces to the initial naive treat- 
ments where polarization transfer and depolarization effects of scattering within ring 
acceptance were not included. 



5.2.1 cr + out , a _ out , 0+ , a _ + and the spin observables 

The spin observables of a spin 1/2 - spin 1/2 scattering process are defined in sec- 
tion 13.51 In spin filtering where the polarization of the recoiled target particle is not 
important one is interested in the polarization transfer, depolarization and double 
spin asymmetry spin observables. These have been calculated for electromagnetic 
antiproton-proton and antiproton-electron elastic scattering in Chapter HI The spin 
transfer observable has been calculated for low energy antiproton-positron scattering 
in Ref. [88]. A large increase of this spin transfer cross-section at very low energies 
is the basis for the proposal to polarize antiprotons by interaction with a co-moving 
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polarized positron beam presented in Ref. [75]. It is claimed in Ref. |107] that the 
polarization transfer cross-section for ep or e + p (like charges) scattering is enhanced 
at very low relative velocities, but by much less than that claimed in Refs. [72J 155] . 
An experiment has been proposed to test, and distinguish between, these claims 

msi. 

The cross-sections cr + out , o~_ out , o~ + _ and cr_ + can be related to the spin observables 
that have been calculated in Chapters ([3] and H]) by the following relations 



< ut = /out + V T A out , (5.2.16) 
a_ out = J out - V T A outl (5.2.17) 

<r+- = L in + ^ (A in - K in ) , (5.2.18) 

(7_+ = L in - ^ (A in - K in ) , (5.2.19) 

where Vt is the polarization of the target, and L in = (J in — D in ) / 2 is a loss of 
polarization quantity. These relations involve integration of the spin observables 
presented in Chapter 0] over the following angular ranges, as seen in Table 15.11 The 
"in" subscript refers to particles that are scattered at small angles < 9 acc remaining 
in the beam, and the "out" subscript refers to particles that are scattered out of 
the beam. Thus the integrals over scattering angle 9 are labeled "in" where the 
range of integration is 6 m - m < 9 < 9 acc , "out" where the range of integration ijf| 
#acc < 9 < 6> max , where # max is the maximum scattering angle for the process in the 
given reference framcl, and "all" = "in" + "out" where the range of integration 
is #min < 9 < #ma X ; as seen in Table I5TT1 I = da / dQ is the spin-averaged differential 
cross-section and A, K and D are the double spin asymmetry, polarization transfer 
and depolarization spin observables respectively as calculated in Chapter HI 



3 While not occurring in the case of antiprotons which we focus on here, an additional effect must 
be accounted for in the case of polarization buildup of a proton beam by spin filtering off a hydrogen 
target [120] . as in the FILTEX experiment. Because the final state particles are identical and hence 
indistinguishable, w-channelpp scattering can contribute; i.e. protons from the hydrogen target can 
be back scattered into the circulating proton beam. This happens when the beam protons are back 
scattered into the angular range (ir — # acc ) < 8 < it in the CM frame. This effect can be accounted 
for in the above formalism by changing the angular ranges in the CM frame to 9 mnl < 9 < # acc plus 
(it — # acc ) < < tt for the "in" integrations and ^acc < 9 < (tt — 6* acc ) for the "out" integrations. 
The physical result of this effect is to lessen the rate of decrease of beam intensity. 

4 In the CM frame = tt corresponding to total backward scattering, but in other frames, 

for example the LAB frame, this extreme value is not reached for some reactions and 9^ x < tt. 
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All cross-sections and spin observables contributing to spin filtering are azimuthally 
averaged, due to the geometry of the scattering, where the scattering plane can be 
at any azimuthal angle. Hence single spin observables, for example the analyzing 
power, do not contribute to the polarization evolution equations because they vanish 
when azimuthally averaged. The cylindrical symmetry of the system also implies 
that the ring acceptance angle has no azimuthal dependence. 
Note the following linear combinations of the cross-sections 

(5.2.20) 



V T A ont = -± g - ^' (5 ' 2 ' 21) 

L in = -^—^ ± , (5.2.22) 

V T (A in ~ K- m ) = a + _ - (5.2.23) 

Again to ensure consistency and to highlight the physical properties we are trying to 
describe mathematically we prove a number of short Lemma's on the above relations 
between the cross-sections and the spin observables. 
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Lemma 4 



// the target is unpolarized (Vt = 0) then one has that cr + out = cr_ out and cr + _ = a_ + 
so no polarization buildup will occur. 



Proof: 



Setting V T = into the eqs. (15.2.161 15.2. 171 15.2.181 and 15.2.190 one immediately 
obtains cr + out = J out = a_ out and a + _ = L in = Once this is satisfied it is 

proved in Lemma 1 that no polarization buildup will occur in this case. □ 



Lemma 5 



The spin-flip cross-sections should depend only on spin observables relating to parti- 
cles scattering within the ring, i.e. only to "in" spin observables which are integrated 
from 9 min to 9 acc . 
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Proof: 



This is immediately satisfied by the relations in eqs. (15.2. 181 and 15.2. 19]) . □ 



Lemma 6 



The cross-section differences a + out 
to the target polarization Vt- 



a 



out 



and 



o~_ + should both be proportional 



Proof: 

This is immediately satisfied by the relations in eqs. (15.2.211 and 15.2.25]) . □ 

While the system of polarization evolution equations involving the variables 
N + (t) and AL(r) presented in eq. (15.2. ip is very transparent, one is more inter- 
ested in the variables N(t) = N+(t) + N_(r) and J(r) = N + (r) - iV_(r) which 
immediately lead to V(r) = J(r) / N(t). We can transform the system of two first 
order ODE's in variables N + (t) and N-(r) presented in eq. (15.2. ip to the following 
system of two first order ODE's in variables N(t) and J(r) [721 H21j : 
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(5.2.24) 



where we have also transformed from the cross-sections to the spin observables which 
have already been calculated. The parameters n and v are the target areal density 
and the beam revolution frequency respectively. 

The systems presented in eqs. (15.2.11 and I5.2.24p are identical and from now on 
we concentrate on the latter as its solution is more illustrative of the underlying 
physical phenomena, and the dependence on the target polarization is explicit. In 
particular one immediately sees that when the target is unpolarized no beam po- 
larization buildup occurs, as when Vt = the system reduces to two uncoupled 
separable first order ODE's as in eq. (I5.2.6P with solutions as presented in eq. (15.2.71) 
showing V{t) = for all r if Vt = 0. The parameters in the matrix of coefficients 
of eq. (I5.2.24p depend on the state of the target polarization, i.e. longitudinal or 
transverse, as seen in Table 15.11 
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Transverse polarization requires 



Longitudinal polarization requires 



' mB f da 
/out = 27r/ { — ) S m6d6 



^4 out — 2 7T 



dn 

Axx + A Y y d a 



K m = 2 ir 



D in = 2tt 



— I sin0d0 
d 12 



A xx + A Y y da\ . 

— smVdv 

2 dQ 



^min 



-Kxx + K YY da 

2 dVl 

Dxx + Dyy da 

2 d£l 



sin d# 



sinfldtf 



out 



2tt 



out 



2tt 



A =11 = 2 7T 



K iri = 2vr 



D«„ = 2tt 



4^ 1 sin 0<10 
dil 



A zz ^ ] sintfdtf 



Kzz ^ ) sin0d0 



D zz ^ Ism.d. 



Table 5.1: The entries in the system of equations for polarization buildup involve angular 
integration over the spin observables presented in Chapter^ X, Y and Z are the coordinate 
axes where the beam is moving in the positive Z direction. The minimum value for 9 (0 m i n ) 
relates to the average transverse electron separation for a pure electron target and to the 
Bohr radius for an atomic gas target, 6> acc is the ring acceptance angle and # max is the 
maximum scattering angle for the process. 

5.3 Solving the polarization evolution equations 

We have shown in the previous section that when circulating at frequency u, for a 
time r, in a ring with a polarized internal target of areal density n and polarization 
Vt oriented normal to the ring plane, (or longitudinally with rotators) 
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(5.3.1) 



describes the rate of change of the number of beam particles N(t) = N+(t) + AL(t) 
and their total spin J(r) = N + (t) — iV_(r) [72]. In this section we solve this system 
of polarization evolution equations. The eigenvalues of the matrix of coefficients are 
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found to be 



Ai = -nv (/ out + L in + L d ) and A 2 = -nv (J out + L in - L d ) , (5.3.2) 
where the discriminant L d of the quadratic equation for the eigenvalues is 



L d = ^/P^A out (A all - K in ) + Ll . (5.3.3) 

Note that J ut, L in and L d are all non- negative. As a consequence the eigenvalues are 
non-positive and Ai < A 2 < 0. When there is no scattering out of the ring all of the 
"out" integrations are zero, and one finds that L d = L- m and hence Ai = —2nvL d 
and A 2 = 0. 

Now enforcing the initial conditions N(0) = iVo the total number of particles 
in the beam initially, and J(0) = iV+(0) = AL(0) = N /2 i.e. initially the 
beam is unpolarized, one obtains the solutions: 

N(t) = — J - , (5.3.4) 

ZL d 

2L d 

The time (r) dependence of the polarization of the beam is given by 

P(r) = JM= -^U all -K m ) 
w JV(r) L in + L d coth(L d nz/r) ' v--; 

The expression for P{r) is proportional to Pt which confirms that if the target 
polarization is zero there will be no polarization buildup in the beam. The approxi- 
mate rate of change of polarization for sufficiently short times, and the limit of the 
polarization for large times are respectively: 

^ « -nvP T (A ail - K in ) , (5.3.7) 
dr 

P max = limF(r) = - P T A ; 11 " K in . (5.3.8) 

In order to compare to earlier treatments of spin filtering notice that in the absence 
of scattering within the ring, when all "in" spin observables are zero, one has that 
L d — Pt ^4 out and eq. (15.3.61) reduces to eq. (15.2. 1 5j) which was the initial treatment 
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of spin filtering proposed in 1993 where scattering within ring acceptance was not 
included. Having said that, the general behaviour of equations 15.2.151 and I5.3/61 are 
similar. 

For pure electromagnetic scattering the double spin asymmetries equal the po- 
larization transfer spin observables [90], thus one can simplify the above equations 
using A in = K in , A out = K out and A all = K^; hence A all - K m = K out . 

At this point one may wish to find expressions for N + (t) and iV_(r) which can 
easily be obtained from eqs. (I5.2.3[ 15.3.41 and I5.3.5[) : 



NJt) 



N(t) + J(r) 



^- {e AlT [L d - L in + (A aU - K in ) V T 



4L d 

+ e A2T [L d + L in - (A M - K in ) V T ] } , (5.3.9) 
N(r) - J{t) 



^- {e AlT [L d - L- m - (A all - K in ) V 7 



4L d 

+ e A2T [L d + L in + (Art - K- m ) V T \ } , (5.3.10) 

which can easily be seen to satisfy N + (0) = A^_(0) = N / 2, the correct initial 
conditions. If the target is unpolarized (Vt = 0) one sees that N + (t) = N-(t) 
therefore no polarization buildup will occur, as was required by Lemma 4. 

5.3.1 Beam lifetime and figure of merit 

The beam lifetime r#, the time taken for the beam intensity to decrease by a factor 
of e ~ 2.718, i.e. N (t* ) = N /e, can be obtained from eq. (15.3.41) . One finds 

r* w — , (5.3.11) 

n v I out 

a more accurate form of which is derived in section 15.3.31 

The Figure Of Merit (FOM) provides a measure of the quality of the polarized 
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Figure 5-6: A schematic graph showing the behaviour of the beam intensity N and beam 
polarization V as time (scaled by the beam lifetime t*) increases. The behaviour of the 
figure of merit FOM is also shown on the graph, being blown up to clearly show it has a 
maximum at twice the beam lifetime. 



beam, and is given by 

FOM(r) = V 2 (t) N(t 
The figure of merit for the above case is 

FOM(r) 



(A an -K in ) 2 N V 2 



2L d 



J\r) 
N(t) ■ 



(e AlT -e A2T ) 2 



(5.3.12) 



. (5.3.13) 



Maximizing the figure of merit gives the optimum polarization buildup time, taking 
into account the trade-off between decreasing beam intensity and increasing beam 
polarization. Solving d FOM / d r = yields 



t, 



optimum 



nvl 



2r* 



(5.3.14) 



out 



approximately twice the beam lifetime. Thus the optimum time for polarization 
buildup is twice the lifetime of the beam, as also found in Ref. [66] . 

The behaviours of the beam intensity, beam polarization and the figure of merit 
as time increases are shown in Figure 15-61 Note the characteristic trade-off of spin 
filtering: as the beam polarization increases the beam intensity greatly decreases. 
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In Chapter [6] we address possible ways of circumventing this drawback, for instance 
by continuously inputting particles into the beam. 



5.3.2 Pure electromagnetic scattering 



For pure electromagnetic scattering Aij = Kij where i,j € { X, Y, Z }, so the system 
simplifies to 



_d_ 

d7 



" N ' 






= —nv 


_ J _ 







' N ' 




_ J _ 



1 out 

V T K out 

^T-K'out I all — D{ n 

the coefficient matrix of which is symmetric, and the solutions become 

J(t) -K out V T 



(5.3.15) 



V{t) 



N(t) L in + L d coth [L A nvr) 
L d = JvSK 2 out + L\ 



(5.3.16) 
(5.3.17) 



The approximate rate of change of polarization for sufficiently short times, and the 
limit of the polarization for large times simplify to respectively: 



dV 
dr 



n v Vt K c 



V 



lim V{t) 



K out V T 



Lin + Ld 



(5.3.18) 



We now show that the maximum polarization achieved cannot exceed one, i.e. 
I^maxl < 1- From the definition of L d for the case of pure electromagnetic scat- 
tering we have L\ = V? K 2 ont + Lf n =>■ L d > [PrKcmtl which can be used to 
obtain 



\v I < hA < h± = i 

I' max l — T , T — J L 



(5.3.19) 



since L m is a non- negative quantity. □ 

Note the upper bound l^maxl = 1 only happens for L in = 0, i.e. no polarization 
is lost. 
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5.3.3 Fraction of antiprotons lost per revolution 

We now calculate the fraction of the antiproton beam that is lost per revolution. 
Define Ar to be the time taken for one revolution, thus u At = 1, and AiV to be 
the change in beam intensity during this time. Manipulating the d N / d r equation 
gives us 

AiV 

— — = -nvI out N - nvV T A out J , 
Ar 

__ = -nuI out Ar - nuV T A out — Ar , 
= ~nI out - nV T A out — , 

= -n(I out + V p V T A out ) . (5.3.20) 

Since AiV is the change in beam intensity during time Ar, and the beam intensity 
decreases with time, the quantity AiV is negative. But the fraction of the antiprotons 
lost per revolution is n (I out + VpVr A out ) which is positive, and thus the fraction 
of antiprotons remaining in the beam per revolution is 1 — n (/out + VpVr A out ). 
The fraction of antiprotons lost per revolution is not constant in time, it depends on 
the polarization of the antiproton beam Vp , which is zero initially. The fraction of 
antiprotons lost per revolution decreases as the polarization of the antiproton beam 
increases. This makes physical sense, since as the antiproton beam polarization 
increases there will be fewer particles in the spin state that is scattered out more 
often. 

This can also be used to obtain the beam lifetime. Since the frequency v is the 
number of revolutions per second, the fraction of particles lost from the beam per 
second is n v ( 7 out + Vp Vt A out ), hence the fraction of particles lost from the beam 
in r seconds is nvr (J ut + VpVx A out ) . The beam lifetime is the time taken for 
the beam intensity to decrease by a factor of 1 / e, i. e. the time taken for a fraction 
1 — 1 / e = (e — 1) / e of the beam particles to be lost due to scattering out of the 
ring. Thus we can calculate the beam lifetime (r*) by solving 

nvn(I out + V p V T A out ) = e —^, (5.3.21) 
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leading to the beam lifetime 



e- 1 



n v e ( J out + V p V T A out ) ' 



1 



1 



(5.3.22) 



nv (J out + V P V T A 

out J 



nue ( J out + V P V T A 

out J 



the first term of which is e « 2.718 time larger than the second. Note that since 
lout > 'Pp'Pr A out , this expression for the beam lifetime limits to the one calculated 
by the other method above are equal to leading order. Of importance here is the 
fact that the beam lifetime is not constant, it increases as the beam polarization 
increases; which makes physical sense as there will be fewer particles in the spin 
state that is scattered out more often. 
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Chapter 6 

Various scenarios of spin filtering 



"Let no one ignorant of Mathematics enter here. " 
Plato, inscription over the entrance to the Academy. 

A major problem with spin filtering is that as the beam polarization increases the 
beam intensity decreases, since particles are being continuously scattered out of the 
beam. While one may obtain a polarized antiproton beam its intensity may be too 
low to be of use in any experiment. In this chapter we investigate the possibility of 
continuously inputting unpolarized particles into the beam to counteract this loss of 
beam intensity. We present a thorough investigation of spin filtering under various 
alternate scenarios of interest to any practical project to produce a high intensity 
polarized antiproton beam. These scenarios are: 1) spin filtering while the beam 
is being accumulated, i.e. unpolarized particles are continuously being fed into the 
beam at a constant rate, 2) unpolarized particles are continuously being fed into 
the beam at a linearly increasing rate, i.e. the particle input rate is ramped up, 3) 
the particle input rate is equal to the rate at which particles are being lost due to 
scattering beyond the ring acceptance angle, the beam intensity remaining constant, 

4) increasing the initial polarization of a stored beam by spin filtering, and finally 

5) the input of particles into the beam is stopped after a certain amount of time, 
but spin filtering continues. The five sections of this chapter each treat one of the 
scenarios of spin filtering listed above, in that order. 
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6.1 Accumulation of antiprotons in the ring 



In the discussion so far we have only considered polarizing an antiproton beam when 
the beam is already accumulated in the storage ring. The VAX Collaboration plans 
to obtain their antiproton beam by collecting the produced antiprotons from high 
energy interactions of protons on targets of light nuclei, such as Beryllium. The 
antiprotons will be continuously fed into the storage ring at a fixed rate and accu- 
mulated, hence increasing the beam intensity, allowing for a greater luminosity in an 
experiment. The VAX Collaboration estimates the production rate of antiprotons 
as being 10 7 per second [23]. Since 10 11 antiprotons are required in the storage ring, 
antiprotons will be fed into the storage ring at a rate of 10 7 per second for 10 4 seconds 



We now consider a system where spin filtering occurs as the antiprotons are being 
fed into the ring. The original system of equations must be amended to account for 
this constant accumulation. The effect will be to add a term to the diV(r) /dr 
equation, where is the constant rate at which antiprotons are fed into the ring; while 
the d J(t) / d r equation remains unchanged. The initial conditions are N(0) = No , 
which will be set to zero in section IB. 1.1 1 and J(0) = 0. The new system of differential 
equations is 



dJV(r) 
dr 

dJ(r) 



-nv [J out N(t) + V T A out J(r)] + 0, (6.1.1) 
-nv [V T (A all -K in ) N(t) + (Jaii-^in) J(t)]. (6.1.2) 



By differentiating eq. (16. 1.2p with respect to r and substituting in eq. (I6.1.ip one 
obtains an inhomogeneous second order linear differential equation with constant 
coefficients for J(r): 

- (A! + A 2 )^^ + X 1 X 2 J(t) = - n vV T (A all -K in ) 0, (6.1.3) 
the solution of which is 

J{r) = ^Fr^ [A 2 (A 1 iV + /3)e Al - - X l ( A 2 N + ) e x ^ 

+ (Ai - A 2 )] . (6.1.4) 

Differentiating eq. f)6.1.4p with respect to r and substituting into eq. (I6.1.2p gives an 



87 



expression for N(r) : 
1 



N(r) 



[A 2 (XtNo + P) (L d -L in ) e Al ' 



2 L d Ai A 2 

+ Ai (A 2 iV + /3) (L in + L d ) e A2T + (3 (I an -D; 



(6.1.5) 
A 2 -Ai)] . 



As a consistency check one can see that these solutions for J(r) and N(r) satisfy the 
initial conditions J(0) = and N(0) = N , and in the particular case when (3 = 
the above expressions reduce to the solution of the homogeneous system presented 



in section 15.31 Dividing J(r) by N(r) we obtain an expression for the polarization 
V(r) as a function of time, 



V(r) 



L in + L . 



1 



A 2 e Mr( AlA r 0+/ 3)-/3 



Ai [e A 2r (A2iVo+( g)_ ( g 



(6.1.6) 



When the particle input rate is zero (i.e. = 0) the above equation simplifies to 

- (Aaii-liffc) P T - Mall - K in ) V T 



V(t) 



L d 



1 



>(Ai-A 2 )t 



1 



L in + L d coth(L d nz/r) 



(6.1.7) 



which is the solution of the homogeneous case presented in eq. f)5.3.6p . 

Using a Taylor Series expansion we find that the approximate initial rate of 
polarization buildup for each of these cases (N ^ with f3 ^ and N = with 
j3 7^ 0) is the same as in the homogeneous case (N ^ with f3 = 0): 



dV 
dr 



-nvV T (A an - Ki 



(6.1.8) 



The maximum polarization achievable is the limit as time approaches infinity 

-V T (A all -K- m ) 



lim V(t) 



all — -Din 



(6.1.9) 



which is independent of both N$ and (3, however note that in taking this limit we 
used the fact that f3 ^ 0. If (3 was equal to zero then the maximum polarization 
achievable would equal that from the homogeneous case; as can be easily seen from 
eq. (I6.1.6P remembering that Ai < A 2 < 0. Thus for the complete case there are just 
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two values of the maximum polarization, one for (3 = and one for all /3 ^ 0. The 
figure of merit for this inhomogeneous case is: 



FOM(r) = V 2 {t)N{t) = 



J_Ht) 
N(t) 



A 



all 



2 L d Ai A2 



x 



[ Cie Air _ C2 gA 2 r + /3 ( Al - A 2 ) ] 



ci (Ld-Lin) e A i- + c 2 (L in + L d ) e A ^ + /3 (I^-D^) (A 2 -Ai) 



(6.1.10) 



where for convenience we have defined the two constants Ci = A 2 ( Ai N + j3 ) and 
c 2 = Ai ( A 2 Nq + /3 ). Note the FOM will not have a maximum in finite time if the 
accumulation rate /3 is high enough to make the beam intensity a constant or increase 
with time. If this happens the FOM will increase monotonically. 



6.1.1 No initial beam 

Of interest is the particular case when N = 0, i.e. there are no particles in the beam 
initially. In this case the above solutions simplify to 



J(r) = 



(3V T {A^-Kj 
2 Ld Ai A 2 



[M 



1) + Ai (1 - e x * T )} , (6.1.11) 



N(r) = 



6 



2 L d Ai A 2 



[A 2 (L d -L in ) e AlT + Ai (Li n + Ld) e 



A2 r 



+ (hn-D in ) (A 2 -Ai)] , 



(6.1.12) 



7>(r) 



- (A all -X in ) P T 



in + < 



(l-e A i T ) A 2 
( l-e A 2^ ) Ai 



(6.1.13) 



FOM(t) 



A an -K hl ) 2 VS, 
2 Ld Ai A 2 



x 



[A 2 (e AlT - 1) + Ai (1 - e A2r )] 



(6.1.14) 



A 2 (L d -L in ) e A ^ + Ai (L in + L d ) e Aar + ( / a ii - D m ) (A 2 - Ai) ' 

Interestingly the (3 dependence of V(t) vanishes in this case, i.e. the polarization 
buildup rate is independent of the rate at which antiprotons are fed into the ring, if 
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there are no particles in the beam initially. But we have used the fact that (3 ^ 
to obtain the above result. We should note the obvious physical fact that if N = 
and P — i.e. there are no particles in the beam initially and no particles are fed 
into the beam, then there will never be any particles in the beam; so measuring the 
beam polarization is meaningless. Notice that in this case where N = 0, the figure 
of merit is proportional to the particle input rate (3. 

Summary 

The results obtained thus far for the polarization buildup in various scenarios can 
be summarized as 



V(t) 







- (A all - K iQ ) V T 
L in + Ld coth (L d n z/ r) 

- (A aU -K in ) V T 



1 - 



l-e A i ' 



I A 2 



- 1 



(A all -K in ) V T 



L in + L d 




2 


1 


1 ^ 2 


"e A i T (A 1 Af +/3)-/3 








e x 2 T (\ 2 N +l3)-l3 





for/3 = & A^ = 



for (3 = & N ^ 



for/3^0 & A^ = 



for/3^0 & iVo^O 



where, as usual, /3 is the constant rate at which particles are fed into the beam and 
Nq is the number of particles in the beam initially. 



6.2 Constant beam intensity 

In this case the accumulation rate is set specifically so that extra particles are fed 
into the beam at such a rate so that the beam intensity is kept constant, i.e. fed in 
at such a rate to balance the rate at which particles are scattered out of the beam. 
The system of equations is much simpler in this case. Here N(t) = N is a constant, 
hence dN(r)/dr = 0, and the J(r) equation becomes a first order linear ODE with 
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constant coefficients 

+ nv (/ a u -D^) J(r) = -nv {A all -K in ) V T N , (6.2.1) 

Cl T 

and imposing the initial conditions N(0) = N and J(0) = one obtains the solution 
J(r) = ~ U ;'r^n )PriV0 [l- e -"^- g ^ T ] ■ (6.2.2) 

( -t all — -Din ) 

Now the polarization as a function of time can be presented 
V(t) = I^_ J M= -(A,n-K m )V T _ {lM _ Din)T 

1 ' N(t) N (/.dl-Din) L J 

To find the maximum polarization achievable, i.e. the limit as time tends to infinity, 
we note that J an > D in thus —nu{ J a n — D in ) < and hence one obtains 

V m = hm P(r) = ~ U ; 11 '^" l)PT , (6.2.4) 

oo i a n — iri 



which is the same as in the inhomogeneous case of section IBTTl when (3 ^ 0. 

The initial rate of polarization buildup can be obtained by expanding V[j) as 
a Taylor expansion in nvr. Assuming nvr is small we neglect terms of second or 
higher order giving 

*L -nv (A all -K in ) V T , (6.2.5) 
a r 

as it was in the homogeneous case presented in eq. ( 15.3.71) . The figure of merit in 
this case is easily obtained 

FOM(r) = Uall-^in) 2 ^^ ^ _ e -n,(/ all -D in ),^ ( g 3 g) 
(-fall - D m ) 

and increases monotonically with time. 

The behaviour of the beam intensity, beam polarization and figure of merit as 



functions of time are plotted in Figure 16-11 One notes that in this case the figure of 
merit increases monotonically and hence does not have a maximum in finite time. 
Therefore in this scenario it would be best to continue the spin filtering process for 
as long as possible before extracting the polarized antiproton beam. 
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Figure 6-1: A schematic graph showing the behaviour of the beam intensity N , which 
is constant in this case, and beam polarization V as time (scaled by the beam lifetime 
T*) increases. The figure of merit is also shown on the graph, with the same scales as 
the other functions. The graph just shows general trends and is not in exact numerical 
correspondence to the equations, which will be presented in Chapter^ The beam lifetime 
is defined as earlier, in the absence of particle input to the beam. 

6.2.1 Approximating the critical input rate 

The accumulation rate needed to keep the beam intensity constant is important, 
as this critical rate divides the solution of the system into two physically distinct 
cases. Smaller accumulation rates than this critical value cause the beam intensity 
to decrease, hence the FOM will have a maximum in finite time. Larger values than 
the critical value cause the beam intensity to increase continuously, hence the FOM 
will increase monotonically. We can see from eq. fl5.3.4p that N(t) does not decrease 
linearly with time r. So the accumulation rate needed to keep the beam intensity 
constant, say /(r), will not be linear in r. We now derive the function /(r) and 
obtain a linear approximation to it, which can be used in the inhomogeneous case 
treated in section I6TT1 We must solve N(r) = N hom {r) + /(r) = N , i.e. 

\e XlT (L d -L iQ ) + e X2T (L d + L in )] N /s 
N{t) = I K — ^— K — + /(r) = N , (6.2.7) 

ZL d 

which leads to 

f(r) = [2L d -e AlT (L d -L in )-e A2T (L d + L in )] . (6.2.8) 

A La 
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A linear approximation /l(t) to /(r) can be found by Taylor expanding the expo- 
nentials to order r (which is valid since nvr is small), to obtain 



f L (r) = nvN I out r, (6.2.9) 

which is in the linear form + j3 c r where (3 C = nv Nq J out is the critical value of (3 which 
when added to the diV(r) /dr differential equation in section 16.11 approximately 
makes the beam intensity constant. 

6.3 The input rate is ramped up 

In this section we investigate a scenario where unpolarized particles are input into 
the beam at a linearly increasing rate, i.e. the input rate is ramped up. This is 
accounted for by the following system of polarization evolution equations 

<lA? ' i = -nv [J out N(t) + V T A out J{r)\ + 0r, (6.3.1) 



dr 

dJ(r) 
dr 



-nv [V T (Aau- K iQ ) N(t) + {I a n-D in ) J(r) ] , (6.3.2) 



where /3 r is the rate at which particles are fed in, the input ramped up at a rate 
proportional to the time elapsed. The initial conditions are N(0) = N which we 
may later set to zero, and J(0) = 0. By differentiating eq. (I6.3.2P with respect to 
r and substituting in eq. (16.3. ip one obtains a second order linear inhomogeneous 
differential equation for J(r): 

- (A 1 + A 2 )^^ + Ai A 2 J(t) = -nvV T (A an -K in ) 0r (6.3.3) 
cl r z dr 

the solution of which is 

J if) = F X2Xl e XlT + F XlX2 e X2T + (3 (A lT + A 2 ) . (6.3.4) 
Where for convenience we have defined the constants 

A 1 = — , (6.3.5) 

Ax A2 

. 2n 2 v 2 V T (A aU -K hl )(L m + I out ) 

2 = x[\l ' ^ ' 
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F \ 2 \i = 



nu ( A a n — K m ) N V T + (3 (A 1 -X 2 A 2 



(6.3.7) 



A2 — Ai 

obtained by imposing the initial conditions J(0) = and N(0) = N thus d J(0) /dr = 
— nv Vt ( A a u — if i n ) Ao- The function F\ 1 \ 2 is F\ 2 \ 1 with Ai and A2 interchanged. 
Differentiating eq. (I6.3.4p with respect to r and substituting into eq. (I6.3.2p gives an 
expression for N(r): 



N(t) 



-1 



\A<a- Kin) V T 



|F A2Al e Air (L in - L d ) + F XlX2 e X2T (L in + L ( 



^4 

— + (J out + 2L iD ) (A 1 r + A 2 
n v 



(6.3.8) 



As a consistency check it can be seen that the inhomogeneous solutions for J(r) and 
AT(t) satisfy the initial conditions, and that when = they reduce to the solutions 
of the homogeneous system eq. (I5.2.24p presented in section 15.31 

Dividing J(r) by N(t) we obtain an expression for the polarization as a function 
of time (r), 



V(r) 



-V T (A m -K 



L in + L d 



T F X2Xl (A 2 -Ai)-/3[Ai (1-A 2 t)-A 2 A 2 ] 



1 - 



e A 2 



' A 2 



(A 1 -A 2 )-/3[A 1 (1-A 1 r)-A 1 A 2 ] 



(6.3.9) 



When /3 = the above equation simplifies to 

-V T (^11 - K h 



V(r) 



L in + L d coth (L d n u r) 



(6.3.10) 



which is the solution of the homogeneous case eq. (15.2.241) presented in section 15.31 
Of interest is the case when N(0) = Nq = 0, i.e. there are no particles in the 
beam initially. To obtain this result we set N Q = in the above equation to obtain 



-V T (A^-Ki 



L in + L ( 



e A i T -l) A 2 A2-A1 (e A i T +A 2 r-1 
'e A 2 T -l) Ai A 2 -Ai (e A 2 T +Ai r-1 ' 



- 1 



(6.3.11) 



where for (3 7^ the /3 dependence vanishes. Using a Taylor Series expansion we 
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obtain the approximate initial rate of polarization buildup 

dV 



dr 



nvV T (A m - K ia ) , (6.3.12) 



identical to that of the homogeneous case presented in section 15.31 The maximum 
polarization achievable is the limit as time approaches infinity: 

P max = km V{t) = = — — . (6.3.13) 

r^oo 1 all — U in 1 out + I L in 

The above expression is only valid for ^ 0, the /3 = expression is presented in 



section 15.31 

For this inhomogeneous case the figure of merit is: 
FOM(r) = V\r)N(r) = £M = (6.3.14) 

-V T (A aU -K in ) [F X2Xl e x ^ + F XlX2 e x ^ + g (At + A 2 )} 2 
F X2Xl e^(L in + L d )-F XlX2 e^(L ia + L d )+/3[AL + (J all - D- m ) (A l t + A 2 )]' 

If the particle accumulation rate t is high enough to make the beam intensity 
constant or increase with time the figure of merit will be a monotonically increasing 
function of time, i.e. it will not have a maximum in finite time. 



6.4 Stored beam with initial polarization 

We now solve the homogeneous system where the initial polarization is not zero. This 
will be used if two methods of polarizing antiprotons are combined, i.e. if antiprotons 
were produced with a small polarization by some other method and one wanted to 
increase that polarization by spin filtering in a storage ring, where the luminosity 
could also be increased. In this section the beam has been stored and there is no 
further input of particles into the beam. 

The system of differential equations, the eigenvalues and eigenvectors are the 
same as section 15.31 but one of the initial conditions is different. The new initial 
conditions are N(0) = No > the total number of particles in the beam initially, and 
J(0) = Jo 7^ =>- iV + (0) 7^ AL(O) i.e. initially the beam is polarized. Note that since 
the number of particles in the beam in one particular spin state must not be greater 
than the total number of particles in the beam the bound |Jo| < iVo is respected. A 
negative value for J simply implies that the antiproton beam is initially polarized 
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in the opposite direction to the polarization direction of the target. Enforcing these 
initial conditions leads to the solutions 

( J V T A out - N L in ) (e x ^ - e x ^) + N L d (e x ^ + e x ^) 
N {V = 7TJ " , 

ZL d 

(6.4.1) 

[N V T (A^-K in ) + J L in ] ( e A ^-e A ^) + J L d ( e x ^ + e x ^ ) 

J(T) = 2L d ' 

(6.4.2) 

which reduce to the solutions of the original homogeneous system presented in 
eqs. (15.3.41 and I5.3.5j) when J ~~ 0. The beam lifetime is the same to leading 
approximation as in the homogeneous case when J(0) = 0. Dividing J(r) by N(t) 
provides the time dependence of the polarization of the beam 

v( v L d J - tanh (L d nvr) [L in J + N V T (A all - K hl )] 
L d N + tanh(L d nz/r) [L in N - J V T A out ] 

Denoting the initial polarization V(0) = J / N = Vo the above can be written as 

v( s _ L d V - tanh(L d nz/r) [L in V + V T ( A aU - K in )] 
[T) ~ L d + tanh(L d nz/r) [L in - V V T A out ] [ ' ' j 

The approximate rate of change of polarization for sufficiently short times is found 
by Taylor expanding to first order in r 

^ « nu { [A ont V 2 - (A aU -K in )} V T - 2V L in } . (6.4.5) 
The limit as time goes to infinity of V(t) in eq. (I6.4.4P is 

hm V(r) = ^^-L m); V T (A all -K m ) 

which of course agrees with the earlier maximum polarization if J = (i.e. Vo = 0). 
The figure of merit for this case is: 

FOM(r) = V 2 (t)N(t) = ^ = (6.4.7) 

N { [V T (A all -K in )+V L in ] (e x ^ -e x ^)+V L d (e x - + e x ^) } 2 
2L d {(V V T A out - L- m ) (e x ^ - e A ^) + L d (e x ^ + e A ^)} 
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Note that in terms of the number of particles in each spin state the initial conditions 
for an initially polarized beam, N(0) = No and J(0) = Jo, are 

N+(0) = ^ (1 + Vo) and AL(0) = ^ (1 - V ) . (6.4.8) 
6.4.1 An unpolarized target 

A special case of this system deserves comment. Given that the beam is initially 
polarized what happens if the target is unpolarized? One would imagine that the 
beam polarization should decrease and eventually reach zero. We now analyze the 
equations of section EH when the target is unpolarized (i.e. Vt = 0) and use the fact 
that Ld = a/ Vt A out (A a \\ — K m ) + L 2 n = L in when Vt = to obtain the beam 
polarization function of time 

V(t) = P e {Al ~ Az) r , (6.4.9) 

which is an exponentially decreasing function of r for Ai — A2 < 0. The beam 
polarization will not decrease in the special case of Ai — A 2 = 0, but this only 
happens when L in = 0, i.e. when there is no depolarization. The special case of 
Ai — A2 = 0, which does not lead to polarization buildup as seen from eq. ( 16.4.2]) . 
would be avoided by any experimental effort, thus is omitted from the rest of the 
discussion. 

The limit of beam polarization for large times when Vt = is 

lim V(t) = 0. (6.4.10) 

T— > OO 

Thus, as expected, if the beam is initially polarized and the target unpolarized then 
the beam polarization will decrease with time and eventually the beam polarization 
will reduce to zero. Thus a beam cannot gain polarization from an unpolarized target 
by spin filtering. 

The figure of merit in this case simplifies to FOM(r) = NqV 2 e ( - 2Xl ~ X2 ^ T which 
is a monotonically decreasing function of r. One can derive a polarization half-life in 
this case, the time taken for the polarization to decrease by a factor of 2, by solving 
V( T )=V /2 to obtain 

ln2 In 2 

Ti = = 

2 A 2 — Ai 2nvL A 



(6.4.11) 
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This scenario occurs in an electron cooler, a device used to focus the beam in many 



storage rings j. The beam passes through a co-moving beam of unpolarized electrons 
with low transverse momentum, in order to dampen the transverse momentum of the 
antiprotons in the stored beam. But the low electron areal densities in cooler beams, 
where typically n w 10~ 18 fm~ 2 = 10 -19 mb _1 , causes the polarization half-life 
to be very large. Thus our work shows that this depolarization effect is negligible, 
in agreement with Ref. |1U6] . Since electron cooling is a necessary part of the spin 
filtering process of polarization buildup it is very important that the depolarization 
caused by electron cooling is negligible. 

It has recently been suggested that the positron-antiproton polarization transfer 
observable is very much enhanced at low energies [881 1 107] . This enhancement is the 
basis of the recent proposal by Th. Walcher et al. to polarize an antiproton beam by 
repeated interaction with a co-moving polarized positron beam in a storage ring [75 J . 
All of the antiprotons remain within the beam in this scenario and one avoids the 
problem of the antiprotons annihilating with protons in an atomic gas target. This 
large enhancement of the polarization transfer observable at low energies is due to the 
unlike charges of the positron and antiproton, and does not occur for the like charges 
case of antiproton-electron scattering. Hence this does not affect the conclusion that 
depolarization of an antiproton beam in an electron cooler is negligible. 

6.4.2 A critical value for the target polarization 

The beam polarization will also decrease for low values of the target polarization. 
In fact there is a critical value of the target polarization Vt which keeps the beam 
polarization constant. If the target polarization is above this critical value the polar- 
ization of the beam will increase, and if the target polarization is below this critical 
value the beam polarization will decrease. The critical value is obtained by solving 

L d V - tanh(L d nz/r) [L in V + V T {A^ - K- m )] , 10 , 

v v> = t — n — T77 m t> t> a — i = ( 6 - 4 - 12 ) 

L d + tanh(L d nz/r) [L ia - VoV T A out \ 
for Vt, where the time dependence will cancel leading to 

v critical _ 2VoL in 



1 Electron cooling is described in section [7TT1 
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For target polarizations below this critical value the maximum beam polarization 
occurs at time r = 0, and for target polarizations above this critical value the 
maximum beam polarization occurs at large times r — > oo. 



6.5 Particles fed in for a limited time 



The Heaviside step function could be used in the system of equations to explain the 
case of particles input into a beam for a certain amount of time after which the 
input is turned off and no more particles are fed into the beam, but spin filtering 
continues f|. This scenario is under consideration by the VAX Collaboration [23J . but 
there has been no theoretical treatment of it to date. The Heaviside function is a 
piecewise continuous function which is zero in one region and one everywhere else, 
it is used in many mathematical modeling problems to describe an external effect 
turned on or off after a certain duration of time. In our case this external effect is 
the input of particles into the beam. The Heaviside function is defined as 

if r < r c 

H(r - t c ) = { (6.5.1) 

if t > t c 

and is used to describe an external effect turned on at time r c , but we require an 
external effect on initially and turned off at time r c , thus we need 

if < r < r c 

H(t) H(t -,.) = { (6.5.2) 

if r > t c 

Note that in the special case when r c = 0, [ H (r) — H (r) ] = 0, as this is in the 
second region. This describes a physical situation where particles are being fed in 
for zero seconds, which is the same as saying no particles are fed in, so the second 
order ODE should be homogeneous in this case, which it is. The extra term to add 
to the d N / d r equation to account for particles being fed in at a constant rate per 
second for r c seconds after which the input is switched off is /3 [H (r) — H ( r — r c ) ] . 

The equations are now broken into two pieces, i.e. discontinuous, but piecewise 
continuous. The solutions will be broken into two regions < r < r c and r > r c , 
where the solutions in the region < r < r c should equal those in the inhomogeneous 



2 This section makes extensive use of Laplace transform methods in solving differential equations. 
All results needed are presented in Appendix [F] 



99 




case presented in section 16.11 The initial conditions will be N(0) = and J(0) = 0, 
thus J'(0) = 0. The Heaviside function above is included in the second order ODE 
for J(r) to obtain: 

d 2 J(r) dJ(r) , . 

1 9 - (Ai + A 2 ) — ; + AiA 2 J(r) (6.5.3) 

dr dr 

= (A all -if in ) /3 [F(r) - if(r - r c )] , 

which can be solved by Laplace Transform methods, as described in Appendix [Fj to 
obtain 

r ) + A 2 (e AlT - 1) if < r < r c 

(6.5.4) 

Ai (e A2 ( r ~ Tc ) - e A2T ) + A 2 (e AlT - e Al(r ~ Tc) ) if r > r c 
where for convenience we have defined the constant factor 

Cj s W^-^) , (6.5.5) 
z Ai A 2 

One sees from eq. (I6.5.4P that J(t) = for all r when r c = 0, which is physically 
reasonable as there are never any particles in the beam if r c = 0. Also the Vt 
factor in C\ indicates that J(r) will always be zero if Vt = (i.e. if the target is 
unpolarized). It can also be seen that the complete solution in the region r > r c 
is the combination of the solution in the region < r < r c and an additional part 
dependent on r c ; which is 

d [Ai ( e x ^ T ~ T ^ - 1) + A 2 (1 - e Xl( - T ~ T ^)] , 

and immediately one sees that when r = r c this additional part vanishes. So when 
t = t c , i.e. at the boundary between the two regions, the two solutions match. 
Therefore the solution for J(r) is continuous as expected. The expression for J(r) in 
the first region < r < r c of eq. (16.5.41) is equal to the solution of the inhomogeneous 
system presented in eq. (16.1.111) . 

A similar analysis as that done for J(r) reveals the second order ODE for N(t) 

d 2 N(r) . ^ xN diV(r) x x Ar/ . 
, I - (Ai + A 2 — ^ + Ai A 2 N(t ) 
dr dr 

= nu (I M -D- m ) (3 [H(r) - H (r - r c )\ . (6.5.6) 
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The initial conditions are N(0) = and iV'(O) = (5 [1- H(-t c )}, the latter of 
which deserves comment. The rate N'(0) should be /3 when r c > and when 
r c = 0, corresponding to the case of particles being fed in for zero seconds. Note that 
H ( — t c ) = 1 when r c = and H ( — t c ) = when r c > 0, and note physically that 
t c , the duration for which particles are fed into the beam, cannot be negative. The 
Heaviside function in this initial condition forces the solution for N(t) to be split 
into three regions. On solving by Laplace Transform methods one obtains 



Ai 



1 



all 



+ A 2 



, A2 T 



/all - 



if < T < T r 



Ai 



, A 2 ( T - T C ) 



+ A 2 



Li 



all 



, A2 r 



g Ai t _ g Ai ( T — To ) 



(6.5.7) 



if r > tv. > 



if r c = 

where again for convenience we have defined a constant factor 



Co = 



-P (lall-^i 

2 L(j Ai A2 



(6.5.8) 



Again one sees that the complete solution in the region r > r c > is the combination 
of the solution in the region < r < r c plus an additional part dependent on r c ; 
which is 

C 2 [Ai (e x ^ T - T ^ - 1) + A 2 (1 - e Xl( - T - T ^)] . 

Immediately we see that when r = r c this additional part vanishes, thus the solution 
for N(t) is continuous. The solution in the first region < r < r c is equal to the 
solution from our inhomogeneous case presented in eq. (16.1.121) . and it satisfies the 
initial condition N(0) = 0. 

We now present results for the polarization V(r) = J(r) / N(t) as a function of 
time, in both regions. The polarization is undefined when there are no particles in 
the beam, thus we need not treat the case r c = 0. As expected in the < r < r c 



101 



region V(t) equals the solution of our inhomogeneous case presented in eq. (I6.1.13p . 
and in the region r > r c > one finds 

V(t) = (6.5.9) 

Pr Mali --frm) [ Ai e A2-r (e~ A2 ^ - l) + A 2 e Al " (l - e~ A ^)] 

Aie A ^[(/ all - J D in )e- A ^ C -(L in + L d )] + A 2 e A i-[(L in -L d )- ( I all - £> in ) e - A i ^ ] " 

The approximate initial rate of polarization buildup and the maximum polarization 
achievable will both reside in the < r < r c region, and thus will be identical to 
those presented in section 16.11 This is because the maximum polarization achievable 
occurs when the input rate is never switched off, i.e. in the < r < r c region. 
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Chapter 7 
Numerical results 



"The whole point of physics is to work out a 
number, with decimal points etc. ! Otherwise 
you haven't done anything. " 

Richard Feynman 

As an application of the theoretical work presented throughout the thesis we now 
investigate a possible method to produce a high intensity polarized antiproton beam 
by spin filtering off an opposing polarized electron beam. It is also outlined how this 
work can be applied to polarizing antiprotons by spin filtering off a polarized hydro- 
gen target. Firstly a description of the electron cooling technique to refocus the beam 
after scattering off the target each revolution in order to maintain high beam density 
is presented in section 17.11 Then the various experimental input parameters, such as 
revolution frequency, target areal density, target polarization and the effective accep- 
tance angle, needed to obtain realistic numerical estimates from our mathematical 
formalism are each described in section 17.21 The benefits of using a lepton target 
are described in section I7.3[ before analyzing the case of spin filtering off an oppos- 
ing polarized electron beam. Finally spin filtering off a polarized hydrogen target is 
discussed in section I7.4[ in the three cases of hydrogen with only electrons polarized, 
hydrogen with only protons polarized and finally hydrogen with both electrons and 
protons polarized. In section 17.4.11 it is shown that electromagnetic effects dominate 
hadronic effects in pp scattering in the region of low momentum transfer of interest 
in spin filtering. 
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7.1 Beam cooling 



After interaction with the internal target in the storage ring many of the beam 
particles do not move exactly along the beam axis, i.e. they have acquired a small 
deflection angle. This causes the beam to spread transversely and this process is 
called beam "heating" , in analogy to the random motion of atoms in a hot thermo- 
dynamic gas. In order to maintain a well ordered beam, where all particles move as 
collinearly as possible and also to increase the beam transverse density, a method 
to counteract this beam spread is required. Fortunately such a method exists and 
has been utilized successfully in many experiments world wide over the past three 
decades. 

The method is called Electron Cooling and was invented in 1966 by G. I. Budker 
at the Institute for Nuclear Physics (INP) laboratory (later renamed the Budker 
Institute of Nuclear Physics (BINP) in his honour) in Novosibirsk |122j . as a way to 
increase the luminosity of pp and pp colliders. Electron cooling was first tested in 
1974 with 68 MeV protons at the NAP-M storage ring of INP Novosibirsk [T23] . and 
is currently operational at over a dozen storage rings world wide. 

The terminology of "cooling" is analogous to thermodynamic cooling, as the 
random transverse motion of the (anti)protons is dampened by cooling the beam. A 
"hot" beam has many particles with large transverse motion, whereas a "cool" (or 
"cold" ) beam has low transverse motion, i. e. all particles move collinear to the beam 
axis. 

An electron cooler is a device inserted into the storage ring, where the antiproton 
beam passes through a co-moving cold electron beam, and on multiple Coulomb 
scattering with the electrons the transverse motion of the antiprotons is reduced, i.e. 
the antiproton beam phase-space density is increased. One immediately thinks of an 
analogy with the temperatures of mixed gases: Gas A with a high temperature is 
mixed with gas B having a low temperature, after some time the combined mixture 
tends to a uniform temperature which is midway between the initial temperatures 
of the individual gases. If the gases could be separated afterwards, one could say 
that gas B has reduced the temperature of gas A. Fortunately the electrons can be 
injected into, and extracted out of, an antiproton beam easily by magnets which 
deflect charged particles at different angles depending on the mass of the charged 
particles. Given that the mass of the antiprotons is approximately 1800 times that of 
the electrons, the electrons can easily be completely removed from the mixture. The 
velocity of the electrons in a cooler is carefully set to equal the average velocity of the 
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antiprotons, to maximize the interaction time. The antiprotons undergo Coulomb 
scattering in the electron "gas" and lose transverse energy, which is transfered to the 
co-moving electrons until some thermal equilibrium is attained. The electrons get 
"heated up" but are discarded after each pass and new cold electrons are injected 
continuously. 

Electron cooling is conventionally used on low to medium energy (anti)protons. 
Many laboratories are now investigating high energy electron cooling [124J and it 
has been shown to work on 8.9 GeV antiprotons in Fermilab |125j . GSI Darmstadt 
are investigating electron cooling the antiproton beam in the HESR at up to about 
8 GeV It is expected that in the near future high energy electron cooling will 
be commonplace in many laboratories. 

In section [63]of this thesis we have proven that the depolarization of an antiproton 
beam due to electron cooling is negligible because of the low areal density of electrons 
in a cooler beam. Therefore the beam can be refocused after interaction with the 
target each revolution, without losing significant beam polarization. 



7.2 Input parameters 

As shown in the system of polarization evolution equations of Chapter to give 
the highest possible antiproton polarization after a given filter time, the maximum 
antiproton revolution frequency, maximum target areal density and maximum target 
polarization are required. We now investigate each of these in turn before computing 
the numerical quantities. 



7.2.1 Maximum revolution frequency 

An upper bound on the velocity of the antiprotons is the speed of light c = 3 x 
10 8 ms" 1 . The revolution frequency (u) is simply the reciprocal of the time taken 
for one revolution, which in turn is the circumference of the storage ring (L) divided 
by the velocity of the antiprotons (v = (3 c), i.e. 

„ = (7.2.1) 

Obviously one can maximize the revolution frequency by using a very high energy 
beam in an extremely small circumference storage ring, but this is limited by the 
power of the magnets to bend high energy antiprotons around such a small circum- 
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ference. A plausible example is treated in [75] with an antiproton polarizer ring of 
circumference L = 75 m and a beam velocity of /3 = v / c — 0.5 giving a revolu- 
tion frequency of 2 MHz. The VAX Collaboration proposes an Antiproton Polar- 
izer Ring (APR) of circumference L = 86.5 m with antiprotons of kinetic energy 
250 MeV which, using the relation p = \JT ( T + 2 M ), corresponds to a momen- 
tum of 729.13 MeV/c and a velocity of (3 = 0.6136, thus they propose a revolution 
frequency of v = 2.12662 MHz [23]. The High Energy Storage Ring (HESR) of the 
proposed Facility for Antiproton and Ion Research (FAIR) at GSI Darmstadt pro- 
vides a high energy example. Antiprotons with momentum 15 GeV/c will be stored 
in a 574 metre circumference ring, giving a velocity of /3 — 0.998 and a revolution 
frequency of 521628 Hz [23]. 

Having more than one target in the ring, or having the electron beam overlap 
with the antiproton beam at more than one point in the ring, has exactly the same 
effect as increasing the revolution frequency, i.e. using two targets has the same 
effect as doubling the revolution frequency, using three targets has the same effect as 
tripling the revolution frequency etc. This is obvious since we are using the revolution 
frequency as a measure of the number of times the beam passes through the target. 
Having R targets in the ring will increase the rate of polarization buildup by a factor 
of R. The ring is limited by space, and such targets would have to be purposefully 
built, so realistically R = 1 is most likely for VAX, anything above R = 5 would be 
very challenging. 

The effect of multiple targets in the storage ring could be included in our system 
of polarization evolution equations simply by multiplying the entire coefficient matrix 
by a parameter R, where R is the number of targets in the ring. This would carry 
directly into each of the solutions presented in this thesis under the substitution 
nv —t- nv R. Note one must set R — 1 to compare to other work in the field which 
assume only one target in the ring. 

One could choose v — 5 MHz as a best case scenario available in the near future. If 
necessary two interaction regions of the opposing electron beam with the antiproton 
beam could be used to achieve this effective revolution frequency. 

7.2.2 Maximum target areal density 

As shown in Chapter [5] the rate of polarization buildup is highly dependent on the 
target areal density. One requires as high as possible a target areal density to achieve 
the highest rate of polarization buildup. In this section we review the maximum areal 
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densities of different types of target currently available. 

The FILTEX polarized hydrogen target, developed in the early 1990's, had an 
areal density of 6 x 10 13 atoms per cm 2 [54J. Since then some advances have been 
made in polarized atomic gas targets, in particular by the HERMES Collaboration. 
The HERMES Collaboration has produced and used polarized hydrogen and deu- 
terium targets with densities of up to 10 14 atoms per cm 2 [5T]. It is expected that 
this maximum areal density of polarized atomic gas targets could be increased by 
a factor of 100 in the near future. For a recent review of polarized gas targets see 
Ref. [129j . The PANDA Collaboration, also at the Facility for Antiproton and Ion 
Research (FAIR) at GSI Darmstadt, aim to produce a hydrogen pellet target of 
areal density up to 10 16 atoms per cm 2 [126J, similar to the target in operation at 
the WASA detector in COSY Julich. 

The maximum areal densities of polarized electron beams is many orders of mag- 
nitude lower than that of atomic targets, because of the electromagnetic repulsion 
felt by the like-charge electrons in the beam. This effect is absent for atoms which 
are electrically neutral. Typical areal densities of polarized electron beams pro- 
duced thus far is 10 s electrons per cm 2 . Given the enormous research and devel- 
opment effort that is currently been afforded to the International Linear Collider 
(ILC) project, which will use electron and positron beams, the maximum electron 
beam areal densities can be expected to be increased in the near future |130[ 1131] . 
Perhaps polarized electron beams with areal densities of up to 10 9 electrons per cm 2 , 
or even 10 10 electrons per cm 2 may be available in the next decade. It has been 
claimed that electrons stored in a Penning trap may soon reach areal densities of 
10 12 electrons per cm 2 [66J. 

For a realistic maximum value of the polarized electron beam areal density avail- 
able in the near future we use 10 12 electrons per cm 2 in our numerical calculations. 

7.2.3 Maximum target polarization 

Beams of electrons and positrons with polarizations of up to Vt = 0.9 have been 
produced and utilized in many laboratories |127| 1128] . Polarized internal targets of 
atomic hydrogen and deuterium with polarizations of nuclei, electrons or both of 
up to Vt = 0.9 have also been constructed, in particular the HERMES polarized 
hydrogen and deuterium targets which have been operated both longitudinally and 
transversely polarized [5TJ 1129] . There is an uncertainty of about 3% to 5% in the 
measurement of these polarizations, thus polarizations of above 0.95 are impractical 
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to produce. As a target polarization that should be state of the art in the coming 
decade before VAX will be realized we pick a target polarization of Vt = 0.9 in the 
following numerical calculations. 



7.2.4 Effective acceptance angle 

The acceptance angle 8 acc introduced earlier, is an idealistic simplification. It assumes 
all beam particles are moving exactly along the infinitely narrow beam axis, hence 
neglecting the following two effects: 

(a) In reality the beam has a finite extent, i.e. a r.m.s. radius of about 8 mm at the 
target [23]. So some particles are moving collinearly with the beam but at a distance 
of a few millimetres from the beam axis. 

(b) Not all particles move collinearly, many particles have a slight angle of motion 
with respect to the beam axis due to scattering off the target. 

Some of these particles that are scattered at or even slightly less than the idealistic 
acceptance angle will be lost. Hence there is an effective ring acceptance angle 
^acc CCtlvc "which is always less than the idealistic ring acceptance angle , i.e. 

^acc ectlvc < ^acc lve • The effect of this is to lessen the region of integration for the 
"in" spin observables thus reducing them, and increase the region of integration for 
the "out" spin observables thus increasing them. This effect also reduces the beam 
lifetime. 

These effects have been investigated by the VAX Collaboration for the polarized 
proton beam scattering off an internal hydrogen target at COSY Jiilich by com- 
paring the theoretically calculated idealistic loss cross-section to the experimentally 
measured beam lifetime |132j ; the results were as follows. 

They calculate the ratio of the effective spin-averaged loss cross-section c^, tive to 
the naive spin-averaged loss cross-section cr°^ ve , the former of which is proportional 
t° (0acc CCtlvc )~ 2 an d the latter to (O^c" ) 2 ^ because scattering is dominated by the 
Coulomb interaction. One has 



out f n naivc\ 2 

'effect ive _ v a< c ' /y 2 2^ 



U out ( a effective ^ 2 

naive V acc / 



therefore the ratio between the effective and naive acceptance angles is 

n effective i 

T . = -=. (7.2.3) 

(j naive /n 
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For the COSY polarized proton beam scattering off an internal hydrogen target it 
was found that F = 1.14 |132j . i.e. l/y/F = 0.94 and one has 

C ective = i C ve = 0.94 Cf e - (7-2.4) 
v b 



In the numerical calculations that follow we shall take these results and use an 
effective acceptance angle that is 0.94 times the idealistic acceptance angle of the 
storage ring. 



7.2.5 Minimum scattering angle 

From quantum mechanics one has the relation, |q| b = h, between the modulus of the 
three- momentum transfer (|q|) in an interaction and the impact parameter (b). Since 
h is set to one throughout this thesis one can relate the modulus of the minimum 
three-momentum transfer of an interaction to the maximum impact parameter by: 

1 

llminl — 7 j (7.2.5) 

"max 

and using —t = |q| 2 one has that the minimum squared momentum transfer is 

- t min = J- . (7.2.6) 

"max 

This minimum value for 1 1 1 ensures that there is no singularity from the \jt depen- 
dence of many of the spin observables presented in Chapter HI 
There are two cases of interest to us: 



1) A polarized hydrogen target 



For the case of a polarized hydrogen target the maximum impact parameter is given 
by the Bohr radius ag. For an antiproton scattering off a hydrogen atom at impact 
parameters greater a# the Coulomb fields of the atomic electron and proton screen 
each other, hence the antiproton sees the atom as electrically neutral and does not 
interact with it. Hence we have the relation 

- t min = — = a 2 m 2 = 0.000013912 (Me V/c) 2 , (7.2.7) 
where a is the fine structure constant and m the mass of the electron, and the relation 
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a B = (am) 1 has been used. 

One can convert this to the minimum laboratory frame scattering angle # m in using 
eq. (17.4. lUp to obtain §7\ 

#min = , (7.2.8) 

Plab a B 

where pi a b is the laboratory frame antiproton momentum. Using the relations p = 
\JT (T + 2 M ) and a B = (am) - , where T is the Laboratory frame kinetic energy 
of the antiprotons, one can rewrite this as 



am , 

(7.2.9) 



T (T + 2 M 



where M is the mass of the (anti)proton. For FILTEX kinetic energies of T = 
23 MeV one finds that 6 m i Q ~ 0.02 mrad, far below both the acceptance angle of 
any storage ring and the maximum angle antiprotons are scattered off stationary elec- 
trons, hence verifying eq. (17.4.131) . For scattering angle less than this, corresponding 
to impact parameters greater than the Bohr radius, the Coulomb fields of the atomic 
electron and proton of the hydrogen atom screen each other, hence antiprotons do 
not interact with the hydrogen atom. 

2) A polarized electron beam 

For an electron beam of areal density n particles per femtometre squared the average 
distance between electrons is 1/ \/n fm. If an antiproton passes exactly equidistant 
from two electrons in the beam it will feel no force as the Coulomb fields of the two 
electrons will cancel each other. The maximum impact parameter is one-half of the 
average electron separation, i.e. 6 max = 1/(2 s/n). Therefore 

- t min = J_ = 4n, (7.2.10) 



max 



where one must convert the areal density n into units of (MeV/c) 2 using the conver- 
sion factor 1 fm' 2 = 38937.9323 (MeV/c) 2 [T34] . 

Bv Taylor expanding the left hand side of eq. (17.3. 2p with respect to the ring 
framqj scattering angle 6 r one obtains the approximate relation, valid for small 9 r 
and small 1 1 1 



(7.2.11) 



Pi 



1 The ring frame is described in detail in section 17.3.31 
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where Pp is the antiproton momentum in the ring frame. Therefore the minimum 
ring frame scattering angle in this case is 

n r ^ V ^min 
min ~ ~~~r 
fp 

Results for the minimum squared momentum transfer and minimum ring frame scat- 
tering angle for various values of the electron beam areal density are presented in 
Table ED 



n [cm 2 ] 


n [fm 2 ] 


t min [(MeV/c) 2 ] 


#min [ mr ad] 


10 8 


10 -18 


- 1.5575 x IO -13 


2.6 x KT 11 


10 9 


io- 17 


- 1.5575 x HT 12 


8.3 x I0~ u 


10 10 


10 -16 


- 1.5575 x HT 11 


2.6 x 1(T 10 


10 11 


IO" 15 


- 1.5575 x 1(T 10 


8.3 x 1(T 10 


10 12 


io- 14 


- 1.5575 x IO" 9 


2.6 x IO" 9 



'An 



Pi 



(7.2.12) 



Table 7.1: Results for the minimum squared momentum transfer t m [ n and minimum ring 
frame scattering angle 0^ in for various values of the electron beam areal density n. The ring 
frame antiproton momentum is fixed at pi = 15 GeV/c. The conversion factor 1 cm -2 = 
10~ 26 /ro~ 2 has been used. 



7.3 Spin filtering off a polarized electron beam 

7.3.1 The advantages of using a lepton target 

There are many advantages of using a polarized lepton target (or beam) over a 
polarized internal atomic target (hydrogen or deuterium) for spin filtering: 

1. There is no loss of beam intensity due to annihilation of the antiprotons with 
protons as there is in the nuclear targets. 

2. The polarization observables for antiproton - electron scattering are calculable 
in perturbative QED (as presented in this thesis), whereas for an atomic target 
currently less known hadronic polarization observables contribute. 



Ill 



3. Residual gas does not build up over time in the storage ring as would happen 
if an atomic target was used. 

The first of these is by far the most important, and has caused many groups to 
investigate methods to polarize antiprotons by spin filtering off lepton beams and 
targets. We now investigate various scenarios of spin filtering off pure lepton beams 
and targets. 

7.3.2 Antiprotons scattering off stationary electrons 

The maximum laboratory frame scattering angle for antiprotons scattering off sta- 
tionary electrons is m / M = 0.54 mrad as shown in eq. f 1 7 . 3 . 5 j) and Figure. 17-11 (a). 
This is below the acceptance angle of any storage ring [231 ES] so that all scattering 
off atomic electrons will be within the ring. Therefore stationary electrons can only 
contribute to the polarization buildup of the antiproton beam by selective spin-flip. 
The Budker and Julich groups claim that spin-flip effects while scattering within the 
ring are small for p e~ — > p e~ scattering; hence stationary electrons, and very low 
energy electrons such as in an atomic target, are not effective in transferring polar- 
ization to an antiproton beam [TU [72J . We propose the use of an opposing polarized 
electron beam of sufficient energy to increase the scattering angles of the antiprotons 
beyond acceptance as seen in Figure I7-T1 (b)-(f ) . This is the subject of the remainder 
of the chapter. 

Note another solution to this problem would be to use a polarized muon target. 
Stationary muons, having much more mass than electrons (m M m 200 m e ), would pro- 
vide a maximum laboratory frame antiproton scattering angle of 6^ ax « m M / M = 
113 mrad. This is far beyond the ring acceptance angles under consideration at 
VAX, hence allows selective scattering out of the beam to contribute to polarization 
buildup. Another positive aspect of muons is that they can be produced automat- 
ically highly polarized, through the decays of charged pions. Charged pions decay 
into muons and (anti)neutrinos: ir + — > fi + and 7r~ — > . Since neutrino's 

have only one possible polarization state, angular momentum conservation forces 
the produced muons to be polarized. Polarized muon beams have been used in many 
experiments, among them the seminal EMC and SMC experiments which ushered 
in a new era of interest in spin physics and in particular in the spin structure of 
nucleons. But it is feared that the density of such polarized muon beams will be too 
low with today's technologies to provide a reliable method of polarizing antiprotons 
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Plot of scattering angle 6 r verses t 
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Plot of scattering angle 6 r verses t 
for 1 MeV/c electrons 
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Plot of scattering angle 6 r verses t 
for 3 MeV/c electrons 




-120 -80 -40 
(c) t (MeV/c) 2 



SE5 



Plot of scattering angle 6 r verses t 
for 5 MeV/c electrons 
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Plot of scattering angle 6 r verses t Plot of scattering angle 6 r verses t 
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Figure 7-1: (aj Ring frame antiproton scattering angle 6 r versus squared momentum trans- 
fer t for antiprotons of momentum 729 MeV/c scattering off electrons in an atomic target, 
where one confirms that the maximum antiproton scattering angle is 0.54 mrad. (b)-(f) 
With a colliding electron beam the maximum ring frame antiproton scattering angle in- 
creases as the electron beam momentum increases in a direction opposite to that of the 
antiproton beam. The plots show the relationship between 9 r and t given in eq. ( 7.3.1\ ). 
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Plot of scattering angle 9 r verses t 
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Plot of scattering angle 6 r verses t 
for 100 MeV/ c electrons 
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Plot of scattering angle 6 r verses t 



for 200 MeV/c electrons 
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Plot of scattering angle 6 r verses t 
for 300 MeV/ c electrons 
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Figure 7-2: .Ring frame antiproton scattering angle 9 r versus squared momentum transfer t 
for antiprotons of momentum 15 GeV/c scattering off an opposing electron beam of increas- 
ing energy. The plots become more skewed towards higher \ t \ for higher opposing electron 



beam momentum. The plots show the relationship between 9 r and t given in eq. \7.3.1\). 



At higher energies the plots of the relationship between 9 r and t become skewed 
towards higher 1 1 1 as shown in Figure 17-21 above for HESR energies, where the 
antiprotons have ring frame momentum 15 GeV/c. 



7.3.3 An opposing electron beam 

A crucial variable for spin filtering is the acceptance angle, as it defines which parti- 
cles are scattered out of the beam and which are scattered at small angles remaining 
in the beam. This angle is with respect to the beam axis, so for a stationary target 
in a storage ring the acceptance angle is the scattering angle in the LAB frame. We 
are now investigating the use of a colliding electron beam instead of a stationary 
target. Here the acceptance angle does not correspond to either the scattering angle 
in the Centre-of-Mass frame or the LAB frame. We want to use a frame in which 
the antiproton scattering angle with respect to the beam axis still corresponds to 
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the angle featuring in our equations. Note the scattering angle of the electrons is 
irrelevant here, and its use will be avoided by conservation of four-momentum. A 
ring frame for elastic antiproton electron scattering 



p(Pi, M) + e-(P 2 ,m) ->• p(P 3 , M) + e~ ( P 4 , m) 
is defined as follows 2 ]: 



Pi 


= (#1,0,0, 


Pi) 


Ps = 


(E 3 , 


p 3 sm9 r , 0, 


p 3 cos 6 r ) 


P2 


= (£2,0,0,- 


-P2) 


Pa = 


(Ei — E 3 + E 2 , 


— p 3 sin8 r , 


Pi -P2-P3 cos6» r ) 



Looking at E 2 = p 2 + m 2 for the final state electron gives 

(E 1 + E 2 -E 3 ) 2 = p 2 sin 2 9 r + ( Pl - p 2 - p 3 cos 9 r ) 2 + m 2 



After substituting out the unknowns E 3 and p 3 this gives j an equation for the 
ring frame antiproton scattering angle 9 r in terms of the Mandelstam variable t = 
(P3 ~ Pi) 2 , 



t - 



cosO 1 



2M 2 + 



E 1 (t~2E?)p 2 - E 1 (2E 1 E 2 +t) Pl 



pi E 2 + P2E1 



2pi 



(2E 1 E 2 +t)p 1 + (2 E{ 



2{ Pl E 2 + P2E1) 



-t)p2 



(7.3.1) 



M< 



A Taylor expansion to 0(t) of the above provides some clarity of the behaviour at 
small \t\: 

t 



COS0 1 



1 + 



2p 2 



(7.3.2) 



Equation ( 17.3. ip . the ring frame analogy of the much simpler Centre-of-Mass frame 
relation eq. ( 14.5 .6p . is used to graph 9 r versus t in Figure IT^Tl There is no backward 
antiproton scattering, i.e. 9^ ax < ir/2, for opposing electron beam momentum of 



P2 < 



Mp x 



(7.3.3) 



E l + Pl + M ' 

using the p 2 ~ E 2 approximation to make the inequality strict. In this region the 



2 If one of the particles is initially at rest the ring frame equals the laboratory frame. 

3 Many of the expressions in this chapter can be re-written in terms of the relativistic velocity 
Pi = Pi /Ei of each particle, where i € {1, 2}. In particular p% E% + P2 E\ — E\ E2 + P2) and 
E 1 E 2 -p lP2 = E x £ 2 (l-/3i #2). 
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maximum ring frame antiproton scattering angle is given, via eq. ( 17.3. ip . by [7T] 

Pi E 2 + p 2 E 1 



sin # max 



M( Pl - p 2 ) 

which, using sin 0^6 for small angles, limits to the correct formula 



(7.3.4) 



6 r « 

"mm 



sin^' 



rn 
M 



0.54 mrad . 



(7.3.5) 



for stationary electrons (i.e. p 2 = and hence E 2 = m). Notice this maximum 
scattering angle for stationary electrons is independent of the antiproton momentum. 
One can raise # max , and hence scatter more antiprotons out of the beam, simply by 
increasing the electron momentum p 2 . Raising it sufficiently beyond the acceptance 
angle will increase K out and the rate of buildup of polarization. 

Equation (17.3.41) can be used to derive a relation for the electron momentum 
needed to scatter antiprotons beyond the ring acceptance angle 8^ cc , hence allowing 
selective scattering out of the beam to contribute to polarization buildup. This 
happens for electron momentum p 2 > p out where on assuming 8[ cc is small and 
p 2 « E 2 , i.e. that the electron mass is small compared to its momentum, one obtains 



P 



out 



Ex + pi 



(7.3.6) 



f) r 

acc 


[mrad] = 


1 


2 


5 


10 


20 


50 


pOut 


[MeV/c] = 


0.468678 


0.937356 


2.34339 


4.68678 


9.37356 


23.4339 



Table 7.2: This table, from eq. ( 7.3.6), shows the opposing electron momentum needed to 
scatter antiprotons of momentum 15 GeV/c out of the ring for various ring acceptance 
angles. The upper-bound for the electron momentum, given this antiproton momentum, is 
found from eq. \1.3. 3 ) to be 454.478 MeV/c. Typical storage rings have acceptance angles 
in the range of 1 mrad to 50 mrad. 

Note that typical storage rings have acceptance angles in the range of 1 mrad to 
50 mrad. 

The Lorentz invariant A in the ring frame is 



A 



4k s 



4 (piE 2 + p 2 E ± 



(7.3.7) 
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and the Mandelstam s variable is 



s = M 2 + m 2 + 2E 1 E 2 + 2 Pl p 2 . (7.3.8) 

The maximum squared momentum transfer \t\, corresponding to total backward 
scattering, is given in the ring frame by 

/ - —4k 2 - - ~4(Pi^2 + P2E,) 2 

4 cm ~ s ~ M 2 + m 2 + 2E 1 E 2 + 2 Pl p 2 ' 1 } 

using the ring frame expressions of eqs. (17.3.71 and 17.3. 8p . 



The spin observables needed for spin filtering presented in section 14.51 can be 
converted into the ring frame by using eq. (I7.3.7P . Figure 17-31 shows how the regions 
of scattering angle defined by the acceptance angle 9[ cc can be converted into regions 
of squared momentum transfer t. Hence the angular regions of integration of the spin 
observables, presented in Table 15.11 can be presented as regions of integration over 
squared momentum transfer t in Table 17.31 Using the expressions of section 14.5.21 
and Table 17.31 one has that for longitudinals polarization: 



8 7r 2 a 2 u„ , n , , 9 N [ ll 1 , 

A J H 1 



.2„2 „ rti 

*3 



s — m 



A v ' ti 



M 2 )ln-. (7.3.10) 



- I6n 2 M 2 a 2 (s + m 2 - M 2 f ( 1 f t3 1 

IQtt 2 M 2 a 2 (s + m 2 — M 2 ) 2 ( , t x , t 4 
= ^— In— + In — 

A V imin ^3 

16vr 2 M 2 « 2 (a + m 2 - M 2 f h t 4 
= T2 ln 7 T ' ( 7 - 3 - U ) 



4 Since for antiproton-electron scattering the longitudinal spin-transfer observable from sec- 
tion 14.5.21 is greater than the transverse spin-transfer observable we do the numerical calculation 
for the longitudinal case. It has been shown at RHIC that one can rotate the polarization of the 
beam from longitudinal to transverse, or vice versa, without any loss of polarization. Since the 
stable spin direction in a storage ring is transverse, it is likely that the beam will circulate in the 
ring with transverse polarization but be rotated to longitudinal directly before the target and back 
to transverse directly after the target. This is how RHIC operates [1331. The antiproton beam will 
eventually be utilized in a transversely polarized state in order to measure Att and hence obtain 
information on the transversity distribution of quarks in the nucleon. 
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9 




in in 
out 

Figure 7-3: Plot of scattering angle 6 versus squared momentum transfer t, in the region 
of opposing electron beam momentum defined by eq. \7.3.ty . The acceptance angle # acc is 
plotted as a horizontal dashed line showing the region where particles are scattered "out" 
of the ring (t\ to t%) and the regions where particles are scattered at small angles remaining 
"in" the ring (t m i n to t\ and t% to t&). Note t m \ n corresponds to the minimum scatter- 
ing angle mm , scattering below which is prevented by Coulomb screening. The squared 
momentum transfer £2 corresponds to the maximum scattering angle $max • 



Transverse polarization requires 



Longitudinal polarization requires 



,' U da , 
/,,„! = 2 7T I — (1/ 



.4 



out 



.4 



all 



27T 



2tt 



/: / v4 xx + A YY \ da 



t 3 



— dt 
dt 



2tt 



Axx + A Y y\ da 



— dt 
dt 



K X x + Kyy\ da 
dt 



dt 



+ 2./" ( K ™ ± M ^ dt 



2tt 



2vr 



2 J dt 
* min f Dxx + D YY \ da 



— dt 
dt 



t3 'D xx + D YY \ da 



— dt 
dt 



I oat = 2ir [ ^ dt 
Jts dt 

f' 1 „ da , 
A ont = 2tt / A zz — dt 



t?, 



dt 



' * min d a 
A M = 2tt / A zz — dt 



dt 



' * min d a 
K m = 2tt / K zz — dt 



ti 



dt 



'* 3 da 
+ 2tt/ K zz — dt 



dt 



' tmin J. a 

D-, n = 2tt / D Z z -7- dt 

tl dt 

2 7T / "d zz 4? dt 



dt 



Table 7.3: The entries in the system of polarization evolution equations, for both longitu- 
dinal and transverse polarization, integrated with respect to squared momentum transfer t. 
Table [3771 shows similar expressions but integrated with respect to scattering angle 9. 
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Now Ld follows directly from the above expressions via eq. (15.3. 17p . The natural 
logarithms are expressed so that they are all positive, i.e. their arguments are greater 
than one, by using 




To present numerical results for the case of section 16.21 where particles are fed into 
the beam at such a rate that the beam intensity remains constant, we also need J ou t : 




(7.3.12) 



Since all antiproton-electron scattering is purely electromagnetic we now have all 
expressions needed to present the polarization buildup as a function of time from 
eq. (IQTTfi]) . 

One can obtain the values of t\ and t% for various acceptance angles by solving 
eq. (17.3. ip . Table [73 presents results for values of pi = 15 GeV/c and^2 = 50 MeV/c, 
where one sees from eq. (17.3.41) that 6^ ax = 107.248 mrad. The maximum elec- 
tron beam momentum for this value of antiproton momentum, in order for the 
antiproton scattering angle to be less than ir/2, is found from eq. (I7.3.3P to be 
p 2 = 454.478 MeV/c. 

Raising the opposing electron beam momentum has the same physical effect, that 
more antiprotons are scattered out of the beam, as lowering the acceptance angle. 
The acceptance angle is a fixed parameter of a storage ring whereas the momentum 
of the opposing electron beam can easily be altered. Therefore it makes sense to pick 
a typical ring acceptance angle and investigate what effect changing the opposing 
electron beam momentum has on the polarization buildup time. This is done in 
Table 17.51 which follows. 

As correctly emphasized in Ref. [75] the important practical parameters for a 
method to polarize antiprotons are: 

- the polarization buildup time, 

- the degree of polarization achieved after this time, 

- the number of antiprotons available after this time, and 

- the phase space of the polarized antiprotons. 

The latter is a measure of how focused the final polarized antiproton beam is, and 
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e Lc [mrad] 


ti [(MeV/cJ J 


t 3 [(MeV/cJ J 


in — 

ti 


m 

^min ^3 


1 


-224.988 


-2.32222 x 10 6 


9.24199 


25.6962 


2 


-899.809 


-2.32219 x 10 6 


7.85584 


27.0824 


5 


-5617.55 


-2.32195 x 10 6 


6.02427 


28.9140 


10 


-22381.3 


-2.32111 x 10 6 


4.64158 


30.2966 


20 


-88135.2 


- 2.31770 x 10 6 


3.26946 


31.6688 


50 


- 498035 


-2.29123 x 10 6 


1.52617 


33.4120 



Table 7.4: Values of t\ andt%, in units of (MeV/c) 2 , and the natural logarithms involving 
them appearing in the integrated spin observables, for various acceptance angles, in units 
of mrad, obtained by solving eq. \7.3.1 ). Results are for values of pi = 15 GeV/c and 



P2 = 50 MeV/c, where one sees from eq. \7.3.J$ that #max = 107.248 mrad and from 
eq. fnnfy that U = -2.32223 x 10 6 (MeV/c) 2 . Here we take the most optimistic i„ 



^min 



1.5575 x 10 9 (MeV/c) 2 corresponding to an electron beam areal density of 10 12 cm 2 as 
discussed in Table\7.1\ 



given that the beam can be focused by electron cooling we do not worry about 
this parameter here. In section 15.3.11 it has been shown that the figure of merit 
FOM (r) = V 2 (r) N (r) has a maximum at twice the beam lifetime. Of interest then 
is the polarization achieved after two beam lifetimes, i.e. the polarization achieved 
when the beam intensity has decreased by a factor of e 2 ~ 7.389. The later provides 
a combined measure of the second and third points listed above. 

Equations (15.3.141 and 15.3. 18|) can be used to provide an estimate of the polariza- 
tion achieved after two beam lifetimes: 

^optimum) = V(2n) « ~ 2 ^° Ut . (7.3.13) 

out 

It is interesting to notice that this optimum polarization is not dependent on the 
electron beam areal density n or the antiproton beam revolution frequency z/, al- 
though the time taken to achieve this polarization is strongly dependent on both 
of these parameters, as shown in eq. (15.3. 14p . In fact to this approximation, which 
is valid at times of twice the beam lifetime, the optimum polarization achieved is 
dependent only on energy and the ring acceptance angle, i. e. 9 acc or t\ and £3 . Using 
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the expressions for K out and J out presented in eqs. (17.3.101 and 17.3. 12[) one can write 



ln *3 

P(r optimum ) = V(2n) « - _^ M2 ~T~^ ± T > ( 7 - 3 - 14 ) 

is *i 

which we now maximize with respect to acceptance angle and energy. One finds 
that V (Toptimum) is maximal for large t\ ~ t 3 , i.e. for large acc ~ fl ma x, as shown 
in Figure 17-41 and also for high energies. Note these equalities cannot be strict 
otherwise V (T optimum ) is undefined, but this is physically reasonable as if 6* acc = # max 
there is no scattering out of the ring and hence the beam lifetime r* is undefined. 




Figure 7-4: A surface plot of the dependence on t\ and t% in V ( r op timum ) • The units of 
each axis are (MeV/c) 2 . 

Since V ( T opt i m um ) is maximal at high energies we perform the following numerical 
calculations at the HESR energies. Thus the antiproton beam will have ring frame 
momentum 15 GeV/c. Stored antiprotons, having already been accelerated to these 
energies, could be polarized directly in the HESR ring. This eliminates the need for a 
purpose built low energy Antiproton Polarizing Ring, providing another advantage of 
this method over spin filtering off a polarized atomic gas target. This will also avoid 
the problem of having to accelerate a polarized antiproton beam past depolarizing 
resonances in the storage ring. 
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Numerical results: Initial treatment 



Taking the discussion of section 17.2.41 into account we fix the effective acceptance 
angle at 9^ cc = 50 mrad, the highest acceptance angle under consideration by the 
VAX Collaboration [23] • We shall investigate spin filtering with HESR parameters 
of antiproton momentum of 15 GeV and revolution frequency of 521628 Hz. For this 
ring acceptance angle and antiproton momentum one finds from eqs. (I7.3.3l and l7.3.6p 
that the minimum and maximum opposing electron beam momentum, in order to 
provide scattering out of the beam but no backward scattering, are 23.4339 MeV/c 
and 454.478 MeV/c respectively. 



p 2 [MeV/c] 


tj [(MeV/c) 2 ] 


f s [(MeV/c) 2 ] 


U [(MeV/c) 2 ] 




In 

^min ^3 


50 


- 498035 


-2.291 x 10 6 


-2.322 x 10 6 


1.52617 


33.4120 


100 


- 520975 


-5.230 x 10 6 


-5.238 x 10 6 


2.30641 


33.4452 


200 


- 539229 


- 1.119 x 10 7 


- 1.119 x 10 7 


3.03262 


33.4782 


300 


- 546436 


- 1.718 x 10 7 


- 1.718 x 10 7 


3.44796 


33.4914 


400 


- 550267 


-2.317 x 10 7 


-2.317 x 10 7 


3.74039 


33.4983 


454 


-551672 


-2.641 x 10 7 


-2.641 x 10 7 


3.86869 


33.5009 



Table 7.5: Values of t\, £3 and £4, in units of (MeV/c) 2 , and the natural logarithms 
involving them appearing in the integrated spin observables, for various opposing electron 
beam momenta, in units of MeV/c, obtained by solving eq. {7.3.1 ). Results are for fixed 
values of 9[ cc = 50 mrad and pi = 15 GeV/c, where one sees from eq. (7.3.6) that p ont = 
23.4339 MeV/c. Again we take the most optimistic t min = -1.5575 x 10~ 9 {MeV/c) 2 
corresponding to an electron beam areal density of 10 12 cm 2 as discussed in Table \7.1\ 



Using the values^ of \n{t^/t\) and ln(ti t 4 /t min t 3 ) presented in Table 1731 one can 
obtain numerical values of K out , L in , L d and J out using eqs. f)7.3.10l to l7.3.12p . Hence 
using the results of section 15.31 one can obtain numerical values of the maximum 
polarization achievable for various opposing electron beam momenta, as presented 
in Table 17.71 The beam lifetime, r*, is the time taken for the beam intensity to 
decrease by a factor of e = 2.718. The figure of merit has a maximum at twice the 
beam lifetime, as described in section 15.3. 1\ hence we are principally interested in 
the polarization achieved after this time. 



5 At these high energies £3 m ti, due to the skewness of the relationship between 9 and t in the 
ring frame, as seen in Figure T7-2I Hence the major contribution to ln(ti i4/i m i n £3) comes from 
ln(£i/t m in). 
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p 2 [MeV/c] 


lout [mb] 


K out [mb] 


L in [mb] 


L d [mb] 


50 


0.00257272 


- 0.00232374 


0.01068 


0.0108828 


100 


0.00282946 


-0.0017559 


0.00267275 


0.00310499 


200 


0.00288982 


-0.00115439 


0.000668854 


0.00123563 


300 


0.00290078 


- 0.000874996 


0.000297386 


0.000841777 


400 


0.00290458 


-0.000711905 


0.000167314 


0.0006622 


454 


0.00290566 


- 0.000648743 


0.00012989 


0.000598142 



Table 7.6: The integrated spin observables, in units of millibarns, for various opposing 
electron beam momenta, in units of MeV/c. Results are for fixed values of'6[ cc = 50 mrad, 
pi = 15 GeV/c and electron beam polarization Vt = 0.9. 



p 2 [MeV/c] 


^max 


t* [seconds] 


V(2n) 


50 


0.096990 


7.45156 x 10 11 


0.096990 


100 


0.273516 


6.77541 x 10 11 


0.269873 


200 


0.545529 


6.63390 x 10 11 


0.424074 


300 


0.691294 


6.60883 x 10 11 


0.412933 


400 


0.772397 


6.60019 x 10 11 


0.372750 


454 


0.801983 


6.59772 x 10 11 


0.350905 



Table 7.7: The maximum polarization achievable, the beam lifetime and the optimum 
polarization achieved after two beam lifetimes for a stored antiproton beam, for various 
opposing electron beam momenta, in units of MeV/c. Results are for fixed values of9[ cc = 
50 mrad, p\ = 15 GeV/c, electron beam polarization Vt = 0-9, electron beam areal density 
n = 10 12 cmT 2 = 10~ 15 mb~ l and antiproton revolution frequency 521628 Hz. The values 
for the integrated spin observables are taken from Table \7.6] and inserted into eq. 115.3.18]) . 

The positive values for the induced antiproton polarization indicate that it is 
orientated in the same direction as the polarization of the electron beam. This is 
a result of the negative sign of K out , the positive signs of the other integrated spin 
observables presented in Table 17.61 and the overall minus sign of eqs. (15.3.161 and 
15.3. 18p . 

The maximum polarization achievable, satisfying the relation of eq. (15.3. 18j) . in- 
creases with increasing opposing electron beam momentum as a consequence of L m 
and Ld decreasing faster than i^out with increasing opposing electron beam momen- 
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turn. 

While the polarization achieved after two beam lifetimes in the above investiga- 
tions are very high one must note that the time taken to reach these polarizations is 
unpractically long. In practice one requires the beam lifetime to be of the order of a 
few hours, as opposed to the 10 11 seconds ~ 3171 years in the above cases! 

One notes that the closer t\ is to £3, i. e. the closer # acc is to 6 max , the lower the rate 
that particles are scattered out of the beam and hence the longer the beam lifetime. 
Since spin filtering must continue for two beam lifetimes to achieve the optimum 
polarization the beam lifetime in practice must not be greater than a few hours. 
Therefore one should maximize V ( 2 ) in eq. ( 17.3. 14p subject to the constraint 
that r* < 5 hours, which will allow for the optimum polarization to be achieved in 
less than 10 hours. This constraint can be presented in terms of £1 and £3 as 

2 , 1 

= — < 5 hours = 18000 seconds 



n v J out 9000 n v 

11 A . . 

.-. > ? > 0. 7.3.15 

£3 k ~ 72000 n 2 a 2 nv(s - m 2 - M 2 ) 2 

The minimum value for this difference is essentially energy independent, but highly 
dependent on the electron beam areal density n and the antiproton beam revolution 
frequency v. We shall present the numerics for two cases: (1) for n = 10 12 cm -2 
the best electron beam areal density that will be available now or in the near future, 
and (2) for an ideal case of n = 10 20 cm" 2 assuming great advances in electron beam 
areal densities. The latter simply provides a verification that the method of antipro- 
ton polarization buildup by spin filtering off an opposing polarized electron beam 
works in principle. One finds the minimum value for the difference of reciprocals 
of £, as presented in eq. (17.3.150 . to be 5.28527 x 10 6 (MeV/c)~ 2 for case (1) and 
0.0528527 (MeV/cT 2 for case (2). 

Unfortunately with today's electron beam areal densities and antiproton beam 
revolution frequencies the rate at which antiprotons are scattered out of the ring will 
be very low, even with the lowest possible ring acceptance angles of about 1 mrad. An 
extreme case provides a lower bound on the antiproton beam lifetime with today's 
technologies: 15 GeV/c antiprotons scattering off a 454 MeV/c opposing electron 
beam in a ring with lowest possible acceptance angle of 1 mrad and antiproton 
beam revolution frequency 521628 Hz, which still gives a beam lifetime of 2.63468 x 
10 8 seconds ~ 8.35 years. Hence with today's electron beam areal densities and 
antiproton beam revolution frequencies the beam lifetime will always be of the order 
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of years instead of the required hours. Thus the constraint presented in eq. (17.3. 1 5|) 
is too strict with current parameters, and one cannot maximize the polarization 
achieved after two beam lifetimes subject to it. 

We now investigate the scenario where unpolarized antiprotons are continuously 
fed into the beam at such a rate to cancel the rate that antiprotons are being scattered 
out of the beam, the beam intensity remaining constant, as presented in section I6T2"1 
The time taken to reach 8% polarization without any loss of beam intensity is pre- 
sented for various opposing electron beam momenta. Inserting the integrated spin 
observables presented in Table I7T61 into eq. (16.2.31) . where A a u — K in = K out in the 
pure electromagnetic case of interest here and J a u — D m = J out + 2L in , gives the 
results presented in Table 17.81 below. 



p 2 [MeV/c] 


Pmax 


r 8 % [seconds] 


50 


0.087385 


1.97920 x 10 11 


100 


0.193311 


1.25265 x 10 11 


200 


0.245759 


1.78586 x 10 11 


300 


0.225285 


2.40580 x 10 11 


400 


0.197800 


3.06730 x 10 11 


454 


0.184451 


3.44399 x 10 11 



Table 7.8: The maximum polarization achievable and the time taken, in units of seconds, 
to reach 8% polarization in the scenario where unpolarized particles are fed into the beam 
at the same rate particles are being scattered out of the beam, the beam intensity remaining 
constant N(t) = Nq, for various opposing electron beam momenta, in units of MeV/c. 
Results are for fixed values of 9^ cc = 50 mrad, p\ = 15 GeV/c, electron beam polarization 
Vt = 0.9, electron beam areal density n = 10 12 cm~ 2 = 10~ 15 mb^ 1 and antiproton 
revolution frequency 521628 Hz. The values for the integrated spin observables are taken 
from Table ["Of and inserted into eq. A6.2.3\) . 

Due to the long times taken to reach significant polarization it is apparent that 
neither of these methods of obtaining a high intensity polarized antiproton beam are 
practical at present. The key parameters limiting the rate of polarization buildup 
are the areal density of the opposing polarized electron beam, and the revolution 
frequency of the antiproton beam in the storage ring. The antiproton beam revo- 
lution frequency and the areal density of polarized electron beams would have to 
increase by many orders of magnitude in order for these methods to provide signifi- 
cant polarization in a few hours, which would be required. However, considering the 
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immense research and development that will take place in the near future on electron 
and positron beams at the International Linear Collider (ILC), advances in electron 
beam areal densities can be expected in the coming years |13Uj 1131] . 

Numerical results: An ideal case 

To conclude this treatment let us now investigate an ideal case, assuming fanciful 
values for the key parameters that are not currently achievable, but may be achievable 
in the future, to show that in principle the method works. Let us assume a ring 
frame acceptance angle of 50 mrad and a very high electron beam areal density of 
10 20 cm -2 which gives t m i n = — 1.5575 x 10 _1 (MeV/c) 2 . We shall investigate spin 
filtering with HESR parameters of antiproton momentum of 15 GeV and revolution 
frequency of 521628 Hz. For this ring acceptance angle and antiproton momentum 
one finds from eqs. (17.3.31 and 17.3. 6p that the minimum and maximum opposing 
electron beam momentum, in order to provide scattering out of the beam but no 
backward scattering, are 23.4339 MeV/c and 454.478 MeV/c respectively. Since 
tmm is much larger in this case compared to the previous treatment one expects 
L in to be much smaller in this treatment. As a consequence in this case both the 
maximum polarization achievable and the optimum polarization achieved after two 
beam lifetime should be larger, and the beam lifetime should be much shorter. A 
similar analysis to that presented above gives the results presented in the following 
tables. 



p 2 [MeV/c] 


h [(MeV/c) 2 ] 


f 8 [(MeV/c) 2 ] 


U [(MeV/c) 2 ] 


ti 


In 

^min ^3 


50 


- 498035 


-2.291 x 10 6 


-2.322 x 10 6 


1.52617 


14.9914 


100 


- 520975 


-5.230 x 10 6 


-5.238 x 10 6 


2.30641 


15.0246 


200 


- 539229 


- 1.119 x 10 7 


- 1.119 x 10 7 


3.03262 


15.0575 


300 


- 546436 


- 1.718 x 10 7 


- 1.718 x 10 7 


3.44796 


15.0707 


400 


- 550267 


-2.317 x 10 7 


-2.317 x 10 7 


3.74039 


15.0777 


454 


-551672 


-2.641 x 10 7 


-2.641 x 10 7 


3.86869 


15.0802 



Table 7.9: Values of t\, t% and t^, in units of (MeV/c) 2 , and the natural logarithms 
involving them appearing in the integrated spin observables, for various opposing electron 
beam momenta, in units of MeV/c, obtained by solving eq. \7.3.1\ ). Results are for idealistic 
fixed values of9[ cc = 50 mrad andp\ = 15 GeV/c. Here we take the very far fetched t m \ n = 
— 1.5575 x 10 _1 (MeV/c) 2 corresponding to an electron beam areal density of 10 20 cm 2 . 
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p 2 [MeV/c] 


hm [mb] 


K out [mb] 


L in [mb] 


L A [mb] 


50 


0.00257272 


-0.00232374 


0.00479191 


0.00522841 


100 


0.00282946 


-0.0017559 


0.00120068 


0.00198469 


200 


0.00288982 


-0.00115439 


0.000300831 


0.00108163 


300 


0.00290078 


- 0.000874996 


0.00013382 


0.000798785 


400 


0.00290458 


-0.000711905 


0.0000753084 


0.000645125 


454 


0.00290566 


- 0.000648743 


0.000058469 


0.000586789 



Table 7.10: The integrated spin observables, in units of millibarns, for various opposing 
electron beam momenta, in units of MeV/c. Results are for fixed values of 9'[ cc = 50 mrad, 
pi = 15 GeV/c and electron beam polarization Vt = 0.9. 



p 2 [MeV/c] 


V 

• max 


[seconds] 




50 


0.208713 


7452 


0.208648 


100 


0.496115 


6775 


0.459284 


200 


0.751524 


6634 


0.517906 


300 


0.844405 


6609 


0.455736 


400 


0.889346 


6600 


0.395042 


454 


0.904861 


6598 


0.367339 



Table 7.11: The maximum polarization achievable, the beam lifetime t* and the polariza- 
tion achieved after spin filtering for two beam lifetimes for various opposing electron beam 
momenta, in units of MeV/c. Results are for idealistic fixed values of 9[ cc = 50 mrad, 
pi = 15 GeV/c, electron beam polarization Vt = 0.9, electron beam areal density 
n = 10 20 cm~ 2 = 10~ 7 mb^ 1 and antiproton revolution frequency 521628 Hz. The values 
for the integrated spin observables are taken from Table \7.1ty and inserted into eq. K5. 3.16\) . 

One emphasizes the best case presented above as achieving an antiproton beam 
polarization of 51.8% after just 1.84 hours when the beam intensity has decreased by 
a factor of e 2 = 7.389, by spin filtering off an opposing electron beam of momentum 
200 MeV/c. 

Since the beam intensity has decreased significantly by the time the polarization 
reaches a high percentage, we now redo the analysis for the scenario where unpolar- 
ized particles are fed into the beam at the same rate particles are being scattered 
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out of the beam, the beam intensity remaining constant N(r) = N , as presented 
in section 16.21 The time taken to reach 15% without any loss of beam intensity is 
presented for various opposing electron beam momenta. Inserting the integrated spin 
observables presented in Table ETU] into eq. (16.2.3p . where A all — K m = K out in the 
pure electromagnetic case of interest here and 7 a n — D in = I out + 2L in , gives the 
results presented in Table 17.121 



p 2 [MeV/c] 


V 

• max 


Ti 5 % [seconds] 


50 


0.172036 


3241 


100 


0.302115 


2515 


200 


0.297568 


3851 


300 


0.248546 


5597 


400 


0.209713 


7882 


454 


0.193168 


9504 



Table 7.12: The maximum polarization achievable and the time taken, in units of seconds, 
to reach 15% polarization in the scenario where unpolarized particles are fed into the beam 
at the same rate particles are being scattered out of the beam, the beam intensity remaining 
constant N(t) = No, for various opposing electron beam momenta, in units of MeV/c. 
Results are for fixed values of 9^ cc = 50 mrad, p\ = 15 GeV/c, electron beam polarization 
Vt = 0.9, electron beam areal density n = 10 20 cm -2 = 10~ 7 mb^ 1 and antiproton revo- 
lution frequency 521628 Hz. The values for the integrated spin observables are taken from 
Table \7Aty and inserted into eq. h6.2. 3\) . 

Notice that in the best case presented in Table 17.121 for an opposing electron 
beam of momentum 100 MeV/c, the antiproton beam polarization builds up to 15% 
in only 2515 seconds ~ 42 minutes, while maintaining constant beam intensity. 

While this idealistic treatment might be far from today's technologies, requiring 
an increase of eight orders of magnitude in the product nv, it highlights that the 
method of polarizing an antiproton beam by spin filtering off an opposing polarized 
electron beam works very well in principle. 



128 



7.3.4 A co-moving lepton beam 

It has recently been proposed by Th. Walcher et al. [751 1136] that a polarized positron 
beam, with very low relative momentum to an antiproton beam, could transfer po- 
larization to the antiproton beam. This method of spin filtering is entirely based on 
selective spin flip while scattering within the ring, as a co-moving positron beam can- 
not scatter antiprotons out of the ring. Coulomb forces dominate the spin transfer 
observables at such low energies, and the antiproton-positron interaction is chosen 
because of the attraction of unlike charges. A positron beam can easily be polarized 
by the Sokolov-Ternov effect as described in sections (12.3.31 and 15. ip . The basis for 
this proposal is a dramatic enhancement of the polarization transfer cross-section in 
the reaction pe + ^p^e + as calculated in Ref. [BE]- Such a dramatic enhancement, 
of over nine orders of magnitude, has been called into question by many groups, who 
point out that, if true, multiple scattering effects should not be neglected. In par- 
ticular Ref. |1U7] claims that there is an enhancement at low energies but by many 
orders of magnitude less than that claimed in Ref. [88J, and that the polarization 
transfer cross-section is still far too low to make the Walcher et al. proposal practi- 
cal at present. It also remains to be seen if the depolarization observable also gets 
enhanced greatly at such low relative velocities. An experiment has been proposed 
to test these claims [119] . As encountered earlier the relatively low areal densities of 
polarized positron beams is a crucial disadvantage of this method. Any advances in 
the technology of obtaining high intensity positron beams would greatly benefit this 
proposal. 

An interesting application of a great increase in the areal densities of polarized 
electron beams is that an electron cooler could utilize a high intensity polarized 
electron beam. This would allow for an antiproton beam (or a beam of any other 
particles for that matter) to be both cooled and polarized after repeated interaction 
with the electron cooler. This would represent a major advance in accelerator physics 
and hadron storage rings. 

7.4 Spin filtering off a polarized hydrogen target 

7.4.1 Electromagnetic and hadronic scattering 

The critical squared momentum transfer t c , below which electromagnetic effects 
dominate and above which hadronic effects dominate in pp scattering, is now derived. 
For both hadronic and electromagnetic scattering, at low \t\, the non-spin-flip 
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amplitudes (f>\ + 3 = + dominate the spin-averaged differential cross-section. The 
leading t imaginary part of the hadronic amplitude is given by the Optical Theorem: 



lm{<f> h + (s,t)} 



2 k cm ^/s (j. 



pp 

tot 



t=0 



i7T 



(7.4.1) 



and the leading t part of the pp electromagnetic amplitude is given by 



a (s - 2M 2 ) F 2 (t) 



. 8 i 



(7.4.2) 



where F\ w 1 for small 1 1 | and e s 1 ~ 1 as 5 is small |137] . The Coulomb phase 
shift, e* 51 , accounts for the small correction to the single-photon exchange amplitude 
coming from multi-photon exchange |138j . The hadronic amplitude is 



2 k cm y/s <jf p 



. 7T 



'■" t (i + p) e 6 *. 



(7.4.3) 



where p = p(s,t) = TZe^^} / Im{<f)+} the ratio of real to imaginary parts of the 
hadronic non-flip amplitude. For small 1 1 1 one has that e bt ~ 1 1105] . Thus to 
leading order in small 1 1 | the spin- averaged pp cross-section is 



da 
~dt 



oc 



a (s - 2M" 
i 



2 fccm a/s cr^' 



tot 



8tt 



i + p\ 



(7.4.4) 



The electromagnetic and hadronic amplitudes are of equal size when 
-a (s - 2M 2 



i.e. when 



tr 



\TTCt (s - 2M' 



pp 



2 fccm \/g gtoi 

8tt 



z + p. 



(7.4.5) 



8tt(x E] 



lab 



2 /Ccm \/s CTtof V 1 + P 2 



Plab O"tot 



(7.4.6) 



where the relations k cm y/s = pi a b M and = (s — 2 M 2 ) / (2M) have been used. 
Finally using the relativistic laboratory velocity (3^ = pi a b / £a a b gives 1105] 



■ 7r a 



(7.4.7) 



the critical squared momentum transfer below which electromagnetic effects domi- 
nate and above which hadronic effects dominate □. 
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At high energies /3i aD ~ 1 and p 2 ~ |139[ 1140] , giving the often used result 



■ re a 



to = --«=-. (7.4.8) 



"tot 

This region of momentum transfer is referred to in the literature as the Coulomb- 
Nuclear-Interference (CNI) region. 

In analogy to t c above, we now derive an expression for the critical antiproton 
laboratory frame scattering angle below which the electromagnetic interaction dom- 
inates the hadronic interaction in pp scattering. It is then shown that the scattering 
angles of importance in spin filtering are below this critical angle, and hence the 
electromagnetic pp cross-sections calculated in Chapters [3] and 0] provide a good 
approximation to the total pp interaction in this region. 

Using eqs. (14.5.51 and I4.5.6[) and the relation between the scattering angles in the 
Centre-of-Mass and LAB frames 

^ = »*, (7.4.9) 
sin 6>i ab k cm 

and Taylor expanding sin 2 9^ for small #i a b gives 

-t « p 2 ah 9 2 ah . (7.4.10) 

So the critical scattering angle below which the electromagnetic interaction domi- 
nates the hadronic interaction is 



*£b - (7-4.11) 

and inserting the expression for t c in eq. (14.2. ip one obtains that electromagnetic 
effects dominate hadronic effects for laboratory scattering angles 



6 < 6" w J ^ ; w 152mrad, (7.4.12) 

\]T(T + 2M)a t ^^ h ^TT7 

where M and T are the antiproton mass and kinetic energy in the Laboratory frame 
respectively [72], and the numerical result is for FILTEX kinetic energies of T = 
23 MeV where ag% w 325 mb H3UHD] and p » 0.1 [HO] PHI] . This gives the strong 
inequality [72] 

Omin < 0e < ^acc < ^ ab , (7.4.13) 
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where 6 e = 0.54 mrad is the maximum angle antiprotons are scattered by stationary 
(or atomic) electrons as shown in eq. (17.3. 5p and Figure 17-11 (a). Since spin filter- 
ing utilizes angles below 152 mrad, electromagnetic (QED) effects provide a good 
approximation to the total pp interaction in this region. 

In conventional particle physics experiments particles must be scattered out of 
the beam pipe into detectors for measurements to be made, hence no direct experi- 
mental observations can be made for scattering at angles below the ring acceptance 
angle. Thus there is very little experimental data on this region of very low an- 
gle scattering, which is of interest only in storage rings. In particular hadronic 
antiproton-proton amplitudes are completely unknown in this kinematical region. 
The LHC very forward detector TOTEM hopes to obtain some data on low angle 
proton-proton scattering in the near future |142j . 

One can derive an expression for the antiproton momentum in the laboratory 
frame that makes the squared momentum transfer for total backward scattering 
equal to the critical squared momentum transfer below which electromagnetic ef- 
fects dominate hadronic effects. Solving the equation t c = —4:k^ m = — 4M 2 p^ h / 's 
one obtains pi a b ~ 31.6 MeV/c, below which all scattering is electromagnetically 
dominated. 

7.4.2 Polarization states of a hydrogen target 

Unpolarized hydrogen atoms in a strong magnetic field equally populate each of four 
hyperfine states: 

I Tp Te ) I Tp Te ) I Tp Te ) | Tp Te ) 

It is explained in section l2.2.1l how these hyperfine states and pairs of hyperfine states 
can be isolated to give polarized hydrogen. We are particularly interested in three 
types of polarized hydrogen target, with all atoms in the hyperfine states as follows 

V p = 1 and V e = (7.4.14) 

V p = and V e = 1 (7.4.15) 

Itpte) = V p = 1 and V e = 1 (7.4.16) 

where we denote the polarization of the electrons in the hydrogen by V e and the 
polarization of the protons in the hydrogen by V p . In practice the atoms are not 



Tp Te ) Tp i-e ) 



Tp Te ) ' Tp Te ) 



132 



perfectly isolated in certain hyperfine states, thus the electron and proton polariza- 
tions in polarized hydrogen are less than one. The HERMES Collaboration have 
utilized polarized hydrogen targets with V e = 0.9 and/or V v = 0.9 [5T]. Spin fil- 
tering off a polarized hydrogen target in each of the polarization states presented in 
eqs. (17.4. 14147.4. 161) can be treated similarly to the treatment presented in section I7T31 
using the spin observables for antiproton-proton and antiproton-electron scattering 
presented in sections (14.5.11 and 14.5.21) respectively. 

In the second case where the electrons in the hydrogen target are polarized but the 
protons are unpolarized one has that a+ nt = <x° ut , i. e. particles in both spin states are 
scattered out of the beam at equal rates, since only the protons are massive enough 
to scatter the antiprotons beyond the ring acceptance angle. Thus while there is 
scattering out of the ring, it is not spin- dependent and does not lead to a buildup of 
beam polarization. Therefore only selective spin-flip in electromagnetic antiproton- 
electron elastic scattering can contribute to polarization buildup in this case, yet 
one also has the negative effects of antiprotons being scattered out of the beam and 
annihilating with the protons in the hydrogen target; decreasing the beam intensity. 
The only advantage this case has over using a pure lepton target, considering it 
has the disadvantages listed above, is that the areal densities of electrons in atomic 
targets is greater than those achievable in pure lepton targets to date. 

Some experimental tests must now be carried out to decide which of the three 
possible states of a polarized atomic target would be most effective in polarizing an 
antiproton beam by spin filtering |74j . The Antiproton Decelerator (AD) storage ring 
at CERN is the only source of antiprotons in the required energy range. Consequently 
spin filtering studies of antiprotons scattering off a polarized hydrogen target at the 
AD ring are planned in the near future [531 H43[I144] . In particular these experiments 
should determine which of the antiproton-proton spin-dependent cross-sections 

Aa[ p = a™ - a™ and Aa| p = a™ - a™ , (7.4.17) 

for longitudinal and transverse polarizations respectively is largest, and at what 
laboratory frame antiproton kinetic energies in the AD range 50 — 200 MeV they are 
maximal. This should provide essential information on the optimal parameters of a 
proposed Antiproton Polarizing Ring (APR) to produce a high intensity polarized 
antiproton beam at the FAIR facility at GSI Darmstadt. 

It is well known that there is a significant spin-dependence of the proton-proton 
total cross-section. This was the basis of the FILTEX experiment, and the re- 
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suiting polarization buildup of the proton beam was a consequence of this spin- 
dependence. At FILTEX laboratory frame kinetic energies of 23 MeV one has that 
| (<Xm P — Vff) — 122 mb [BB]- Any measured polarization buildup of the AD 
antiproton beam during the proposed spin filtering experiments will provide direct 
evidence of, and the polarization buildup rate will provide a measure of, a spin- 
dependence in the antiproton-proton total cross-section. 

Measurements of the spin-dependent antiproton-proton cross-sections at the AD 
ring will also provide the first experimental results to test and distinguish between the 
current models [73j H40[ 11451 11461 1147] of the hadronic antiproton-proton interaction 
in the non-perturbative regime. 

We have pointed out that a pure lepton target has many advantages over an 
atomic target, but has the disadvantage of lower target densities because of the 
electromagnetic repulsion of the leptons in a pure lepton target which is less severe 
in the electrically neutral atomic targets. Consequently it is now appropriate to 
outline some possible, albeit far fetched, solutions to this problem: 1) Positronium, 
an electrically neutral electron-positron bound state could in future be used as a 
high density polarized pure lepton target. Being electrically neutral it should allow 
for similar densities as atomic targets. 2) Similarly muonium, an electrically neutral 
electron-antimuon bound state could possibly be used as above. 3) A polarized muon 
target/beam would allow for spin-dependent scattering out of the ring, as discussed 
in section 17.3.21 Because muons have approximately 200 times the mass of electrons 
using a muon target will enhance the polarization transfer cross-sections Kxx and 
Kyy as seen in eq. (I4.3.3I) . 4) A polarized tau lepton target or beam would be 
even more favourable to our needs, because of the very large mass of the tau lepton 
(m T w 3477 m e « 17 m M ). We stress that these ideas are far fetched, and very far 
from today's technologies, in particular because of the very low mean lifetimes of the 
particles discussed. However they may be practical in future, and in the mean time 
they might prompt other solutions to the problem of low target areal densities. 
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Chapter 8 
Conclusions 



"I was born not knowing and have had only 
a little time to change that here and there. " 

Richard Feynman 

There has been much recent research into possible methods of polarizing an an- 
tiproton beam, the most promising being spin filtering, the theoretical understanding 
of which is currently incomplete. The method of polarization buildup by spin filtering 
requires many of the beam particles to remain within the beam after repeated inter- 
action with an internal target in a storage ring. Hence small scattering angles, where 
it is shown that electromagnetic effects dominate hadronic effects, are important. 

All spin-averaged and spin-dependent cross-sections and spin observables for elas- 
tic spin 1/2 - spin 1/2 scattering, for both point-like particles and non-point-like 
particles with internal structure defined by electromagnetic form factors, have been 
presented to first order in QED. Particular attention is paid to spin transfer and 
depolarization cross-sections in antiproton-proton, antiproton-electron and positron- 
electron scattering, in the low 1 1 | region of momentum transfer. Of the spin-averaged 
formula derived we highlight that a generalization of the Rosenbluth formula has been 
presented in a new compact Lorentz invariant form. It is a two-fold generalization in 
that the masses of both particles are included and both particles are taken to have 
internal structure determined by electromagnetic form factors. While these results 
are eventually applied to spin filtering later in the thesis they are not limited to this 
application. The complete set of spin 1/2 - spin 1/2 helicity amplitudes and spin 
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observables should prove useful to many other areas in particle physics. 

The complete set of spin - spin 1 electromagnetic helicity amplitudes have 
also been presented to first order in QED. These are useful in describing the spin- 
dependent scattering of deuterons off carbon nuclei for example. 

A thorough mathematical treatment of spin filtering has also been presented, 
identifying the two key physical processes involved: (a) selective scattering out of 
the ring and (b) selective spin flip while remaining in the ring. The dynamical 
properties of the physical system have also been highlighted and analyzed. Sets of 
differential equations which describe the buildup of polarization by spin filtering have 
been presented and solved in many different scenarios of interest. These scenarios 
are: 1) spin filtering of a stored beam, 2) spin filtering while the beam is being 
accumulated, i.e. unpolarized particles are continuously being fed into the beam at 
a constant rate, 3) unpolarized particles are continuously being fed into the beam 
at a linearly increasing rate, i.e. the particle input rate is ramped up, 4) the input 
rate is equal to the rate at which particles are being lost due to scattering beyond 
the ring acceptance angle, the beam intensity remaining constant, 5) increasing the 
initial polarization of a stored beam by spin filtering, 6) the input of particles into 
the beam is stopped after a certain amount of time, but spin filtering continues. 

The depolarization of a polarized beam on interaction with an unpolarized target 
or beam, as in the important case of electron cooling, has also been investigated and 
shown to be negligible. 

There are advantages of using a lepton target instead of an atomic gas target for 
spin filtering, principal amongst them that antiprotons will not annihilate with the 
target as they do with the protons in the atomic targets, leading to a loss of beam 
intensity. Since electrons in an atomic target are not massive enough to scatter an- 
tiprotons beyond the acceptance angle of any storage ring we have proposed using an 
opposing polarized electron beam, of momentum large enough to provide scattering 
of antiprotons beyond ring acceptance, as a possible method to polarize antiprotons 
by spin filtering. This is presented as a practical application of the theoretical work 
presented throughout the thesis. The areal density of the polarized electron beam is 
identified as the key parameter limiting the rate of antiproton polarization buildup 
in this proposal. After analyzing this proposal it is concluded that the areal densities 
of electron beams currently available would have to increase significantly in order for 
this method of polarizing an antiproton beam to be practical. 

While this thesis is devoted to investigating spin-dependent antiproton interac- 
tions and the theoretical background to spin filtering in light of a possible method 
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to polarize antiprotons, we emphasize that much of the work presented is applicable 
to many other areas of particle physics. 

Possible extensions to this work would include redoing the analysis for the case 
of antiprotons repeatedly interacting with a polarized deuterium target. Polarized 
deuterium targets have successfully been utilized by the HERMES experiment, and 
would be available for spin filtering studies in future. In the kinematical regime of in- 
terest in spin filtering antiproton-deuterium scattering consists of antiproton-electron 
and antiproton-deuteron scattering. The helicity amplitudes and spin observables 
for electromagnetic antiproton-electron scattering have been presented in this thesis, 
and those for antiproton (spin 1/2) - deuteron (spin 1) scattering could be derived 
in analogy to those for spin 1/2 - spin 1/2 and spin - spin 1 scattering presented 
here. 

We hope that the treatment of spin filtering presented in this thesis will clarify 
some of the confusion in the theoretical literature, and perhaps play some part in 
the eventual achievement of a high intensity polarized antiproton beam. Measure- 
ments obtained using such a beam should lead to a better understanding of the spin 
structure of the protons and neutrons. Given that the proton is the nucleus of the 
hydrogen atom, the most abundant element in the Universe, one cannot overstate 
the importance of a better understanding of its internal structure. 
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Appendix A 
Dirac algebra 



The 4x4 Dirac gamma matrices, in the Dirac-Pauli representation, are 



7 



I 

-I 



7 



I 

1 







for j G {1, 2, 3}, where 



0"! = 



1 

1 



<?2 



-i 

1 



0"3 = 



-aj 



1 
-1 



(A.l) 



(A.2) 



I = 



1 

















= 









1 











are the 2x2 Pauli spin matrices, with i = v— T, an d the 2x2 identity and zero 
matrices respectively. Feynman slash notation i> = 7 M p M , and the usual conven- 
tion of Greek characters representing four dimensional space-time indices {0, 1, 2, 3} 
and Latin characters representing three dimensional space indices {1,2,3} are used 
throughout the thesis. The Dirac equation is 

(i7 M ^ - m)tp{x) = 0, (A.3) 

where the plane wave solutions are 

ip(x) = u(p)e~ ip ' x and ip(x) = v(k) e +ik ' x , (A.4) 
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for particles and antiparticles respectively, where u are the spinors for particles and 
v are the anti-spinors for antiparticles. Applying the operator — ( i^d^ + m ) to the 
left of both sides of the Dirac equation yields the Klein-Gordon equation: 

+ m 2 )i){x) = 0. (A.5) 

The hermitian conjugates of spinors u — u(p, A) are as follows 





= m 7 ° 


(7°) f 


= 7° 




= 75 




= (7°) f - 




= 7° 7, 7° 


7° 


= 7o 


(7o) 2 


= 7°7o = 



7 ° M 



(A.6) 



14x4 • 

The completeness relations for the sum of spinors are 

^2u(p, A) u(p, A) = p + m, (A.7) 

A 

5^(MXM') = ft - m , (A.8) 

A' 

«(p, A) u(p, A) = 1 + m ) [l + 75 ${p, A)] , (A.9) 

u(M')v(M') = ±0 - m) [I + 7 5#(M')] , (A.10) 

where the polarization four- vector S tl (p, A) of a particle is orthogonal to its momen- 
tum four vector, i.e. S^ip, A)p M = 0, and is normalized such that S^S^ = —1. The 
spinors satisfy the Dirac equation as follows: 

(f-m)u(p) = 0, 

(t + m)v(k) = 0, (A.ll) 
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and conjugating these equations yields 

u(p)(f-m) = 0, 

v(k)(# + m) = 0. (A.12) 
Gamma matrix calculations can be simplified using: 

Yl» = 4, 

yy i, = -2Y, (A.i3) 

YY Yi» = ^Y p , 
YY-ffi* = -2YYY- 

Traces of products of gamma matrices are as follows 

Tr [ odd # of 7 matrices ] = , 

Tr[YY] = 4rf v , 

Tr [YYYY] = k(rf v T) pc - Tf p r) va + rf c r) vp ) , 

Tr[ 7 5 ] = Tr [YY] = Tr [ 7 " Y ' Y] = 0, (A.14) 

Tr[YYYYY] = -4:ie^ upa . 

The totally antisymmetric permutation tensor, also known as the Levi-Civita symbol, 
is defined as 

(+ 1 if pL v p a is an even permutation of 0123 
— 1 if pv pa is an odd permutation of 0123 (A. 15) 

otherwise 
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and it satisfies the following contraction identities: 



e^^e^xr = 6i6 p x 6Z + W6i + 6 v x 6£,6Z-6Z6Z6t-6Z6£,6l-6 v x 6t6Z, 

e^'e^xr = -2{5 P X 5° - 5^1), 

e^^e^.r = -65% (A.16) 

f ^ v P u f — -24 

where the Kronecker delta is defined, not using the Einstein summation convention, 

as 

{1 if pL — v 
(A.17) 
\S.\i±v 

and when using the Einstein summation convention is used to sum over repeated 
indices // £ {0,1,2,3} one has 5^ = 6$ + 6{ + 6% + 6l = 4. 
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Appendix B 

Relations between Mandelstam 
variables 




Figure B-l: The Mandelstam variables are defined from the general two particle to two 
particles scattering process, where p\ and P2 are the momenta of the incoming particles 
and pz and p^ are the momenta of the outgoing particles. 

The Mandelstam variables 179] are defined as follows 



s = ( Pl + p 2 ) 2 = ( p 3 + p A ) 2 , 

t = (P4-P2) 2 = (P1-P3) 2 , (B.l) 

u = {P3-P2) 2 = (Pi ~Pa) 2 , 
where we have used the law of conservation of four momentum 

Pi + P2 = P3 + Pi ■ (B.2) 
Initial and final state particles, which in elastic scattering are the same, are on-shell 
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thus Pi • Pi = P3 ■ P3 = M 2 and p 2 ■ P2 = Pa • Pi = m<2 , where M and m are the 
masses of the two particles in the elastic process. One can square the above relations 
to obtain 



s = m 2 + M 2 + 2 Pl -p 2 = m 2 + M 2 + 2p 3 -p 4 , 

t = 2m 2 -2p 2 -p 4 = 2M 2 -2p 1 -p 3 , (B.3) 

u = m 2 + M 2 -2p 2 -p 3 = m 2 + M 2 - 2p l ■ p 4 , 

hence one sees that 

Pi ■ P2 = P3- Pi and Pi-Pa = P2-P3, (B.4) 

and adding, using conservation of four momentum, gives the defining relation for 
Mandelstam variables 

s + t + u = 2m 2 + 2M 2 . (B.5) 
The above relation is the special case for elastic scattering of the general relation 

4 

s + t + u = J^m 2 . (B.6) 



Rearranging the t equation gives 

Pi-P3=(m 2 -0 and p 2 ■ p 4 = (^m 2 - . (B.7) 

Similarly for s and u 

Pi -Pa = P2-P3 = \ {m 2 + M 2 - u) , (B.8) 

Pi-P2 = P3-PA = ~(m 2 + M 2 - s) . (B.9) 

For convenience define 

R» = p»+p£ and r v = p 2 + p 4 , 
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thus we have 



R 2 = R^R^ = vl + 2pi -P3+P3 = 2M 2 + 2M 2 -t = 4M 2 -t, 



r = r r v 



(Am 2 -t) 



(B.10) 



Now 



R. r = R^ r ^ = (pf + pg) (p 2fl +P4 M ) = Pi -P2+P1 -P4+P2 -P3+P3-P4, 



= 2 



+ 2 



- (m 2 + M 2 - u) 



= s-u, (B.ll) 



similarly 



r • pi = r • p 3 = R-p 2 = R-Pa = -(s-u) 



(B.12) 
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Appendix C 

Derivation of the Gordon 
decomposition identities 

The Dirac gamma matrix structure of the most general proton electromagnetic cur- 
rent can greatly be simplified using the Gordon Decomposition Identity j82j 183] . 
which we now derive. 

The commutation and anticommutation relations 

= and = 2rT, (CI) 

where 

[Y,Y] = YY -tY and {Y,Y} = YY + YY, (C.2) 

and 7]^ y = diag(+l, —1, —1, —1) the Minkowski metric tensor, can be used to write 

[ 7 m 7 "j = YY - YY = - - i^Y) = 2j ti Y - ■ 

Hence 

ia"" = -(7" 7" - Y v ) = Y v ~ l"Y, 

but equivalently 

io^ = r]^ u - (2r^ - Yl") = Yl" - 77"". 
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Using the above relations one can compute 



u(p')ia^ u (p' u - p u ) u(p) 



u(p') \io^p' v - ia^p u ]u(p), 

u(p')[(Yi» - rnp'u - w -rr)Pu]u( P ), 

U{P')[1 U P' V 1" - P"*-P" + l"l U Pu]u{p), 

u(p')[|>V - {p + p'Y + 7"jfU(p), (C3) 



where the Feynman slash notation i> — 7 /i p jtl has been used. Now simplify by using 
the Dirac equation and its conjugate 



(f - M) u(p) = =► 
u(p') (/ - M) = =^ 

where M is the mass of the particle, to obtain 



f u{p) = M u(p) , 
u(p') f' = u{p') M . 



(C.4) 



u(p')ia» u (p' u - p v ) u{p) = u(p') [M 7 ^ - (p + p'Y + 7 "M]u(p), 

= w(p')[2M 7 ^ - (p + u(p). 

Rearranging gives the general form of the Gordon Decomposition identity: 

'(p + pT %a^{p'-p) u 



u(p') 7 M = ■u(p') 



2M 



+ 



2M 



u(p) . (C.5) 



Another Gordon decomposition identity, which we now derive, can be used to sim- 
plify the gamma matrix structure of the antiproton current using the v anti-spinors. 



Equation (1C.3j) for anti-spinors gives 



v(p')ia» v {p' v - p v ) v(p) = v(p')[f'-i^ - (p + p'T + iy]v(p). 



But now we must use the Dirac equation for antiparticles and its conjugate 



{f + M) v(p) = 
v(p') (/ + M) = 



f v ip) — ~ Mv(p) , 
v{p')f' = v{p') (-M) 



(C.6) 
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on the above equation to obtain 



v{p')i(T llv (p' v - p v ) v(p) = v(p') [-M 7 " - (p + p'Y + 7 M (-M)]v(p), 

= v{p')[-2M 1 tl - (p + p'Y] v(p) . 

Now rearranging gives the Gordon Decomposition identity for anti-spinors: 

(p + p'Y ia» u (p' - p) l 



v ip') 7 M v ip) — v (p') 



2M 



2M 



v(p) ■ (C7) 
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Appendix D 

Feynman rules for QED 



The Lagrangian for Quantum Electrodynamics (QED) is 

C = - m)V - \ (F^f - e^ 1 ^A lx . (D.l) 

Momentum is conserved at each vertex, fermion loops receive an additional factor of 
(— 1) and undetermined loop momenta are integrated over by: 




(D.2) 



We work in the Feynman gauge where the gauge parameter is set to £ = 1. 

The charge factor Q = —1 for an electron or antiproton, and Q — +1 for 
a positron or proton. Time increases from left to right in all Feynman diagrams 
throughout the thesis. 
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Photon propagator: 



QED vertex: 



External fermions: 



External antifermions: 



External photons: 



/ WW\y 




p 2 + it 

= u(p, A) (initial) 

= u(p, A) (final) 

= v(p, A) (initial) 

= v(p, A) (final) 

= e A*(p) (initial) 

= e«(p) (final) 
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Appendix E 

Sample Mathematica code 



The Mathematica code to derive the generic depolarization equation from eq. ( 13.2. lip 
is included in this appendix. The Mathematica add on package for High Energy 
Physics Tracer. m [87] is used, which defines G5 = 7 s and uses Spur [] for Trace. All 
other QED calculations in the thesis are done analogously. Only input commands 
are included. All parameters are as defined and used throughout the thesis. 

(* Import the file Tracer. m *) 

« Tracer/Tracer .m; 
VectorDimension [4] ; 
Spur[l] ; 

(* Spin projection operators *) 

SP1 = G[le, U] + G[le, G5]**G[le, SI]; 
SP3 = G[le, U] + G[le, G5]**G[le, S3]; 

(* Spin trace *) 

SpinTrace = G[le, pi + M U]**SP1**(H G[le, {nu}] + F (pl.jnu} 
+ p3.{nu}))**G[le, p3 + M U]**SP3**(H G[le, {mu}] 
+ F (pl.jmu} + p3.{mu})) ; 

SpinTrace2 = FullSimplify [SpinTrace/ . le -> 1]; 
(* Spinless trace *) 

SpinlessTrace = G[le2, p4 + m U]**(h G[le2, {mu}] + f (p2.{mu} 

+ p4.{mu}))**G[le2, p2 + m U]**(h G[le2,{nu}] + f (p2.{nu} + p4.{nu})); 
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NoSpur [1] ; 
Spur [2] ; 

SpinlessTrace2 = FullSimplif y [SpinlessTrace/ . Ie2 -> 2]; 
(* Contract all Lorentz indices *) 

msql = Simplify [ContractEpsGamma[SpinTrace2 SpinlessTrace2] ] ; 
msq2 = FullSimplif y [msql] ; 

(* Take the depolarization observable away from the spin-averaged *) 

SpinAveraged = Simplify [msq2/ . {SI -> 0, S3 -> 0}]; 
SpinTerml = SpinAveraged - msq2; 

(* Now make algebraic simplifications *) 

SpinTerm2 = Simplify [SpinTerml/. {p2.p4 -> m~2 - t/2, pl.p3 -> M~2 
- t/2, p4.p4 -> m~2, p3.p3 -> M~2, p2.p2 -> m~2, pl.pl -> M"2, 
p3.p4 -> s/2 - M"2/2 - nT2/2, pl.p2 -> s/2 - JT2/2 - nT2/2, 
p3.p2 -> M~2/2 + m~2/2 - u/2, pl.p4 -> M~2/2 + m"2/2 - u/2, 
Sl.pl -> 0, S3.p3 -> 0}] ; 

SpinTerm3 = Simplify [SpinTerm2/ . {p3 - pi -> q, p2 - p4 -> q, 
u -> 2m~2 + 2M~2 - s - t}] ; 

SpinTerm4 = Simplify[SpinTerm3/.{p3.Sl -> q.Sl, pi. S3 -> - q.S3}] ; 

SpinTerm5 = FullSimplif y [SpinTerm4 /.{F -> - F_2/2M, f -> - f _2/2m}] ; 

GenericDepolarisationEquation = FullSimplif y [SpinTerm5 /.{H -> F_l 
+ F_2, h -> f_l + f_2, pi. S3 - p3.S3 -> - q.S3, - pi. S3 + p3.S3 -> 
q.S3, p2.Sl - p4.Sl -> q.Sl}] 



The output is now the generic depolarization equation, for the electromagnetic scat- 
tering of two non-identical spin 1/2 particles with internal structure defined by form 
factors, used in the thesis. All helicity amplitudes and spin observables for elastic 
spin 1/2 - spin 1/2 scattering to first order in QED can be obtained from eq. (13.2. lip 
in a similar fashion. 
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Appendix F 

Laplace transform methods 



Laplace transform methods are used to solve the differential equations involving the 
Heaviside step function in section 16.51 All Laplace transform results needed are 
presented here. 

The Laplace transform £ { /(r) } = F(s) of a function /(r) is defined as 

P oo 

C{f(r)} = F(s) = / e-/(r)dr. (F.l) 

Jo 

It is a linear operator, which can be proven as follows: 



r oo 

£{ Cl A(r) + c 2 / 2 (r)} = / e- ST [a h(r) + c 2 / 2 (r)] dr, 

Jo 

P OO f oo 

= Cl / e~ ST /i(r)dr + c 2 / e -^/ 2 (r)dr. 
Jo Jo 



Ci£{/i(r)} + c 2 £{/ 2 (r)} 



where ci and c 2 are constants and /i(r) and / 2 (r) are arbitrary functions. 

Differential equations involving the Heaviside step function can be solved by 
taking the Laplace transform of the entire equation, applying the relations which 
follow, solving for £{ /(r) } and then taking the inverse Laplace transform of what 
is left to obtain the solution f(r) to the differential equation. 
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The Laplace transforms used in solving the differential equations in section 16.51 

are: 



C{f'(r)} = sC{f(r)} - /(0), (F.2) 

C{f"{r)} = s 2 C{f(r)} - sf(0) - /'(0) , (F.3) 

C{H(r-r c )} = 6 -^, (FA) 

C{H(r)} = - = £{!}, (F.5) 



where /(r) is an arbitrary function, ' denotes first derivative and " denotes second 
derivative. One obtains /(0) and /'(0) from the initial conditions of the differential 
equation. The Heaviside step function, H (r — r c ), is defined as 

f if t < T c 

H(t -t c ) = I (F.6) 
[ 1 if r>r c 

The Inverse Laplace Transforms go from right to left in the above list, for example 

/r X {^} = H(r-r c ) . (F.7) 
The Laplace inverse is defined, as always, such that 

CiC-'Vir)}} = f(r) = C-'iCifir)}}. (F.8) 
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